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lonen- und Atomquerschnitte 
im Plasma verschiedener Gase. 
Von 
H. MAECKER, TH. PETERS und H. ScHENK. 
Mit 13 Figuren im Text. 


(Eingegangen am 15. September 1954.) 


Die bekannten Zusammenhange zwischen der Temperatur und der Leitfahigkeit 

in einer Bogenentladung werden benutzt, um die Wirkungsquerschnitte der den 

Elektronenstrom bremsenden Ionen und Atome aus Stromstirke, Feldstarke, 

Elektronendichte- und Temperaturverteilungen verschiedener Bogenentladungen 

zu ermitteln. Die Ionenquerschnitte stimmen recht gut mit den von SPrITzER 

berechneten iiberein. Atomquerschnitte fiir die Elemente H, C, N, O und Fe werden 
mitgeteilt. 


I. Einfitthrung. 


Bei der Beschreibung der Leitfahigkeit eines Plasmas muB man bis 
auf die Wirkungsquerschnitte der Ionen und Atome gegeniiber Elek- 
tronenstoB zuriickgehen, weil diese die Zahl der ZusammenstéBe zwi- 
schen Elektronen und den iibrigen Plasmapartnern und damit die 
mittlere Driftgeschwindigkeit der Elektronen bestimmen. Die Ge- 
schwindigkeit bei der Feldstarke Eins, die Beweglichkeit, ist neben der 
Elektronenkonzentration die entscheidende Gr6Be fiir die Leitfahigkeit. 

Zur Messung der Wirkungsquerschnitte kann man sich hinreichend 
zylindersymmetrischer Bogenentladungen bedienen, fiir dieman zunachst 
aus dem Verhialtnis von Stromstarke zur Feldstarke den Leitwert je 
Langeneinheit ermittelt. Die gleiche GréBe sollte sich ergeben, wenn 
man mit Hilfe der gemessenen radialen Temperaturverteilung die spezi- 
fische Leitfahigkeit, in die die Wirkungsquerschnitte eingehen, berechnet 
und diese iiber die Querschnittsflache des Bogens integriert. So sind 
jedenfalls die Wirkungsquerschnitte der Ionen und Neutralteilchen ein- 
zusetzen. Damit ist die Bestimmungsmethode fiir die spezifische Leit- 
fahigkeit und fiir die gesuchten Wirkungsquerschnitte kurz skizziert. 

In dieser Abhandlung soll nun die bereits mitgeteilte Messung des 
Ionenquerschnittes am wasserstabilisierten Hochleistwngsbogen einer 
kritischen experimentellen Nachpriifung unterzogen und mit neuen theo- 
retischen Berechnungen von L. Spitzer und Mitarbeitern verglichen 
werden. Sodann wird iiber die Bestimmung des RAMSAUER-Querschnitts 
der neutralen H- und O-Atome im normalen wasserstabilisierten Bogen, 
wie er von G. JURGENS eingehend untersucht worden ist, berichtet. 
Diese MeBergebnisse werden gestiitzt durch Untersuchungen am 
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wasserstabilisierten Hochdruckbogen, der nach den Arbeiten von 
Tu. PETERS noch wesentlich geeigneter fiir diesen Zweck ist. Weiter 
werden die gleichen Methoden auf die von H. MAEcKER durchgemessene 
Sdule des Hochstromkohlebogens angewendet, um daraus die Wirkungs- 
querschnitte von C-, N- und O-Atomen zu gewinnen. SchlieBlich werden 
die von F. BurHorRN ausgefithrten Messungen am Fe-Pfund-Bogen heran- 
gezogen, um auch fiir die Fe-Atome die RAMSAUER-Querschnitte zu 
finden. 


II. Grundlagen der MeBmethode. 


Bei einer zylindersymmetrischen Bogenentladung kann man aus den 
gemessenen Werten fiir die Stromstarke J und die Feldstarke E den 
elektrischen Leitwert je Langeneinheit 2 nach dem OxnMschen Gesetz 
finden zu 


J 
rat. (1) 
Zum gleichen Ergebnis muB man kommen, wenn man die spezifische 
Leitfahigkeit o tiber den Bogenquerschnitt integriert: 
R 
X= 22 fa(r)rdr. (2) 
0 
Da nun die lokale Leitfahigkeit nicht direkt der Messung zuganglich ist, 
wohl aber die radiale Temperaturverteilung, muB tiber die gaskinetische 
Vorstellungen vom Ladungstransport der Elektronen im Plasma — der 
Anteil der Ionen ist wegen ihrer sehr geringen Beweglichkeit zu vernach- 
lassigen — ein Zusammenhang zwischen Temperatur und Leitfahigkeit 
gesucht werden. 

Auf das Elektronengas wirkt einmal die Kraft des elektrischen Feldes 
eE und zum anderen die durch ZusammenstéBe der Elektronen mit den 
ibrigen Plasmapartnern erzeugte Reibungskraft, welche als abgegebener 
Impuls in der Zeiteinheit definiert ist. Bewegt sich das Elektronengas 
mit der mittleren Driftgeschwindigkeit v, so gibt ein Elektron pro StoB 
den Impuls m,v ab, wobei als StoB alle die Vorgange gezihlt werden 
sollen, durch die die Impulskomponente in Feldrichtung im Mittel voll- 
stindig abgegeben wird. Finden in der Zeiteinheit y derartige StéBe 
statt, dann ist die Impulsabgabe pro Zeiteinheit und damit die Reibungs- 
kraft ym,v. Die Driftgeschwindigkeit v stellt sich nun so ein, daB sich 
antreibende Feldkraft und Reibung die Waage halten: 


gl See} 6) 


Die reziproke StoBzahl ist die StoBzeit rt, die wir einfach aus der ther- 
mischen Geschwindigkeit v,, und der freien Weglange 4, der Elektronen 


bestimmen: 
1 Vth 


oma fe A, (4) 


————— 


> 
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Die sich einstellende Driftgeschwindigkeit berechnet sich so zu 


(a Ae E=b,E, (5) 


2 = 
Me Vth 


worin der vor der Feldstairke stehende Faktor als Beweglichkeit b, be- 
zeichnet ist. Die strenge Ableitung dieser Beziehung nach der statisti- 
schen Mechanik fiigt noch einen Zahlenfaktor hinzu, der aber nahe bei 
Eins liegt, und den wir uns in das A, hineingezogen denken wollen. Aus 
der Elektronengeschwindigkeit folgt die Stromdichte durch Multipli- 
kation mit der Ladungsdichte des Elektronengases e+, zu 

p=en, bE, (6) 
was dem OxMschen Gesetz entspricht, wenn man als Leitfahigkeit 

ao=e-n,b, (7) 
definiert. Fiihrt man in (5) noch die gaskinetischen Beziehungen 


"(03 22 (8) 


me, 


und 
1 


A, — (9) 


My Qo +n, Q. 
ein, dann ist nach (7) die Leitfahigkeit auf die Temperatur T, auf die 
Elektronenkonzentration , und auf die beiden Wirkungsquerschnitte 


der Atome Q, und der Ionen Q, gegeniiber ElektronenstoB zuriick- 
gefihrt: 


ne 1 
o = ————*— - = 10 
\3m,kT No Qo +2, Q4 e) 


Hierbei ist 7, noch mit TJ durch die SAHA-Gleichung 


; ni - 24% 

NZ+1 come ae 2Uz+1— (2amk ) e RT (11) 
€ 7, 3 
nz bi h 


(Z = Ionenladungszahl, U = Zustandssumme, Ez = Ionisierungsener- 
gie) und durch die Forderung nach Quasineutralitat > Znz=n, Ver- 
bunden. Aus (1), (2) und (10) folgt dann 

R 


Say nate i ic ae : vy dr. (12) 


It == Pee $e et 
E [ 2 VeT 9% Qo + 249+ 
3m. jaz 0 ‘0 
Auf diese Gleichung werden wir uns bei den folgenden Auswertungen 
stiitzen. Im Bedarfsfall kann der Term 7 Q, bei Anwesenheit mehrerer 
Atomarten in ,Q,-+,Q,+--+: und der Term n,Q, bei Anwesenheit 
mehrfach geladener Ionen in 1, Q,+,.,Q.,+°°° "zs Q,, aufgeteilt 


werden. 
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Bei der Betrachtung der Genauigkeit wollen wir zwei Extremfalle 
unterscheiden: Im ersten Fall moge n, 0, >>) Qo sein. Hierzu braucht 
der Ionisationsgrad gar nicht besonders hoch zu sein, weil der Ionen- 
querschnitt Q, selbst schon durch das weitreichende CouLomB-Feld der 
Ionen im allgemeinen viel groBer als der RAMSAUER-Querschnitt Q, ist. 
Ist also obige Ungleichung erfiillt, so kann man naherungsweise ”, gegen 
n, kiirzen und es bleibt nur eine geringe Temperaturabhangigkeit tibrig. 
Ist im zweiten Fall »,Q,>",@Q., dann bleibt das nach der SAHA- 
Gleichung (11) sehr temperaturempfindliche », im Integral stehen, und 
es ist unerlaBlich, m, direkt zu messen, woraus anschlieBend T mit Hilfe 
Gl. (41) berechnet werden kann. Wenn man so verfahrt, ist auch in 
diesem Fall die Methode recht genau. 


III. Der Ionenquerschnitt. 
a) Die MeBmethode am Hochleistungsbogen. 

Zur Messung der Ionenquerschnitte empfiehlt es sich, einen Bogen 
zu verwenden, bei dem der Gesamtquerschnitt )Q, der Neutralteilchen 
zu vernachlassigen ist oder nur eine kleine Korrektion darstellt, also 
einen Bogen moglichst hohen Ionisationsgrades. Diese Bedingung er- 
fiillt in hohem MaBe der wasserstabilisierte Hochleistungsbogen (GER- 
DIEN-Bogen [1}), der in einem durch rotierendes Wasser gebildeten 
Kanal von 2,3mm @ brennt und bei Stromstarken bis zu 1500 Amp 
Achsentemperaturen bis tiber 50000° K aufweist [2], [3). In diesem 
Bogen sind sowohl die Feldstarken [4) als auch die radialen Temperatur- 
funktionen fiir eine Reihe von Stromstarken gemessen worden. Damit 
kann nach Gl. (1) zunachst der Leitwert/em = J/E angegeben werden. 
Ferner k6énnen bei jeder Temperatur die verschiedenen Teilchen- 
zahlen/cm? fiir H- und O-Atome, fiir ein- und mehrfach geladene O-Ionen 
und fiir die Elektronen unter Verwendung der SAHA-Gleichung (11), der 
chemischen Formel fiir Wasser H,O und des Gesamtdruckes berechnet 
werden. SchlieBlich benétigen wir noch den yon GvOSDOVER [5] u. a. 
theoretisch ermittelten Ionenauerschnitt 


ee to iS ez 
7 2 (kT)? tn e my (13) 


der sich nach der RuTHERFORDschen Streuformel auf Z-fach geladene 
Ionen mit Hilfe des Zusammenhangs 


Qz, =2°Q, (14) 

libertragen laBt. 
Mit diesen Beziehungen kénnen wir das Integral in (12) ausfiihren 
und die rechte Seite der Gleichung mit dem fiir jede Stromstirke ge- 
messenen Leitwert/cm auf der linken Seite von (12) vergleichen. Diese 


Ionen- und Atomquerschnitte im Plasma. 123 


Auswertung ist in einer friiheren Arbeit von uns [6] ausgefiihrt worden 

und fihrt nur dann zur Ubereinstimmung beider Seiten, wenn man fiir 
: “tt Cth ei 

den Gvospov ER-Querschnitt Q“ einen um den Faktor 0,55 verkleinerten 

Wert einsetzt: 


Q. = 0,55 Qt. (15) 


Obwohl dieses Ergebnis nach wie vor als richtig anzusehen ist, waren 
doch gewisse Unstimmigkeiten vorhanden, zu deren Aufklirung eine 
Priifung der gemessenen Temperaturver- 
teilungen notwendig war. 


Die seinerzeit verwendete und auch hier 
wieder herangezogene TemperaturmeBme- 
thode stammt von R.W. LarENnz [2] und 
geht auf die bekannte Abhiangigkeit der 
Intensitat einer Spektrallinie von der Tem- 
peratur zuriick: 


tp ~~ n,(T)e ear. (16) 


Temperatur T 


Danach steigt die Intensitat der Spektral- 
linien zunachst nach MaBgabe des Bortz- 
MANN-Faktors an, bei hohen Temperaturen pig 4, Temperaturverteilungen im 
jedoch tiberwiegt die Abnahme derTeilchen- _wasserstabilisierten __Hochleistungs- 
zahl m,(T) durch die Jonisierung dieser = oa er ee 
Teilchen, so da8 dann die Intensitat nach urven sind ans dem relativen Inten- 
Durchlaufen eines Maximums trotz steigen- ee eat ae pats Un. 
der Temperatur wieder abnimmt. Man kann — the oes Bogens im Kanal. Die Be- 
nun die Temperatur, bei der das Maximum _ korrigierte Verteilung bei 200 Amp. 
liegt, aus Gl. (16) berechnen, und dem- 

jenigen Bogenradius, in dem die maximale Intensitat der Spektrallinie fest- 
gestellt wird, diese Temperatur zuordnen. Weiter kann man von diesem 
Fixpunkt aus nach innen und nach auBen durch Vergleich der relativen 
Intensitaten den radialen Temperaturverlauf in der Umgebung des 
Maximums ermitteln. Bei der Messung am Hochleistungsbogen zeigte 
sich, daB nicht nur die Linien der Neutralteilchen H und O, sondern auch 
die Linien der einfachen Sauerstoffionen ein Maximum und eine Ein- 
senkung in der Bogenachse aufweisen. Im Extremfall waren sogar noch 
die Linien des zweifach ionisierten Sauerstoffs eingesenkt. All diese 
Maxima dienten als Temperaturfixpunkte, mit deren Hilfe die radialen 
Temperaturverteilungen gezeichnet werden konnten. Bei der Bestim- 
mung des Temperaturverlaufs in der Umgebung der Maxima aus ein 
und derselben Spektrallinie traten jedoch Abweichungen auf dergestalt, 
daB der Temperaturverlauf immer flacher verlief als die durch die Fix- 


punkte gelegte Kurve hatte erwarten lassen (Fig. 1 punktiert). Diese 


0 05 jomm 
Bogenradius p 
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Unstimmigkeit wurde seinerzeit zuriickgefiihrt auf Verwaschungen der 
Spektrallinien in der Langsrichtung durch die mdéglichen Unebenheiten 
des Wassers an seiner inneren Kanaloberflache, denn die Intensitats- 
messung geschah von der Seite durch ein Quarzfenster und durch die 
den Kanal bildende Wasserschicht. Die Verwaschung andert jedoch 
nicht den durch mehrere Fixpunkte festgelegten Temperaturverlauf, 
weil die Lage der Maxima durch sie praktisch nicht beeintrachtigt wird. 
Nur bei den mittleren Stromstarken (@ 200 Amp), bei denen auBer der 
horizontalen Tangente in der Bogenachse und dem Randwert nur noch 
ein Fixpunkt zur Festlegung des radialen Temperaturverlaufs existiert, 
sind die Verteilungen unsicher. Dieses war ein Grund, der die Uber- 
priifung der Messung bei 200 Amp veranlaBte. 


b) Uberpriifung friiherer Messungen. 

Um eine mégliche Streuung des Lichtes an der inneren Kanalober- 
flache zu umgehen und die bei seitlicher Beobachtung notwendige Um- 
rechnung von der beobachteten Querverteilung der Intensitat auf die 

PreBlofzvtir 


4 ustrerendes Wasser 


4 t-O1R  Srevlichtachats 
t o SF 
Endplatte A Plexiglas e pektrographenspa. 


‘ tangentiale Wasserzufir 
Fig. 2, Anordnung zur end on-Beobachtung des 200 Amp-Hochleistungsbogens. 


radiale Verteilung zu vermeiden, wurde in Verlangerung der Bogenachse 
(end on) beobachtet. Dabei ist allerdings der Umstand von Nachteil, 
daB aus dem Bogenkanal nach beiden Seiten ein intensiver Plasmastrahl 
austritt, in dem die Temperatur stetig abfallt, so daB bei end on-Be- 
obachtung nicht nur das Leuchten im eigentlichen Bogenkanal, sondern 
auch in den inhomogenen Schichten zwischen den Kanalenden und den 
Elektroden wirksam wird. Um diesen Ubelstand zu beseitigen, sind die 
durchbohrten Kohleelektroden ganz dicht an den Kanal herangefiihrt 
(Fig. 2). Mit Hilfe eines PreSiuftstrahles durch die Elektrodenbohrung 
in Richtung auf den Bogen wird dafiir gesorgt, daB das austretende 
Plasma unmittelbar nach seinem Austritt aus dem Kanal zur Seite ge- 
drangt wird. Dadurch wird die inhomogene Schicht so stark verkiirzt, 
daB sie klein bleibt gegeniiber der Gesamtlinge des Bogens von 2 cm. 

Trotz dieser VorsichtsmaBnahme sind doch noch Médglichkeiten 
gegeben, die eine Verwaschung des Intensititsverlaufes entlang einer 
Spektrallinie hervorrufen kénnen. Diese sind: 


1. Mangelhafte Scharfentiefe trotz starker Abblendung der Optik, 


2. mangelhaftes Auflésungsvermégen wegen der starken Abblendung 
der Optik, 
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3. Ablenkung der aus dem Bogen austretenden Lichtstrahlen durch 
die Schlierenbildung in den inhomogenen Schichten vor den Elektroden. 


Fur die Abblendung der Optik wurde ein solcher Offnungswinkel ge- 
wahlt, daB die Tiefenscharfe und das Auflésungsvermégen einander gleich 
sind, weil dann die Abbildungsschirfe optimal ist. Fiir die gewahlte 
Anordnung betrug der Offnungswinkel 1:70, entsprechend einem Auf- 
l6sungsvermégen von 0,1mm am Ort des Bogens. Dieses Auflésungs- 
vermogen wurde nicht nur durch die Rechnungen gefunden, sondern 
auch durch Messung bestatigt, indem ein diinner Platindraht an der 
Stelle der Bogenachse zum Gliithen 
gebracht und sein Bild in der Ebene 
des Spektrographenspaltes aufge- 
nommen und vermessen wurde. 


Um die Wirkung der Ablen- 
kung durch die Schlieren zu priifen, 
wurde das Licht einer Quecksilber- 
lampe, aus der durch Zwischenab- 
bildung ein sehr kleiner Teil mit 
Hilfe einer feinen Lochblende her- 
ausgeblendet war, durch den bren- 
nenden Bogen hindurchgeschickt 
und die GréBe der punktférmigen 
Quecksilberlinien im Spektrum 
des Hoch leistungsbogens gemessen. Fig. 3. Zur Entzerrung gemessener Intensitiatsver- 

i . = i teilungen. Der Spalt ist in der a-Achse zu denken, 
Dabei ergab sich, daB diese 
Punkte zwar etwas gréBer waren 
als bei den Aufnahmen ohne brennenden Bogen, daB aber das Auflésungs- 
vermégen der gesamten Apparatur dadurch nicht wesentlich beeintrach- 
tigt wird, zumal wenn man bedenkt, daB das Quecksilberlicht beide 
inhomogene Schichten durchdringen mub, wihrend das Licht aus dem 
Bogen nur eine dieser Schichten durchstrahlt. (Von der Wiedergabe einer 
solchen Aufnahme muBte abgesehen werden, weil die auBerordentlich 
feinen Hg-Linienpunkte in der Reproduktion untergehen wiirden.) 


Trotz der geringen Beeintrachtigung des Auflésungsvermégens durch 
die drei oben genannten Vorginge wurde die gemessene IJntensitits- 
verteilung mit der gemessenen Apparatefunktion auf folgende Weise 
entzerrt: Mége in Fig. 3 der groBe Kreis ein Symbol fiir das fehlerfreie 
Bild des Bogens in der Spaltebene bei end on-Beobachtung sein und moge 
i(r) df die entsprechende wahre Intensitit im Punkte O sein, Dann wird 
durch die Streuvorgiinge ein Teil dieser Intensitat zum Punkt P des 
Spaltes hingetragen. Die GréBe dieses Bruchteils hangt von der Ent- 
fernung O— P=r' ab. Dieser Bruchteil als Funktion vor 7’ ist die 
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Apparatefunktion 4 (7’). Von O nach P gelangt also die Intensitat 
>) 4 (r’) df. Summiert man die Beitrage aller Flachenelemente 


aj=Ffarae (17) 


zur Intensitit in P, dann resultiert die gemessene Intensitat auf dem 


Spalt oe 
| 


ai. r'\i(r)r' dr’ dg’. (18) 


Da A(r’) durch Messung und Rechnung bekannt ist, 1aBt sich prinzipiell 
aus der gemessenen Intensitatsverteilung 7, (a) die wahre Intensitatsver- 
teilung i(r) ausrechnen. Die praktische Lésung dieser Integralgleichung 

seschieht durch Probieren, und zwar so, 


y44 


daB man fiir 7(v) verschiedene, etwas steiler 
als 7.(a) verlaufende Funktionen versuchs- 


Hy 48601 A 


Hy 


OIL46 


weise in das Integral von (18) einsetzt. 
Kommt dann aus der graphischen Integra- 
tion die entlang des Spaltes a und damit 
entlang einer Spektrallinie gemessene Ver- 


teilung heraus, dann ist die eingesetzte 


Fig. 4. End on-Aufnahme 1 
Anordnung Fig 


: Verteilung 7(r) die ungestérte, also auf 
Streuvorgange korrigierte Verteilung. 
Die Durchfiithrung einer Messung soll an dem folgenden Beispiel 
erlautert werden. Fig. 4 zeigt die end on-Aufnahme vom 200 Amp- 
Bogen mit 2,3mm Kanaldurchmesser und 2cm Lange. Man sieht 
die verbreiterten und vor allem eingesenkten Linien H, und H,, und 
die in der Bogenmitte auftauchenden O II]-Linien. Die Photometrierung 
von H,; entlang einer Kuppe und die Umrechnung mit der Schwar- 
zungskurve ergibt die relative Intensitatsverteilung in Fig. 5. Das 
gleiche gilt fiir die Sauerstofflinie O II 4651 A. An Hand der theoreti- 
schen Intensitatskurven als Funktion der Temperatur in Fig. 6 gemaB 
Gl. (16) laBt sich zunachst die Temperatur bei dem Bogenradius fest- 
legen, bei dem das Maximum von H, beobachtet wird. Von hier aus legt 
man dureh Vergleich der relativen Intensitaten den vorlaufigen Tem- 
peraturverlauf nach auben und innen fest. Bei der O II-Linie verfahrt 
man ahnlich, indem man ihre Intensitat auf die beim Maximum von H, 
vorhandene bezieht. Die so gefundenen radialen Temperaturkurven a 
in Fig. 7 divergieren sehr stark und passen auch keineswegs mit der nach 
ig. 1 zuerwartenden tiberein. Den Grund hierfiir wird man zuniachst darin 
suchen, daB die Intensitatsverteilung in den Linien durch die Apparate- 
funktion verwaschen ist. Man fiihrt also das Entzerrungsverfahren nach 
Gl. (18) fiir jede Linie durch und wiederholt die Temperaturbestimmung 
(IXKurven 6 in Fig. 7). Wenn auch dadurch der Temperaturverlauf ein 
wenig steiler geworden ist, so klafft immer noch eine groBe Liicke 
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zwischen der aus H, und der aus der O I-Linie gewonnenen 7-Kurve, 
und beide legen immer noch viel flacher, als es die aus den Intensitiits- 
maximis gewonnene 7-Verteilung in — 76 
Fig. 1 verlangt. Da also trotz der vor- 
genommenen Korrektionen auf Streu- ™ 
effekte diese Diskrepanz bestehen 
bleibt, kénnen Verwaschungen durch 
Streueffekte nicht allein dafiir verant- 
wortlich gemacht werden. Dieses Er- 
gebnis der experimentellen Uber- 
priffung zeigt, daB alle fritheren und 
jetzigen Intensitatsmessungen in Ord- 
nung sind, zumal die Lage des H,- 
Maximums reproduziert werden konnte. 
Ohne uns nun zuniachst davon 
Rechenschaft zu geben, wodurch sonst 
die Verwaschungen hervorgerufen sind, 
k6nnen wir ihre Wirkung dadurch i, < Gemessene Intensitatsverteilungen fir 
eliminieren, daB wir bei der Entzerrung = 4%. und 011 4651 A mit der Anordnung Fig. 2. 
neben der durch Streuung bedingten 
Apparatefunktion A (r’) eine weitere zusatzlich einfiihren mit der 
Forderung, da die damit entzerrten Intensitatsverteilungen von Hg, 
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Fig. 6. Nach Gl. (16) berechnete, relative Intensitat Fig. 7. Unkorrigierte (a) und nur mit Apparate- 
als Funktion der Temperatur fiir Hg und O II 4651 A. funktion korrigierte (b) Temperaturverteilungen 


aus den Kurven Fig. 5 und 6. 


und OI 4651 A ein und dieselbe Temperaturverteilung ergeben. Als 
zusatzliche Apparatefunktion wahlen wir eine Rechteckfunktion, deren 
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Breite in der vorgeschriebenen Weise zu bestimmen ist. indimensit: 
nales Beispiel: eine Spektrallinie werde mit einem relativ breiten Spalt 
photometriert ; man bestimme aus der wahren Linienform und der ver- 
waschenen Registrierkurve die Breite des Photometerspaltes.) Durch 
das oben beschriebene Probieren kommt man schlieBlich zu einer korri- 
gierten Intensitatsverteilung, die in Fig. 8a und b stark ausgezogen ist 

Die Richtigkeit dieser Kurve veri- 
fiziert man dadurch, daB man sie 
wieder zunachst mit der bekannten 
Apparatefunktion und dann nach- 
einander mit Rechteckfunktionen 


relative Intensitol ¢ (ry 


relative Intensital t(r) 


0 02 OF 06 Q&mm ~~ vReeseyve 
Bogenradius 7 Sogesrodus 7 mm 
a b 


Fig. 8a u.b. Die richtige, entzerrte Intensitatsverteilung geht durch Verzerrung mit Cer Apparatefunktion 
A(r) und einer Rechteckfunktion der Breite 0,38 mm in die gemessene Verteilung ther. a Hg. b 4651 A. 


verschiedener Breite verzerrt, bis die gemessene Verteilung heraus- 
kommt, was fiir beide Linien bei einer Rechtecksbreite von 0,38 mm 
der Fall ist. Die hier angenommenen Ausgangskurven sind also auf 
wahre Intensitatsverteilungen richtig korrigiert, wenn sie noch beide 
denselben Temperaturverlauf ergeben, der dann als der richtige an- 
zusehen ist. Auch diese Forderung wird von den beiden Ausgangs- 
funktionen erfiillt. Der so gefundene wahre T-Verlauf ist in Fig. 1 als 
ausgezogene Kurve bei 200 Amp mit eingezeichnet. Obwohl uns diese 
Kurve etwas zu steil ausgefallen zu sein scheint, fiigt sie sich doch 
einigermaBen in die friiher gemessenen Kurven ein. 


Die Ursache fiir die Verwaschung sehen wir in einer gewissen Un- 
tuhe des Bogens im Kanal, die sich entweder in einem Pendeln um etwa _ 
+ 0,2 mm, in einer schlanken Schraubenbahn um die Achse oder in einer 
gewissen Turbulenz von gleicher GréBenordnung auBert. 


| 
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Die Auswertung der korrigierten Temperaturverteilung beziiglich des 
lonenquerschnittes gema8 Gl. (12) fiihrt zu einem Q,/Q"-Verhiltnis von 
0,56 gegeniiber friiher von 0,53 und schlieBt sich damit noch besser als 
vorher dem Mittelwert von 0,55 an. Die geringe Anderung des Zahlen- 
wertes fiir den lonenquerschnitt trotz des scheinbar viel gré®eren Unter- 
schiedes der Temperaturverteilung mit und ohne Korrektion in Fig. 14 
beweist die Unempfindlichkeit der ganzen Methode auf irgendwelche 
Fehler. 


c) Neue Bestimmung des Ionenquerschnittes 
in Verbindung mit der Sp1tzERschen Theorie. 


In der neueren Zeit sind nun weitere theoretische Arbeiten tiber die 
elektrische Leitfahigkeit von ionisierten Gasen erschienen, von denen 
besonders die von L. Spitzer und seinen Mitarbeitern [7] zu nennen ist. 
Hier ist auch die Wechselwirkung der Elektronen untereinander mit- 
beriicksichtigt, die zwar insgesamt keine Impulsdnderung der Elektronen 
ergibt, aber dennoch wegen der Anderung der Geschwindigkeitsvertei- 
lung die elektrische Leitfahigkeit beeinfluBt. Ferner wird in der Arbeit 
von SPITZER, genau wie in unserer friiheren Arbeit [6] die effektive 
Ionenladungszahl Z 

-e D122 is L172? 

Ne Lz Z 

fiir ein Gasgemisch aus mehrfach geladenen Ionen und Elektronen ein- 

gefiihrt. Diese Erweiterung ist tiir unseren Hochleistungsbogen deswegen 

besonders wichtig, weil hier bis zu dreifach geladene Ionen vorkommen. 

Die elektrische Leitfahigkeit eines vollstandig ionisierten Plasmas be- 
schreibt Spitzer durch die Gleichung 


(19) 


a 
pela dds Ye, (20) 
al/ > m, €? Z \n (q’ C?) 
worin schon 
C2 BAT (21) 


Mz 


gesetzt ist. y, stellt den Korrektionsfaktor gegentiber dem LoRENTZ-Gas 
dar. (Beir Lorentz-Gas wird die Elektronenwechselwirkung unter- 
einander nicht beriicksichtigt und die Ionen werden als ,,in Ruhe befind- 
lich’ betrachtet.) 

Fir Z =1 ist yg = 0,5816, 

fir Z=2 ist yg = 0,6833 und 

fir Z = oo wird yg =1 (LORENTZ-Gas). 

In der Spitzerschen Arbeit erfolgt die Ableitung der Leitfahigkeit 
auf so allgemeiner Grundlage, da sowohl der Zahlenfaktor als auch die 
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Temperaturabhangigkeit als zuverlassig angesehen werden kénnen. Da- 
gegen besteht eine gewisse Willkiir in der Festlegung des Gliedes 
In (q'C?). SPITZER verwendet fiir ein Plasma hoher Elektronendichte, 
das gerade bei dem von uns durchgemessenen Hochleistungsbogen vor- 
liegt, den von CHAPMAN und CowLinG [8] tbernommenen Ausdruck 


In (* es dj, worin d bei alleiniger Betrachtung von Zweierst6Ben 
é 

zwischen Elektron und Ion die effektive Reichweite des COULOMB- 

Feldes bedeutet. CHAPMAN und CowLING setzen dafiir den mittleren 


Abstand zweier benachbarter Ionen, namlich 
d=n-5, (22) 


Dieser Wert ist sicher zu hoch gegriffen, weil spatestens in der Mitte 
zwischen zwei benachbarien Ionen das Feld des Nachbarions wirksam 
wird. Aber auch schon vor der Mitte wird sicher ein EinfluB des Nachbar- 
feldes vorhanden sein, so da8 man fiir die effektive Reichweite des 
CouLoMB-Feldes eher 

d<in-3 (23) 


erwarten sollte. Jedenfalls lhegt die gréBte Unsicherheit der Theorie in 
dem Zahlenfaktor des In-Gliedes. Deswegen ersetzen wir ihn durch einen 
> ; : 

- =] und bestimmen ihn 
Ze ns 
so, daB die nach Gl. (2) integrierte Leitfahigkeit (20) den gemessenen 
Leitwert pro cm ergibt. Da wir es im Hochleistungsbogen mit einem 
vollstandig ionisierten Plasma zu tun haben, kénnen wir direkt den von 
SPITZER angegebenen Ausdruck fiir die Leitfahigkeit mit dem Unsicher- 
heitsfaktor « verwenden, ohne schon jetzt auf die Wirkungsquerschnitte 
zuriickgehen zu mitissen. Um bei der Bestimmung des Faktors « vdllig 
sicher zu gehen, sind fiir die Integrationen nur diejenigen Temperatur- 
verteilungen herangezogen worden, die sich auf mindestens zwei Tem- 
peraturfixpunkte aus entsprechenden Linienmaxima stiitzen, und auBer- 
dem die oben erwahnte sehr genau nachgepriifte Verteilung bei 200 Amp. 
Das Ergebnis der Integrationen ist in der folgenden Tabelle 1 zusammen- 
gefaBt und ergibt den mit Riicksicht auf die MeBgenauigkeit abgerun- 
deten Wert «1. 

Die Abweichungen der Einzelmessungen vom Mittelwert fiir « er- 
scheinen teilweise recht hoch, sind aber unerheblich, wenn man bedenkt, 
da « unter dem In steht. Der tatsichliche EinfluB auf die Leitfahigkeit 
kommt in der letzten Spalte von Tabelle 1, in der das Verhaltnis der 
gemessenen Leitfahigkeit zu der nach Gl. (20) mit «=4 berechneten 
steht, und in der Fig. 9 zum Ausdruck, in der die Leitfihigkeit nach 
SPITZER mit «=1 als ausgezogene Kurve in Abhingigkeit von der 
Temperatur aufgetragen ist, wahrend die einzelnen Kreise die gemessenen 


Unsicherheitsfaktor «, schreiben also In [x 
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Tabelle 1. 
ee 

7 RE | Pe > o= 

Amp | V/cm ee pled 
202 106 0,86 26000 1,04 
250 114 0,78 26000 1,07 
300 121 0,89 26 500 1,03 
350 128 0,90 27 500 | 1,03 
400 136 1,00 28 500 1,00 
600 166 1,02 32 500 0,98 
800 196 1,09 35500 0,97 

1000 226 1,38 40000 0,92 

1450 294 2,15 44.000 0,84 


Mittel 1,1 


(Die uber die Leitfahigkeit gemittelte Temperatur T ergibt sich aus der 


Aol de: : ~ [r@)a(T)rar \ 
Definitionsgleichung T = gn 
fordr 
0 


o-Werte wiedergeben. Ferner ist gestrichelt die mit «=4 (nach 
CHAPMAN-COWLING) berechnete Leitfahigkeit eingetragen, die etwa 30% 


unterhalb der Kurve mit « =1 verlauft % 300 Siemens/em. 
und somit auBerhalb der MeBgenauig- 

keit liegt. SchlieBlich sind noch zwei &” 

MeBpunkte als Kreuze eingetragen, = e 

die sich aus der Auswertung der von 3 ¥% 

LARENZ gemessenen Temperaturver- 2 3 

teilungen ergeben haben. Leider haftet S$ , 

0 0 20,5 30 ~~ -y0-107"K 50 


aber diesen Punkten eine ziemliche epee E 
Unsicherheit an, weil die Feldstarke an Fig. 9. Leitfahigkeit eines reinen Plasmas 
dem gemessenen Ort (2 mm vor dem (H,0) nach L. SPITZER (7] mit «=1 (aus- 
Kanal) nicht bekannt ist. Wir haben ret rr Gena 
aus der angegebenen Temperaturver- 

teilung einen effektiven Kanaldurchmesser von 5 mm extrapoliert und 
nach [4] eine Feldstarke 50 V/cm eingesetzt. Ob der systematische 
Gang der Abweichung mit der Stromstarke reell ist, kann im Rahmen 
dieser Untersuchung nicht festgestellt werden, da die MeBgenauigkeit 
von etwa +15% (Fehler in der Temperaturbestimmung sowie Strom- 
starke- und Feldstarkemessung) nicht dazu ausreicht. 

Der Mittelwert «=1 bedeutet, daB das CouLtoms-Feld bei 1/, des 
Gesamtabstandes zweier benachbarter Ionen abzuschneiden ist, wenn 
man mit den Messungen in Ubereinstimmung kommen will. 

In diesem Zusammenhang mag noch interessieren, daB G. BuRK- 
HARDT [20] aus Dispersionsbetrachtungen fiir das reine Plasma einen 
a-Wert von «=1/1,44=0,7 fordert, was einem Abschneiden bei etwa 


1/, Ionenabstand entsprechen wiirde. 
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Bei kleineren Ionendichten ist noch der Einflu8 der entfernteren 
Ionen mit zu beriicksichtigen. Dann hat an Stelle des hier benutzten 
In-Gliedes ein anderes zu treten, das aus den DEByeEschen Vorstellungen 
entstanden ist, namlich In (3/2! n} e*) (kT/2)?. Da sich aber aus unseren 
Messungen ein Abschneideradius von 1/4-Ionenabstand ergeben hat, 
sind offenbar die nahen, wenn auch seltenen StéBe wirksamer als die 
hiufiger stattfindenden fernen Vorbeigénge der Elektronen, wodurch 
nachtraglich die von uns getroffene Wahl des In-Gliedes seine Berech- 
tigung findet. 

Wenn nun ein Plasma vorliegt, in dem auBer den Ionen auch noch 
Neutralteilchen an der Bremsung der Elektronen beteiligt sind, dann 
muB man sich auch fiir die Ionen wieder der Hilfsvorstellung des Wir- 
kungsquerschnittes bedienen, um aus beiden Querschnitten die Gesamt- 
leitfahigkeit des Plasmas nach Gl. (10) berechnen zu kénnen. Um also 
aus Gl. (20) einen effektiven Ionenquerschnitt zu ermitteln, identifi- 
zieren wir sie mit der Leitfahigkeitsgleichung (10) fiir ein vollstandig 
ionisiertes Gas (m)Q)*0) und erhalten fiir Z-fach geladene Ionen 


Cus 


0,567 Ze kT 

ve OTF | Zag)” uy 
Fir Z=1 mit ye =0,5816 folgt daraus innerhalb der MeBgenauigkeit 
von +15% der einfache Ausdruck 


‘2 s mys n 
Q. ae (A T)? ane a |. (25) 


der etwa halb so gro8 ist, wie der von GVOSDOVER berechnete Wert in 
Ubereinstimmung mit Gl. (15). Uns scheint Gl. (25) die geeignetste For- 
mel fiir Gasentladungsrechnungen zu sein. 


IV. RAMSAUER-Querschnitte fiir Wasserstoffatome. 
a) Messungen am normalen wasserstabilisierten Bogen. 


Nachdem nunmehr die Wirkungsquerschnitte der Ionen gegeniiber 
Elektronensto8 bekannt sind, kann man aus der Leitfaihigkeit eines 
Wasserdampfplasmas, in dem sich Ionen und Atome in vergleichbarem 
MaBe an der Bremsung beteiligen, die RAMSAUER-Querschnitte der Atome 
ermitteln. Hierzu eignet sich der von H. MAECKER [9] vorgeschlagene 
normale wasserstabilisierte Bogen, der in einem verhaltnismaBig weiten 
Rohr mit nicht sehr hohen Stromstarken brennt. Einen solchen Bogen 
mit 8mm Rohrdurchmesser, 5 cm Lange, 50 Amp und 35 V/cm hat 
G. JURGENS [10] sorgfialtig durchgemessen und in der Achse eine Tem- 
peratur von 12650° Kund eine Elektronenkonzentration von 8,4+101%em~3 
gefunden. 
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An diesem Bogen haben wir nun die radiale Temperaturverteilung 
aus dem Intensitatsabfall der Linien H, und H , und des Kontinuums 
zwischen diesen beiden Linien in Verbindung mit dem theoretischen 
Intensitatsabfall Gl. (16) gemessen. Dabei wurde sowohl die side on- 
Beobachtung mit anschlieBender Umrechnung der gemessenen Quer- 
verteilung auf radiale Intensitatsverteilung mit Hilfe der ABELschen 
Integralgleichung als auch die end on-Beobachtung herangezogen. In 
beiden Fallen wurde besondere Sorgfalt auf die Ausschaltung von 
Fehlern durch Verwaschungseffekte, also unebene Wasser- oder Quarz- 
oberflache, Tiefenscharfe und Beugung verwendet. Wenn diese Ein- 
fliisse auch durchwegs gering waren, so wurden 0K 
sie doch korrigiert. Die Verwendung eines a 
Spektralapparates eriibrigte sich, weil die frag- 
lichen Linien und erst recht das Stiick Kon- 
tinuum durch geeignete Filter ausgeblendet 
werden konnten. Das Ergebnis einer groBen 
Zahl solcher Messungen ist in Fig.10 darge- 
stellt. Mit dieser Temperaturverteilung kénnen Seat 
wir nun alle temperaturabhangigen GréBen in 0 7 am 
der Leitfahigkeitsgleichung (12), namlich 1,, 1 4 Sco a 
und m, Q, berechnen und Q, so bestimmen, daB RCE ee 
aus der Integration der gemessene Leitwert  wsserstabilisierten Bogen mit 
pro cm //E herauskommt. Fiir Q, ergibt sich =outtt aie anette ar 
der iiberraschend h ohe Wert von tensitatsabfall und der wesent- 


lich breitere Leitfahigkeits- 
abfall, beides relativ. 


72 


Oia 02740 4c, (26) 


Dieser ist so viel gr6Ber als der nach der ersten Bornschen Naherung [11] 
berechnete Wert von 42 a?=3,5- 101% cm?, daB diese Diskrepanz der 
Hauptgrund fiir die sorgfaltige Uberpriifung aller Messungen war. In- 
zwischen haben allerdings quantenmechanische Rechnungen von L. Mo- 
WER [12] und von H. S. W. Massey und B. L. Motserwitscu [13] den 
Wirkungsquerschnitt von H-Atomen fiir den in Frage kommenden 
Temperaturbereich auf 25 - 10-16 cm? heraufgesetzt, wodurch der Unter- 
schied um nahezu eine Zehnerpotenz verringert worden ist. 


b) Messungen am wasserstabilisierten Hochdruckbogen. 


Noch genauer als im vorgenannten Bogen kénnen die Wirkungsquer- 
schnitte der Neutralteilchen im Wasserdampfplasma des wasserstabili- 
sierten Hochdruckbogens ermittelt werden, der von Tu. Peters [J] 
eingehend untersucht worden ist. Bei diesem Bogen entsteht der Stabi- 
lisierungskanal durch Rotation eines teilweise mit Wasser gefiillten 
DruckgefaBes. Der Lichtbogen erzeugt dann durch Wasserverdampfung 
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selbsttitig den Druck, der bei einem maximalen Leistungsumsatz von 
100 kW/cm (200 Amp, 500 V/cm) auf 1000 at ansteigen kann. 

Die aus der Linienform von H, direkt 
gemessene Elektronendichte fiihrte in dem 
durchgemessenen Druckbereich von 50 bis 
300 at zu einem Ionisationsgrad von 1 bis 2%. 
Die aus der Elektronendichte tiber die SAHA- 
Gleichung berechneten Achsentemperaturen 
iH +H per ergeben einen geringfiigigen Abfall der Tem- 
(a-s060A) peratur mit steigendem Druck bei konstanter 
[ (p-120at) Stromstarke und liegen bei etwa 12000° K. 
Unter diesen Umstanden kann der Wirkungs- 
querschnitt der Ionen in Gl. (10) vollstandig 
vernachlassigt werden. Die Leitfahigkeit wird 
allein durch die neutralen Atome bestimmt. 

ahead Deswegen ist der an diesem Bogen gemessene 
Fie asmllite,ecdedhae fom RAMSAUER-Querschnitt der Neutralteilchen 
H--Ion erzeugte Kontinuumsin- frei von Fehlern in den lJonenquerschnitten 
titicierten ‘Hockdeacktesea, Und daher zuverlassig. Der relative radiale 
als Funktion der Temperatur Temperaturverlauf wird durch Messung der 
iia radialen Strahlungsdichteverteilung im Kon- 
tinuum zwischen H, und H,; nach Anschlu8 an die Achsentemperatur 
ermittelt. Da das Kontinuum zum gréBten Teil vom H~-Ion aus- 
gestrahlt wird, mu8 die Berechnung der Kontinuumsintensitat als 
-797°K Funktion der Temperatur in erster Linie auf 
diesem ProzeB aufgebaut werden (Fig. 11). Mit 
dieser Kurve erhalt man aus der gemessenen 
radialen Intensitatsverteilung den Temperatur- 
verlauf, wie er fiir einen Fall in Fig.12 dar- 
gestellt ist. Die Integration gemaB Gl. (12) 
ergibt dann direkt die Wirkungsquerschnitte 
der neutralen Atome. Die Ergebnisse dieser 
Auswertungen sind in Tabelle 2 wiederge- 
geben. 


Emissionskoefizent des Kontinuums 


Ss 8 


Temperatur T 
%& RD ® 


7 
ogenradius 7 
Fig. 12. Relativer radialer Ver- 
lauf des Kontinuums und absolute Der Mittelwert von Oo = 94; 10°48 cm2 


Temperatur im Hochdruckbogen, : “a : 
stimmt tberraschend gut mit dem Wert Q,= 


92+410°' cm® des vorigen Abschnittes iiberein, so daB als Mittel fiir 
beide Bogenarten resultiert: 


5H O = 93 + 10716 cm® + 15%, (27) 
Wegen dieser guten Ubereinstimmung ist es berechtigt, auch aus den 


Messungen am normalen wasserstabilisierten Bogen einen «-Wert fiir die 
lonen zu entnehmen; er ist in Fig. 9 mit eingetragen. 
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Tabelle 2. 
a De ee 
E n } * 
Atm Amp | Vem 407%7 om-* | °K hota 
| | 

50 100 300 6,0 | 42050 106 

80 100 380 aya 11850 99 
100 100 420 8,8 12000 110 
120 100 460 72% 11 500 80 
150 100 500 8.8 11600 88 
300 150 450 20,0 12 500 | 79 


Mittel ~ 94 


Nun ist noch zu bedenken, daB im Wasserdampfplasma ein Gemisch 
von Wasserstoff- und Sauerstoffatomen im Verhiltnis 2:1 vorliegt. Aus 
Messungen am Hochstrombogen in Luft, iiber die im nachsten Kapitel 
berichtet wird, folgt, da8 der Wirkungsquerschnitt von neutralen Sauer- 
stoffatomen gegeniiber ElektronenstoB mit etwa Q° = 20+ 10716 cm? an- 
zusetzen ist. Teilen wir also den Gesamtquerschnitt auf, entsprechend 


ty Oy = n# QE + 190, (28) 
so ergibt sich mit 
Mt, = nt + 9 
und 
ni — 2n° 


ein Wasserstoffquerschnitt von 
OM = $0, — $ O° = 130-10 cm? + 25%. (29) 


Wenn auch dieses Ergebnis noch wesentlich héher liegt als die bisher 
berechneten Werte, so sind die iibereinstimmenden Aussagen der obigen 
Untersuchungen doch zu eindeutig, als daB an ihrer Richtigkeit gezweifelt 
werden kénnte. 


V. RAMSAUER-Querschnitte fiir C-, N- und O-Atome 
aus Messungen an der Sdule des Hochstromkohlebogens. 


Kiirzlich hat H. MAECKER [15] die Elektronendichte in der Sdule des 
Hochstromkohlebogens von 200 Amp aus der STARK-Effektverbreiterung 
von H, und der Absolutintensitat des Kontinuums als Funktion des 
Bogenradius gemessen und die absolute radiale Temperaturverteilung 
nach verschiedenen Methoden bestimmt (Achsenwerte: 1, = 7 - 10'*cm~°, 
T =11000° K). AuBerdem ist fiir diesen Bogen von W. FINKELNBURG 
und S. M. SEGAL [16] die Saulenfeldstarke zu 10 V/cm angegeben worden. 
Diese MeBergebnisse gestatten es, die Leitfahigkeitsgleichung (12) durch 
iterierte Integration nach Q, aufzulésen, so daB man einen mittleren 
Wirkungsquerschnitt fiir die im Bogenplasma enthaltenen Atome C, N 
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und O erhalt. Diese Auswertung fiihrt zu dem Resultat: 


C,N,O — 99-40-16 cm? + 15%. (30) 


0 


Da diese drei Atomarten nebeneinander im Periodischen System legen, 
kann man versuchsweise den gefundenen Querschnitt jeder Atomart 
einzeln zuordnen, wovon im vorigen und im nachsten Kapitel Gebrauch 
gemacht worden ist. 


VI. RAMSAUER-Querschnitte fiir Fe-Atome 
aus Messungen am Pfundbogen. 


F. Burnorn [17] hat im Eisen-Pfundbogen, der zwischen FeO- 
Elektroden mit einer Stromstarke von 5 Amp bei einer Lange von 1,5 bis 
2mm brennt, Feldstarke, Elektronendichte und radiale Temperatur- 
verteilung gemessen. Die Feldstarke hat er aus der Anderung der Bogen- 
spannung mit der Langenanderung zu 40 V/cm bestimmt. Die Elek- 
tronendichte ergab sich aus der absoluten Intensitat des Kontinuums in 
der Achse — nach Umrechnung der beobachteten Querverteilung auf 
die radiale Intensitatsverteilung mit Hilfe der AxBeEitschen Integral- 
gleichung — gemaB der Beziehung [/8) 


txont = 6,36 hi 10-% n° | kT (31) 


zu 3,5 -410'6 cm~* bei einem maximalen IJonisationsgrad von 3%. Aus 
diesem Wert in Verbindung mit der SAHA-Gleichung und aus DoPpPLER- 
Breitenmessungen ergab sich eine Achsentemperatur von 6400° K. Die 
radiale Verteilung der Temperatur wurde aus dem relativen Intensitats- 
abfall geeigneter Spektrallinien gema&B Gl. (16) ermittelt. Damit sind 
uns alle Unterlagen gegeben, um auch am Fe-Bogen RAMSAUER-Quer- 
schnitte nach Gl. (12) zu ermitteln. Fiir die neutralen Atome im Pfund- 
bogen findet man so den Wert 


Qr.G = 27:107'% cm?, (32) 


Zerlegt man diesen Querschnict nach dem Vorbild von Gl. (28) in die 
Anteile des Fe und des O, die beide in gleicher Zahl im Bogenplasma 
vorhanden sind, und setzt wieder fiir Q° etwa 20 - 10-* cm? ein, dann ist 


ore = 35-10% cm? + 30%. (33) 


VII. Graphische Darstellung der RAMSAUER-Querschnitte. 

Zum SchluB dieser Abhandlung haben wir einmal die Bereiche, in 
denen die bisher bekannten RAMSAUER-Querschnitte von Atomen fiir 
den fraglichen Geschwindigkeitsbereich [79] liegen, gegen die Valenz- 
elektronenzahl aufgetragen (Fig. 13). Wie man sieht, haben die 
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Elemente gleicher Valenzelektronenzahl auch dieselbe GréBenordnung 
des RAMSAUER-Querschnitts, wobei nur H und He etwas herausfallen. 
Die Ubergangselemente (z.B. F e) lassen sich natiirlich nicht in dieses 
Schema einordnen, Vielleicht kénnte man in dieser groben Darstellung 
sogar einen gewissen Gang mit der Valenzelektronenzahl vermuten. 
In diesem Zusammenhang waren Messungen des RAMSAUER- Querschnitts 
an den Halogenen mit der Valenzelek- 


7 16 2 
7032 om 
tronenzahl 7 von groBem Interesse. 


Zusammenfassung. Alkalien 


Wie in jedem inhomogenen Leiter 


muB auch bei einer Bogenentladung der & jf rg net 

aus dem Verhaltnis Stromstarke zu %$ 

Feldstarke berechnete Leitwert pro cm 8 

ubereinstimmen mit der iiber den Quer- in apn 9 ph % 
schnitt integrierten Leitfihigkeit. Da Cama 

auf Grund kinetischer Vorstellungen iS 10 

die Leitfahigkeit wber die Wirkungs- 

querschnitte der bremsenden Ionen und ie 
Atome mit der Elektronendichte und ~ lke 
der Temperatur verbunden ist, kann 5000°K<T< 12000 °K re 
man aus der gemessenen Elektronen- 7 r 7 ; 3 
dichte- und Temperaturverteilung die Anzahl der Valenzelektroden 


Ionen und Atomquerschnitte ermitteln. Fig. 13. Bisher bekannte Ramsaver-Quer- 
Sains schnitte ber der Valenzelektronenzahl 
Messungen am wasserstabilisierten aufgetragen. 


Hochleistungsbogen (GERDIEN-Bogen), 
die wegen gewisser friiherer Unstimmigkeiten iberpriift wurden, 
fiihrten zu folgender Gleichung fiir den Ionenquerschnitt: 
P LAeT a Ng o/ 
¢. = a7 oad EY a _ 
Dieser Querschnitt stimmt zahlenmaBig weitgehend mit theoretischen 
Berechnungen von L. SPITZER iiberein. 
Untersuchungen am normalen wasserstabilisierten Bogen und am 


wasserstabilisierten Hochdruckbogen ergaben fiir den RAMSAUER-Quer- 
schnitt der H-Atome 


G*= 130-410 cm?’ 25%. 


Ahnliche Auswertungen am Hochstromkohlebogen brachten fiir die 
Elemente C, N, und O den Mittelwert 


BO =: 20.4054. om? 15%, 
10* 
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Die entsprechenden Messungen und Berechnungen fiir den Fe-Pfund- 
Bogen hatten den Wert 


8 bi = 35 F 409716 cm2 + 30% 
zum Ergebnis. 

Die Messungen am GERDIEN-Bogen, am normalen wasserstabilisierten 
Bogen und am wasserstabilisierten Hochdruckbogen wurden in Kiel 
durchgefiihrt, wahrend die Untersuchungen am Hochstromkohlebogen 
und am Eisenbogen in Erlangen angestellt wurden. 
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Lumineszenzuntersuchungen 
an Uran-aktivierten Erdalkaliwolframaten 
und -molybdaten. 


Von 
H. GOBRECHT und WERNER WEIss. 
Mit 4 Figuren im Text. 


(Eingegangen am 20. Oktober 1954.) 


Das verschiedenartige Auftreten und Verhalten der Lumineszenz Uran-aktivierter 

Erdalkaliwolframate und -molybdate wird untersucht. Man beobachtet zwei von- 

einander getrennte Lumineszenzbanden im grunen und im roten Spektralbereich, 

die durch chemische Zusatze zu den Wolframaten und Molybdaten vor dem Gliihen 

einzeln in ihrer Strahlungsstarke, die griine Bande auch in ihrer spektralen Lage, 
verandert werden k6énnen. 


Einleitung. 

Uber die an und fiir sich bekannte Lumineszenz von Uran-aktiviertem 
Calcium- und Strontiumwolframat durch Anregung mit langwelligem 
ultravioletten Licht sind in der Literatur beziiglich der Lumineszenz- 
farbe keine einheitlichen Angaben zu finden. So soll z.B. CaWO,-U 
nach KROGER | J) bei — 180° C rot und griin, nach GARLICK [2] nur griin 
lumineszieren. Der Grund dieses Verhaltens ist bereits von KROGER [1] 
in der verschiedenen Zusammensetzung der Ausgangssubstanz vermutet 
worden. 

In der vorliegenden Arbeit wird diese Annahme bestatigt. Es zeigen 
sich systematische Abhangigkeiten der spektralen Verteilung der Lu- 
mineszenz, ebenso der Strahlungsstarke der roten und griinen Bande von 
der Zusammensetzung der Grundsubstanzen. Ferner hangen die rela- 
tiven Strahlungsstaérken der beiden Banden in so verschiedener und 
charakteristischer Weise von den Gliihbedingungen, von der Aktivator- 
konzentration und von der Beobachtungstemperatur ab, dab mit hoher 
Wahrscheinlichkeit ein gemeinsamer Ursprung der roten und der griinen 
Bande ausgeschlossen werden kann. Die Versuche wurden nicht nur 
mit Ca— und SrWO,-U angestellt, sondern auch auf BaWO,-U und die 
entsprechenden Uran-aktivierten Erdalkalimolybdate ausgedehnt, deren 
Lumineszenzverhalten unseres Wissens noch nicht beschrieben worden ist. 


Apparatur zur Messung der Strahlungsstarke. 
Das Licht einer Quecksilberhochdrucklampe passiert einen Konden- 
sor und drei geeignete Filter zur Aussonderung der Linie 3650 A und 
dient dann zur Anregung. Die lumineszierende Leuchtstoffflache wird 
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in einer Irisblende abgebildet, die wahrend aller Messungen die gleiche 
Offnung beibehalt und einen Teil der Leuchtstoffflache ausblendet. Das 
durch sie hindurchgehende Licht gelangt auf die Kathode einer Photo- 
zelle mit Sekundirelektronenvervielfachung (Multiplier). Das ange- 
schlossene Galvanometer liegt in einer Kompensationsschaltung zum 
Ausgleich des unvermeidlichen Dunkelstromes des Multipliers. Zur 
Trennung des Lumineszenzlichtes dienen zwei Filterkombinationen: Fir 
das griine Leuchten die Scuorr-Filter VG 6, VG 8, BG 19; fiir das rote 
Leuchten RG 5 und BG 19, je 1 mm stark. 

Die Messungen erfolgten auch bei tiefer Temperatur. Dazu wurde 
die Probe vorher in fliissige Luft getaucht. Zur Messung der Abhangig- 
keit der Strahlungsstarke von der Temperatur wurde ein kleines Thermo- 
element in den Boden des Leuchtstoffbehalters eingeschmolzen und ge- 
eicht. Sollte die Probe bei hherer Temperatur untersucht werden, dann 
wurde sie erhitzt und die Messung wahrend der natiirlichen Abkiihlung 
durchgefiihrt. Die Messungen wurden je 5mal wiederholt. Der mittlere 
Fehler der Photostromablesung ergab sich zu 20 bis 30%. Wegen der 
Selektivitat des Empfangers kann in den Diagrammen die gemessene 
Strahlungsstarke der roten Bande nicht mit der der griinen Bande ver- 
glichen werden. Als Anhalt kann jedoch dienen, daB sich die Empfind- 
lichkeit der Anordnung fiir griin:rot etwa wie 10:1 verhalt, wie eine 
Empfindlichkeitsmessung ergab. Die Strahlungsstarke ist dem Photo- 
strom etwa proportional, trotz geringer spektraler Verschiebungen inner- 
halb der einzelnen Banden bei den verschiedenen Stoffen. 


Die Leuchtstoffpraparationen. 

Die Wolframate wurden sowohl durch Fallung als auch aus Oxyden 
hergestellt. Die Ausgangssubstanzen bei den Fallungsreaktionen waren 
einerseits NagWO, + 2H,0O p.a., andererseits Sr(CH,COO), +1/, H,O p.a., 
Ba(NOs)y p.a. und Ca(NOs), reinst. Molybdate wurden nicht gefallt. 
Das aus heiBen Lésungen gefallte Produkt wurde Smal mit heiBem bi- 
destilierten Wasser ausgewaschen, abfiltriert, getrocknet, in einem Ofen 
1 Std lang bei 1150°C im unglasierten Porzellantiegel gegliiht und dann in 
einer Reibschale fein zerstoBen. In dieser Form bildete es die Grund- 
substanz fiir alle weiteren Praiparationen. Es ist namlich unzweck- 
maBig, das gefallte, getrocknete Salz bereits fiir die Aktivierung und Sen- 
sibilisierung zu verwenden, da es undefinierte Wassermengen enthalt, 
die mitgewogen werden und die Konzentrationsangaben verfalschen. Da 
keine Kenntnisse dariiber bestehen, ob gefallte Wolframate eindeutige 
Zusammensetzung haben (3), wurde bei einigen Proben die Reihenfolge 
bei der Fallung gewechselt, d.h. es wurde sowohl z. B. Na,WO,-Lésung 
zu Ba(NOs).-Lésung gegeben als auch umgekehrt. Beide Produkte wur- 
den dann gleichartig weiterbehandelt und mit gleichen Mengen Uran 


ts ae Po 


Lumineszenzuntersuchungen, 141 


aktiviert. Da die Lumineszenz von Uran-aktivierten Wolframaten nach 
unseren Versuchen stark von der Zusammensetzung der Grundsubstanzen 
abhangt, hatten sich hierbei Unterschiede bemerkbar machen kénnen. 
Dies war aber nicht der Fall. 

Bei der Fallung bildet sich wegen der doppelten Umsetzung ein 
anderes, lésliches Salz, das der Niederschlag trotz mehrmaligen Wa- 
schens fest adsorbiert. Die Anwesenheit dieses Zusatzes bemerkt man 
auch aus der Sinterung beim Gliihen, die bei den aus Oxyden gewonnenen 
Verbindungen bei 1150° C kaum auftritt. Das verbleibende Salz nimmt 
einen entscheidenden Einflu8 auf die Lumineszenz des mit Uran-akti- 
vierten Erdalkaliwolframats. 

Um auch zusatzfreie Grundsubstanzen zur Verfiigung zu haben, 
wurden diese auch aus Oxyden bzw. Carbonaten hergestellt. Die Aus- 
gangssubstanzen waren H,WO,p.a. bzw. MoO, reinst und CaCO,, 
SrCO,, BaCO,. Die Erdalkalicarbonate lagen in besonders reiner Form 
vor!. Dies erwies sich spater als vorteilhaft, weil die Lumineszenzfarbe 
in den untersuchten Leuchtstoffen besonders stark von Zusdtzen ab- 
hangig ist. Die Oxyde wurden im richtigen Mengenverhiltnis in gla- 
sierten Porzellanschalen mit bidestilliertem Wasser gut vermischt, zu 
einem Brei angeriihrt, im Trockenschrank getrocknet und dann in un- 
glasierten Porzellantiegeln bei 1150° C 1 Std lang gegliiht. Die Gliihung 
wurde 2 bis 3mal unterbrochen und der Tiegelinhalt in einer Reibschale 
kraftig durchmischt. Bei den meisten Versuchen mit CaWO, ist statt 
eines selbsthergestellten Produktes zu Vergleichszwecken CaWO, von 
MERCK als Grundsubstanz verwendet worden. 

Die Gliihungen erfolgten bei allen Versuchen grundsatzlich in un- 
glasierten, offenen Porzellantiegeln. Der benutzte Ofen war ein ,,Silit‘’- 
Stab-Muffel-Ofen. Seine zulassige Maximaltemperatur betrug 1400°C. 
Die Temperaturmessung erfolgte mit einem Pt/Pt-Rh-Thermoclement. 

Das Uranpraparat war in einer Tropfenflasche als Uranylacetat- 
lésung vorhanden. Es eignen sich dazu aber auch alle anderen Uran- 
salze, die in der Hitze zu UO, zerfallen. Die Lésung wurde so angesetzt, 
daB in einem Tropfen 2,5-10°4g Uran vorhanden waren. Grund- 
substanz (gefallt oder aus Oxyden hergestellt) und Zusatze wurden nach 
Wagung mit der Uranldsung und bidestilliertem Wasser vermischt und 
wie oben angegeben behandelt. Die Wagegenauigkeit betrug etwa 1%. 


Der Verlauf der Strahlungsstarke beider Banden 
bei Anderung der chemischen Zusammensetzung. 

Zwei verschiedene Gruppen von Versuchsreihen wurden durchge- 
fiihrt : Einmal wurden zu den aus Oxyden hergestellten Grundsubstanzen 
verschiedene Salze verschiedener Konzentration gegeben und die 

1. Sie wurden uns freundlicherweise von der Auer-Gesellschaft zur Verfiigung 


gestellt. 
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Strahlungsstarke des Rot- und Griinleuchtens bei konstanter Anregung 
gemessen. Die Grundsubstanzen waren folgendermaBen erhalten worden: 


CaWO, aus CaWO, (Merck), 
SrWO, aus 1 mol SrCO; + 1 mol H,WO,, 
BaWO, aus 1,02 mol BaCO, + 1 mol H,WO,. 


Die Urankonzentration betrug stets 0,5 mol-%, die Gliihdauer 15 min 
bei 1150° C. Die Zusammensetzung der Leuchtstoffe bezieht sich stets 
. _ auf das Produkt vor dem Glihen. 

: *——— | | Vonden vielen Versuchsreihen sei 
ein Ergebnis als Beispiel heraus- 
gegriffen. Es ist in Fig.1 dar- 
100, Pf : | gestellt. Man sieht deutlich ein 
| Uberwiegen der roten Lumines- 
SPT ‘ zenz bei Zugabe von Na,S. Gibt 
2 1 man nun statt dieses Salzes z. B. 

Na,P,0, hinzu, dann erhalt man 
4 8 2 6mo/% ein Uberwiegen des griinen Leuch- 
; a hy EP ae tens. Die Tabelle 1 zeigt iiber- 
ae eee OC von dew’ Gehalt aa Nag sichtlich, welchen Einflu8 eine 
(vor dem Glihen) in: J CaWO,-U; 2SrWO,U;  Reihe von Salzen auf die Lumi- 
bb ti Pnkeaininasar teh es, pee ee neszenzfarbe der Uran-aktivierten 

Erdalkaliwolframate haben. 

Die andere Gruppe von Versuchsreihen war folgender Art: Es wurden 
Erdalkalicarbonate der Grundsubstanz in wachsender Menge zugesetzt. 
In Fig. 2 sind die MeBreihen fiir SrcWO,-SrCO,-U und CaWO,-CaCO,—-U 
dargestellt. Bei der Wahl des zusiatzlichen Erdalkalicarbonats ist es 
gleichgiiltig, ob die gleiche Art der bereits vorhandenen Erdalkaliionen 
oder eine fremde Art hinzugefiigt wird, z. B. zu SrWO,-U BaO usw.: Der 
Verlauf der Kurven ist stets ahnlich. Die Abweichungen werden dann 
groBer, wenn auBer einem Erdalkalioxyd oder -carbonat noch ein sensi- 
bilisierendes Salz fiir die eine oder andere Farbe hinzugegeben wird, 


S 
S 


10 


Strahlungsstarke(willkiirlicheLinhe/t) 


oa 


Tabelle 1. 
Na SO, (R) Na,HPO, (G) | KNO, (R+G) | Na,SO, (—) 
NaHSO, (R) Na,P,0, (G) | NaNO, (R+G) | Na,CO, (—) 
LigSO, (R) Na,B,O, (G) | CH,;COONa (R+G) | NaOH (—) 
Na,S (R) | BeSO, (G) | (CH,COO),Sr (R+G) | NaCl (—) 
Mgso, (R) Al,(SQ,)3 (G) | Ba(NO,), (R+G) | KCl (—) 
MgBry (R) | Na,WoO, (G) | CaF, (R+G) | Naf (—) 
CaS (RR) | CaSO, (R+ G) | CsCl (— 
BaS (R) | SrSOQ, (R+G) | (CH,;COO),Pb (—) 
| | BaSO, (R+G) | 


Die als Zusitze zu den Erdalkaliwolframaten verwendeten Salze. (R) bedeutet: 
rot sensibilisierend; (G) = griin sensibilisierend; (—) nicht sensibilisierend. 
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da dann die Anionen und Kationen dieser Salze eine zusiitzliche Wirkung 
auf die Lumineszenzfarbe ausiiben. 

Man hat somit, um griin leuchtende Uran-aktivierte Erdalkali- 
wolframate zu erhalten, einen hohen Uberschu8 an Erdalkalioxyd bzw. 
-carbonat zu geben. Um rot leuchtende zu erhalten, miissen ein geringer 
UberschuB an Erdalkalioxyd und ein rot sensibilisierendes Salz hinzu- 
gefiigt werden. Das gleiche gilt fiir die entsprechenden Molybdate. 


Strahlungsstorke (willkiirliche Ginhert) 


e@, —— 


mole 6 202 70 i 18 22 26 30 34mol% 
——4We ClO, bow. Sr C05-—~ 


Fig. 2. Abhangigkeit der Strahlungsstarke der Lumineszenz bei — 180°C von der Zusammensetzung bei 
1 CaWO,-U und 2 SrWO,-U. Rote Lumineszenz ausgezogen; griine Lumineszenz gestrichelt. 
Auf Grund dieser Versuche kénnen fiir die Préparation von Uran- 
aktivierten Erdalkaliwolframaten und -molybdaten, die in verschiedenen 
Lumineszenzfarben optimal leuchten, genaue Angaben gemacht werden, 
wie sie in Tabelle 2 zusammengestellt sind. 


Tabelle 2. Zusammenstellung giinstigster Pyrdparationen einiger Uvan-aktivierter 
Evrdalkaliwolframate und -molybdate. 


CaWo,-U 
Tiefgriin (20°, — 180°): 
CaWO, + 40 Mol-% CaCO, + 0,05 Mol-% U 30min 1200°C 


Rot (— 180°): 
CaWO, + 3 Mol-% CaCO, + 0,1 Mol-% U + 6Mol-% Na,S 60min 1100°C 


SrWO,-U 
Griin (20°, — 180°): 
SrWO, + 40 Mol-% SrCO, + 0,05 Mol-% U 30min 1200°C 


Rot (20°, — 180°): 
SrWO, + 4 Mol-% SrCO, + 0,1 Mol-% U + 6Mol-% Na,S 60min 1100°C 


BaWO,-U 
Griingelb (20°, — 180°): : 
BaWO, + 40 Mol-% BaCO, + 0,05 Mol-% U 30 min 1200°C 


Rot (20°, — 180°): an . 

BaWO, + 4 Mol-% BaCO, + 0,1 Mol-% U +6Mol-% Na,S 60min 1100°C 
CaMoO,-U 

Rot (— 180°): , : 

CaMoO, + 2 Mol-% CaCO,+ 0,1 Mol-% U +6Mol-% Na,S 60min = 1200 Cc 
SrMoO,-U 

Rot (20°, — 180°): < 

SrMoO, + 4 Mol-% SrCO,+ 0,1 Mol-% U +6Mol-% Na,S 60min 1200 c 
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EinfluB der Urankonzentration. 

Eine groBe Zahl verschiedener Praparationen dieser Uran-aktivierten 
Erdalkaliwolframate hat tibereinstimmend gezeigt, daB die maximale 
Strahlungsstarke der roten Bande bei einer Urankonzentration von etwa 
0,1 Mol-% auftritt, wahrend die maximale Strahlungsstarke der griinen 
Bande nur eine Urankonzentration von 0,02 bis 0,05 Mol-% bendtigt. 
Dies gilt unabhangig von der Zusammensetzung und der Gliihbehandlung 
der Leuchtstoffe. Gleiche Verhaltnisse ergeben sich auch bei den Uran- 
aktivierten Ca- und Sr-Molybdaten. Dieses deutet auf einen verschieden- 
artigen Einbau des Urans hin. Wenn das Uranleuchtzentrum fiir die 
rote und die griine Lumineszenz den gleichen Platz einnehmen wiirde, 
dann ware der verschiedenartige Einflu8 der Urankonzentration auf die 
Strahlungsstarke beider Banden nicht einzusehen. 


Anderung der Gliihtemperatur. 

Die Versuche haben weiter gezeigt, daB bei diesen Leuchtstoffen fiir 
die Erzeugung einer maximalen Strahlungsstarke der roten Bande eine 
um 100 bis 300° héhere Gliihtemperatur erforderlich ist als zur Erzeu- 
gung der maximalen Strahlungsstarke der griinen Bande. Die giinstig- 
sten Temperaturen liegen — bei einer Gliihdauer von 1 Std — bei 
1150 bis 1250° C fiir rotes, und bei 950 bis 1150° C fiir griines Leuchten. 


Temperaturabhangigheit der Lumineszenz. 
Fig. 3a und b zeigen als Beispiele den Verlauf der relativen Strah- 
lungsstarke der roten und griinen Bande zweier verschieden zusammen- 
gesetzter Uran-aktivierter Erdalkaliwolframate in Abhangigkeit von 
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Fig, 3a—d. Abhangigkeit der relativen Strahlungsstirke der roten (R) bzw. griinen (G) Lumineszenz 
von der Temperatur. 

der Temperatur, wobei die maximale Strahlungsstirke willkiirlich je 

gleich 1000 gesetzt ist. Bei allen hier untersuchten Leuchtstoffen nimmt 

die Strahlungsstirke mit wachsender Temperatur ab. Charakteristisch 


i 
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jedoch ist der verschiedenartige Verlauf der Abnahme beider Banden: 
Wahrend dieser fiir die griine Bande bei allen Uran-aktivierten Erdalkali- 
wolframaten nahezu der gleiche ist, ist der Verlauf fiir die rote Bande 
bestimmt durch die Art des Grundgitters. Der Abfall der roten Strah- 
lungsstarke ist bei CaWO,-U am starksten, bei BaWO,-U am schwachsten. 
Die Kurve fiir SrWO,-U geht etwa parallel zur _,,griin‘“-Kurve. Die 
Fig. 3c und d zeigen die Kurven zweier entsprechender Molybdate. Auch 
hier ist die ,,griin“-Kurve unverandert, wahrend die ,,rot‘““-Kurve vom 
Kation der Grundsubstanz bestimmt wird. Beide Kurven sind in ihrem 
Verlauf stets unabhangig von der Art des zugefiigten Erdalkaliions. Man 
sieht also, daB die Temperaturabhangigkeit der griimen Bande vom 
Grundgitter unbeeinfluBt ist im Gegensatz zur Temperaturabhingigkeit 
der roten Bande. 


Abklingmessungen an SrWO,-U. 


Anregung mit Licht der Wellenlange 3650 A ergab im Phosphoroskop 
bei einer Anregzeit von 16-10 ~4sec exponentielles Abklingen beider 
Banden. Abklingmessungen bei Elektronenanregung ergaben ebenfalls 
exponentielle Abklingung. Als Abklingkonstanten wurden fiir beide 
Banden fir Ultraviolettanregung 10*sec 4, bei Elektronenanregung 
3 -10*sec! gefunden mit einem Fehler von etwa 10%. 


Lettfahigkettsuntersuchungen. 

Hier ergab sich an vorwiegend rot und vorwiegend griin leuchtenden 
Proben kein Einflu8 des Uranaktivators auf den Photostrom. Der 
Photostrom rot und griin leuchtender Proben zeigte keinen Unterschied. 
Naheres iiber Versuchsanordnung und Ergebnisse findet sich in der Arbeit 
von GOBRECHT, HAHN und KOsEL (4). 


Anregharkett der Uranlumineszenz. 


Versuche iiber die Anregbarkeit der Uranlumineszenz in Uran-akti- 
vierten Erdalkaliwolframaten und -molybdaten zeigten verschiedene 
Ergebnisse fiir die rote und griine Bande. Das spektral zerlegte Licht 
einer Hg-Hochdrucklampe erregte den Leuchtstoff, dessen Lumineszenz 
mit einer der obengenannten Filterkombinationen betrachtet wurde. 
Dabei ergab sich folgendes: Die rote und die griine Emissionsbande lassen 
sich im Gebiet von 265 bis 391 my gut anregen. Dabei regen die Linien 
Hg 365 und Hg 391 mu beide Banden besonders stark an. Die rote Bande 
allein lat sich dariiber hinaus noch von der Linie Hg 405 my gut an- 
regen, die griine dagegen nicht. 

Elektronen von einigen Tausend Elektronenvolt erregen beide Banden 
gleichmaBig kraftig. Eine Anregbarkeit durch Réntgenstrahlung von 
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51) kV besteht nicht. Werschiedenheiten der Reflexionsfarbe im Sicht- 
baren sind nicht vorhanden: Alle Substanzen sind (bei héchstens 0,1 Mol- 


Urankonzentration) rein well. 


Die spektrale Verteilung der Lumineszenz 
Besonders auffallend ist die Veranderung der spektralen Lage der 
Emission bei Variation der Zusammensetzung. Um beide Banden gut 
photographieren zu konnen, 
= wurden diese Versuche nur 
- - : mit Wolframaten durchge- 


2 fiihrt, welche stark (8 Mol-%o 
TE 


= Sere mit Erdalkalicarbonat ver- 


setzt waren. Fig. 4 zeigt neun 
uy 1 4 


as Spektralaufnahmen Man 


sieht deutlich, daB sich bei 


, erschiedener Erdalkalicar 
ie nate uberhaupt nicht ver 
ae schiebt. Die griine Bande ver- 

schiebt sich bei CaWO,-U 

nicht. Bei SrWO,-U_ ver 

schiebt sie sich bei Ersatz 

Ss ene von SrCOQ, durch CaCQ, 
‘ kiirzeren Wellen, beim 

BaWO, l bei Ersatz von 

; BaCQ, durch SrCO, eben 

wie) falls und durch CaCO, noch 

. weiter nach kiirzeren Wellen 

. Man sieht also, daB eine Vet 
i ax.enk tani schiebung der griinen Bande 


nur dann erfolgt, wenn ein 
Erdalkalion hinzugefiigt wird, welches kleiner ist als das Erdalkaliion 
des Woltramats. Nur das kleinere Erdalkaliion bestimmt die Lage der 
ertinen Uranlumineszenz, unabhangig vom Grundgitter. In der Tat be- 
wirken nicht nur Erdalkalicarbonate und -oxyde diese Verschiebungen 
der griinen Bande, sondern alle Salze, die ein entsprechendes Erd- 


alkaluon besitzen, z.B. CaF,, Ca(NO,.),, CaS usw. Dabei kénnen die 


Zusatzkonzentrationen schwanken zwischen 2 und 30 Mol-%. Diese 


Ergebnisse sind unabhaingig von Urankonzentration und Gliihbe 
dingungen 


Auch zusitzliche Alkalionen ergeben eine charakteristische Beein 


flussung der griinen Bande. Besonders auffallig ist in diesem Fall das 
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Auftreten typischer Uranylbanden, wie man sie z.B. von der Lumines- 
zenz von Uranylsalzen her kennt. 


Um die Struktur der roten Bande zu erkennen, wurden die Spektral- 
aufnahmen photometriert. Wahrend bei CaWO,-U, BaWO,-U und 
CaMoO,-U mit verschiedenen Erdalkalizugaben keine Struktur zu er- 
kennen war, bemerkt man bei SrWO,-U und SrMoO,-U bei tiefer Tem- 
peratur eine deutliche Strukturierung der roten Bande. 


Die Verschiebung der griinen Bande durch Einwirkung von Fremd- 
ionen ist nicht auf Gitteranderungen zuriickzufiihren: DEBYE-SCHERRER- 
Aufnahmen an diesen Praparaten haben keine Anderung der Gitter- 
konstanten iiber 0,03% durch Zugabe von 8 Mol-% Erdalkalioxyd an- 
gezeigt. (Die mittlere Anderung der Gitterkonstanten bei Ubergang von 
CaWO, zu BaWO, betragt etwa 4%.) Es tritt auch keine Anderung der 
roten Bande ein, was ebenfalls Gitteranderungen ausschlieBt. McKEAG 
und Ransy [5] fanden allerdings réntgenographisch Gitterdnderungen 
bei Zugaben von SrO und MgO zu CaWO,, jedoch ohne die Oxydkonzen- 
trationen angegeben zu haben. Bei zu hoher Oxydkonzentration wird 
sich natiirlich etwa aus 


CaWO, + SrO = SrWO, + CaO 


bilden und damit auch die Lage der roten Bande verschieben, was unsere 
Versuche auch bestatigt haben. 


Diskussion der Ergebnisse. 


Abkling-, Leitfahigkeitsmessungen und die zum Teil strukturiert auf- 
tretenden Banden machen es wahrscheinlich, daB es sich bei der Lumi- 
neszenz beider Banden um ,,monomolekulare“ Vorgange handelt, der 
Ursprung der Lumineszenz also wohl in geschlossenen Komplexen des 
Urans zu suchen ist. Als solche Komplexe kommen nur die des 6-wertigen 
Urans, also das Uranyl-Ion UO;* und das Uranat-Ion UO, ’, in Frage. 
Beide unterscheiden sich voneinander in der ,,Umgebung“‘ des Uran- 
atoms. 

Da reine Uranylsalze in ihrer griinen Fluoreszenz der griinfluores- 
zierenden Varietat der hier behandelten Stoffe sehr ahnlich sind, wird 
man auch in den Uran-aktivierten Wolframaten die UO} *-Gruppe fiir 
die griine Lumineszenz verantwortlich machen diirfen. Hingegen hat 
die rote Lumineszenz dieser Stoffe médglicherweise in den Uranatkom- 
plexen ihren Ursprung. Diese Vermutung gibt auch Krocer [J] an. 
Ein Beweis fiir die Richtigkeit dieser Vermutung ist aber nicht vorhanden. 

Wiahrend die griine Lumineszenz in vielen mit Uran aktivierten 
Stoffen und in reinen Uranylsalzen beobachtet wird, existiert die rote 
Lumineszenz — soweit bis heute, einschlieBlich dieser Arbeit, bekannt — 
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nur in den Uran-aktivierten Wolframaten der Metalle der zweiten Haupt- 
und Nebengruppe einschlieBlich des BaWO, und in den Molybdaten von 
Ca und Sr. Auch in Alkali- und Erdalkaliuranaten ist sie nicht zu finden. 
Es scheint also die Existenz der roten Bande an das Wolframat- bzw. 
Molybdatgitter der Scheelitform gebunden zu sein. Die Lumineszenz 
von Zn-, Cd-, Mg-Wolframat, die im Wolframitgitter kristallisieren, hat 
andere Eigenschaften, auf die hier nicht eingegangen wird. 

Dieser Sachverhalt wird durch unsere Versuche bestatigt: Das 
Temperaturverhalten der Strahlungsstarke, die Spektrailage und die 
Bandenstruktur der roten Lumineszenz werden nur durch das Grund- 
gitter, nicht aber durch irgendwelche Zusatze zum Grundgitter bestimmt. 

Demgegeniiber ist die Existenz der griinen Lumineszenz nicht an 
das Grundgitter gebunden. Das zeigen die Unabhangigkeit des Tempe- 
raturverhaltens der Strahlungsstéarke vom Grundgitter sowie der Spek- 
trallage und der Bandenstruktur vom Grundgitter, wenn diesem be- 
stimmte Kationenzusatze hinzugefiigt werden. 

Auf Grund dieses verschiedenen Verhaltens der beiden Banden k6n- 
nen Aussagen iiber die nahere ,,Umgebung‘ des Aktivators gemacht 
werden: Da namlich die griine Bande in ihrer spektralen Verteilung stark 
von Kationenzugaben beeinfluBt wird, die rote dagegen aber véllig un- 
empfindlich ist, kbnnte man das ,,griine“‘ Zentrum in der ,,Nahe“‘ eines 
hochgeladenen Kations, das ,,rote“’ — gut gegen elektrostatische Krafte 
abgeschirmte — etwa im WO, bzw. MoO,” — Tetraeder lokalisieren. 
Eine verschieden groBe Einbauwahrscheinlichkeit beider Zentren liegt 
offensichtlich vor, wie man aus dem verschiedenen EinfluB der Uran- 
konzentration und der Gliihtemperatur auf die Strahlungsstarke beider 
Banden entnehmen kann. Die hdhere Einbauwahrscheinlichkeit hat 
hiernach das ,,griine“‘, die geringere das ,,rote‘‘ Zentrum. Man kénnte 
sogar zufolge dieser — allerdings recht geringen — Unterschiede an- 
nahmen, daB das ,,rote‘‘ Zentrum isomorph, das griine aber in Zwischen- 
gitterplatzen angeordnet ist. Die starke Wechselwirkung des roten 
Zentrums mit dem Gitter, die besonders gut die Temperaturcharakteristik 
der Strahlungsstarke zeigt, befiirwortet iibrigens diese Deutung. 

Da Zugaben von Fremdkationen, die das griine Leuchten in seiner 
spektralen Verteilung beeinflussen, nach unseren Versuchen keine Ande- 
rung der Gitterdimensionen zur Folge haben, ist die obige Annahme, 
dem griinen Zentrum einen Zwischengitterplatz zuzuordnen, nicht ab- 
wegig. Sie geht mit der beobachteten Temperaturcharakteristik der 
Strahlungsstirke konform. Das beeinflussende Kation kann dabei offen- 
bar isomorph (z.B. in CaWO,-BaO-U) oder in Zwischengitterplatzen 
(z.B. in BaWO,-CaO-U) angeordnet sein. MaBgebend fiir die Beein- 
flussung ist wahrscheinlich nur das Ladungsdichteverhaltnis von Uber- 
schuBkation zu Grundgitterkation, 
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Der verschiedene Charakter beider Banden macht sich auch in der 
verschiedenartigen Sensibilisierbarkeit bemerkbar. Wahrend fiir die 
griine Lumineszenzintensitat insbesondere Erdalkaliionen (unabhangig 
von der Anionenart) sensibilisierend wirken, wird die rote Intensitit 
durch Zugabe gewisserAnionen (O--, S~-, SO;~ usw.) (unabhangig von 
der Kationenart) verstarkt. 

Da das ,,griine‘‘ Zentrum als in der ,,Nihe‘‘ der Kationen hegend 
gedeutet wurde, ist die Verstaérkung des Leuchtens mit wachsender Uber- 
schuBkationenzahl verstandlich. Der gleiche Effekt schlieBt einerseits 
nach dem Massenwirkungsgesetz den isomorphen Einbau der Urany]l- 
(UO;")Komplexe in das Gitter, also die Bildung von UO,WO,, aus. 
Andererseits kénnte mit der Méglichkeit der Bildung von Erdalkali- 
uranat gerechnet werden, was wiederum die (Uranat) UO; -Komplexe 
als ,,griine“’ Leuchtzentren markieren wiirde. 

Ahnlich wird sich die Sensibilisierung der gittergebundenen roten 
Lumineszenz durch Anionenzugaben kaum deuten lassen, wenn man 
nicht gerade annimmt, daB diese Anionen die Grundgitterplaitze der 
WO, -Ionen einnehmen. Vielleicht aber beruht dieser Sensibilisierungs- 
effekt zum Teil nur auf einem Wegfangen der fiir die griine Lumineszenz 
notwendigen Erdalkalionen durch die Zugabeanionen, so daB die rote 
Lumineszenzintensitat teilweise nur auf Kosten der griinen anwachsen 
kann. In der Tat scheinen unsere Ergebnisse (Fig.1) dies zu bestatigen. 
Daneben kann allerdings noch ein anderer Vorgang — etwa eine Ko- 
aktivierung — stattfinden, der das Entstehen der roten Lumineszenz 


begiinstigt. 
Herrn Dr. D. HAHN danken wir fiir zahlreiche Diskussionen, der 


Deutschen Forschungsgemeinschaft fiir finanzielle Unterstiitzung dieser 
Arbeit. 


Literatur. 
[1] Krécer, F.A.: Some Aspects of the Luminescence of Solids, S. 151. 
New York, Amsterdam u. London 1948. — [2] Garticx, G. F. J.: Luminescent 


materials. Oxford 1949. — [3] GmELIn: W-Bd., S. 289ff. — [4] Gosrecut, H., 
D. Haun u. H.-J. KésEx: Z. Physik 136, 57 (1953). — [4] McKeaa, Hy tit ba We 
Ranspy: USA. Pat. Spec. 2, 257, 699. 


Berlin-Charlottenburg, 11. Physikalisches Institut der Technischen 
Universitat. 


Zeitschrift fiir Physik, Bd. 140, S. 150—155 (1955). 


Zur Theorie der inneren Feldemission. 
Von 
F. SAUTER und J. WEISSE. 
Mit 1 Figur im Text. 
(Eingegangen am 5. November 1954.) 


Es wird eine besonders iibersichtlich erscheinende Methode zur Ableitung der ZENER- 
Formel fiir die innere Feldemission angegeben, welche im PErERLSschen Grenzfall 
schmaler verbotener Energiezonen im Material anwendbar ist. 


Unter der inneren Feldemission in einem Kristall versteht man _be- 
kanntlich den Ubergang eines Elektrons aus dem Valenzband in das 
Leitungsband des Kristalls unter der Wirkung eines von auBen an das 
Material angelegten elektrischen Feldes. Dieser Vorgang kann nach 
CL. ZENER}, in Analogie zur 4uBeren Feldemission an einem Metallrand, 
als ein wellenmechanischer Tunneleffekt durch ein fiir die Elektronen 
verbotenes Gebiet aufgefaBt werden, namlich durch die im Feld gekippte 
verbotene Zone zwischen Valenz- und Leitungsband. Zur Ermittlung 
der (nicht streng berechenbaren) Wahrscheinlichkeit fiir diesen Elek- 
tronendurchtritt wurden verschiedene Naherungsverfahren herangezogen, 
die im wesentlichen stets zu den gleichen Formeln fiihrten, so daB man 
das Problem der inneren Feldemission heute als im wesentlichen geklart 
ansehen kann®. Wenn hier dennoch erneut auf dieses Problem einge- 
gangen wird, so geschieht dies, um eine bisher anscheinend noch nicht 
benutzte Methode mitzuteilen, die in einem bestimmten Grenzfall eine 
wohl besonders einfache und wibersichtliche Behandlung der inneren 
Feldemission gestattet und die unter Umstinden auch bei anderen 
ahnlichen Problemen angewandt werden kann, 


Die Hauptaufgabe bei der Berechnung der Durchtrittswahrschein- 
lichkeit besteht im Aufsuchen einer méglichst zweckmaBigen Naherungs- 
lésung der fiir das Problem zustandigen ScHRODINGER-Gleichung 


ne . ee h op 

- 2m A+U (t) + ey ( t)} p mE — ine , (1) 
in der U(r) das Gitterpotential und § die als raumlich und zeitlich 
konstant angenommene duBere Feldstarke bedeutet; die Elektronen- 
ladung wurde dabei mit — ey bezeichnet. 


+ ZENER, Ct.: Proc. Roy. Soc. Lond. A 145, 523 (1934). 
* Vel. den zusammenfassenden Bericht von W. FRANz in den Erg. exakt. 
Naturwiss. 27 (1953), insbes. S. 13—20. 
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Im feldfreien Fall (%} = 0) kann man zur angenaherten Lésung von (1), 
je nach der Starke der Potentialschwankungen im Gitter, entweder nach 
PEIERLS in erster Naherung mit freien Elektronen rechnen und die 
Potentialschwankungen im Gitter als Stérung behandeln, oder aber 
nach BLocH von den an die einzelnen Gitterbausteine gebundenen Elek- 
tronen ausgehen. 

Auch im Fall §=+ 0 kénnte man zunichst an ahnliche Naherungs- 
verfahren denken. Doch la8t sich nur im Fall schwacher Potential- 
schwankungen im Gitter ein der PEreRtsschen Behandlung nachgebil- 
detes Verfahren durchfiihren, indem man in erster Naherung mit Wellen- 
funktionen rechnet, welche das Verhalten von Elektronen allein im 
#-Feld beschreiben, und das Gitterpotential U(r) als Stérung betrachtet ; 
dieses Verfahren wird im folgenden dargestellt. Eine entsprechende 
Ubertragung der Brocuschen Behandlung scheitert an der Unkenntnis 
der Wellenfunktionen von an Gitterionen gebundenen Elektronen im 
o-Feld. 

Liegt also der PEIERLSsche Grenzfall schwacher Potentialschwan- 
kungen vor, so ist es zweckmaBig, in erster Naherung mit den von 
C. G. DARwrin! angegebenen Impulseigenfunktionen 


we (t, t) = ay exp fi i (R(t) sh fo R2(r (2) 
zu rechnen, wobei zur Abkiirzung 
() =1— GF 3) 


gesetzt wurde. Ersichtlich stellen diese Funktionen Wellenvorgange dar, 
bei denen sich die Impulseigenwerte =A durch das Feld im Sinn 
der klassischen Mechanik zeitlich andern. 

Durch das Gitterfeld werden die reinen Zustande (2) gestort und 
gehen dadurch in eine Uberlagerung aus solchen Zustanden tiber, be- 
schrieben durch Wellenfunktionen der Form 


p(t, t) = [ aKea(t jexp fi (ae i (St( a hfe an). (4) 


x ‘1 Darwin, C. G.: Proc. Roy. Soc. Lond. A 154, 61 (1936). Zu véllig analogen 
Funktionen, namlich 


we (t, t) = by ox ih fe )as| (3') 


kommt man als Energie- und Impulseigenfunktionen, wenn man das ¥-Feld nicht, 
wie in (1), durch ein skalares Potential, sondern durch ein Vektorpotential YU = — cht 
beschreibt. Dabei gehen die Funktionen (3’) aus den Funktionen (3) durch eine 
Eichtransformation hervor. Vgl. hierzu F. SAUTER, Ann. Physik (5) 43, 404 (1943). 


Zeitschrift fiir Physik. Bd. 140. 
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Dabei besteht fiir die Amplituden a,(¢) wegen (1) die noch strenggiiltige 
Gleichung 


£24 Ya snl jexp| 2 fixe ~(8(e)—208)*} 47}, 6) 


u 


wenn man die aut Potentialfunktion U(r) durch die Fou- 
RIER-Reihe 
U(t) = >) Uexp {2 x1 (6; r)} (6) 
1 
darstellt und unter 6; die Ortsvektoren der Netzpunkte im reziproken 
Gitter versteht. 

Nach (5) induziert also die Gitterstérung Ubergange zwischen den 
einzelnen Elektronenzustanden, die auch im auBeren Feld genau den 
LavgEschen Interferenzbedingungen im Kristallgitter entsprechen!. Geht 
man namlich von einem Anfangszustand mit dem scharfen Ausbreitungs- 
vektor (0) =f, aus, so treten in der Streuwelle erster Naherung nach 
(5) nur Elektronenzustande auf, deren Ausbreitungsvektor 8 (¢) sich vom 
Ausbreitungsvektor St, (¢) der Primarwelle gerade um das 22-fache eines 
reziproken Gittervektors b; unterscheidet, fiir die also wegen (3) 


R(t) — Ry (t) =t—f, = 22}, (7) 


gilt. Auf die Impulsbilanz bei den Elektroneninterferenzen hat somit 
die durch das Feld bedingte Driftbewegung im &-Raum keinen EinfluB, 

Anders bei der Energie! Ohne Feld wird der durch (7) charakteri- 
sierte Interferenziibergang im allgemeinen wegen der Verschiedenheit 
der zu f, und zu f gehérigen Energien nicht méglich sein. Im Feld andern 
sich aber diese Energien in verschiedener Weise, so daB die Méglichkeit 
eines Gleichwerdens derselben in einem bestimmten Zeitpunkt durchaus 
gegeben ist. Und zwar folgt aus 


RS (¢) = R2 (2) (8) 


1 In der Metallelektronentheorie ist es vielfach iiblich geworden, statt von 
LAvE-Interferenzen von BraGeschen Reflexionen zu sprechen, auch wenn die 
Interferenzbedingung nicht in der von BraGcG gegebenen Form 2d sin #/2=nd 
vorliegt, sondern in der auf M. v. Lave zuriickgehenden, mit (7) identischen Form 
(t — f, a;) =22n,. Es bedeutet sicherlich keine Schmalerung der allgemein an- 
erkannten Verdienste BRaGGs um die Klarung der Interferenzerscheinungen an 
Kristallen, wenn hier dafiir eingetreten wird, daB auch in der Metallelektronen- 
theorie von Laus-Interferenzen gesprochen wird, sobald die ganze Blickrichtung 
mehr den Interferenzbedingungen (7) als der Vorstellung einer Reflexion an Netz- 
ebenen entspricht. —- In diesem Sinn der Wahrung historischer Prioritat miiBte 
beispielsweise in der von dem einen von uns angeregten und betreuten Dissertation 
des Herrn H. GERSTENKORN [Ann. Physik 10, 49, 80 (1952)] iiber den elektrischen 
Widerstand von Metallen in der Nahe des Schmelzpunktes statt von BraGc-Re- 


flexionen und BraGcG-Faktoren meist von Laur-Interferenzen und Lave-Faktoren 
die Rede sein. 
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mit (7) und (3) fiir diesen Zeitpunkt 
é: me h 

eo (¥ b;) 
Liegt er innerhalb der Beobachtungszeit, so wird man mit dem Auftreten 
. der durch den Index) gekennzeichneten Interferenz zu rechnen haben. 

Man kann sich die Verhaltnisse an einer zweidimensionalen Dar- 
stellung des &-Raumes klarmachen (vgl. Fig. 4 mit den fiir ein quadrati- 
sches Flachengitter eingetragenen BriLLoutnschen Zonen). Geht man 
von dem Ausbreitungsvektor §,(0) =f eines bestimmten Elektronen- 
zustandes zur Zeit ¢=0 aus und verfolgt die geradlinige Wanderung 


(b;, & + 2 b,). (9) 


des Endpunktes des Vektors 8, (¢) mit der Zeit, so trifft man zu einer 
bestimmten Zeit ¢; auf eine der BriLLourNschen Grenzlinien, die etwa 
dem Vektor b, entsprechen mége; dabei sind ¢, und b, miteinander durch 
die Beziehung (9) verkniipft. Beim Durchwandern von &, durch diese 
Grenzlinie besteht die Méglichkeit einer zum Vektor 6, nach (7) gehérigen 
Interferenzreflexion. Diese Reflexion, charakterisiert durch eine ge- 
strichelte Linie in Fig.1, wird mit einer bestimmten, unten berechneten 
Wahrscheinlichkeit eintreten. Der nicht reflektierte Anteil der Primar- 
welle, oder besser gesagt der Endpunkt des S,-Vektors trifft beim Weiter- 
laufen auf eine weitere BRILLOUINsche Grenzlinie, an der wiederum eine 
Interferenzreflexion méglich ist, usw. Durch das %-Feld werden eben 
die Wellenlangen der einzelnen Elektronenzustande in solcher Weise 
stetig verandert, daB ab und zu eine Lave-Interferenz bzw. eine BRAGG- 
Reflexion eintreten kann. 

Zur Berechnung der Reflexionswahrscheinlichkeiten aus der Gl. (5) 
hat man als Anfangsbedingung etwa a, ({ =0) = 6(f— f,) zusetzen. Dann 
kommt man bei Vernachlassigung der zwei- und mehrfachen Reflexions- 
prozesse zur Beziehung 


h Galt 
ava ote = >) Uo(t—t, — 20h) x | 
7 


x exp{ 24!” (by, fy+ 2b) ¢— “2% (8, 6) ar 
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die sich unmittelbar integrieren laBt. Dabei treten die aus der Optik 
bekannten FrEsNELschen Integrale C und S auf, die durch 


C(s) + 12S (s) = J doexp (i072 (14) 


definiert sind und sich in Form der Cornv-Spirale darstellen lassen. Da C 
und S als Funktionen von s fiir |s| >41 um die Konvergenzpunkte + } 


und — 4 oszillieren, schwankt das Integral [ do exp {t2107/2} in Abhan- 


gigkeit von den Grenzen s, und s, nur relativ wenig und besitzt kleine 
Werte, solange s, und sy das gleiche Vorzeichen und nicht zu kleine Be- 
triage haben; es erhalt jedoch angenahert den Wert 1+ 7 fiir —s,>1, 
Sy>1. Bei den aus (10) folgenden Zeitintegralen liegt der dem Punkt 
s=0O entsprechende Zeitpunkt, in dem die Phasenfunktion des Zeit- 
integrals stationar wird, im j-ten Summenglied genau bei dem durch (9) 
gegebenen Wert ¢;. Liegt also dieser Zeitpunkt auBerhalb des zeitlichen 
Integrationsintervalles, so fiihrt die Zeitintegration nach (10) nur zu 
einer schwachen, schnellen Oszillation von a, .9.»,(¢), entsprechend den 
bei jeder nichtstationaéren Stérungsrechnung auftretenden Schwan- 
kungen. Fallt aber ¢; in die Beobachtungszeit hinein, so wachst 4y, . 25, (#) 
in der unmittelbaren a dieses Zeitpunktes dem Betrag nach auf 


(//2/%) |G: |/m/2e0 | ( (b,)| an. Dabei ist die Wurzel durch die Trans- 
formation des ¢- Titegtals in die Form (11) bedingt und stellt (mit ihrer 
GréBenordnung von 10°! bis 10sec) ungefahr die Zeitdauer dar, 
innerhalb welcher der Integralwert wesentlich anwdchst, innerhalb 
welcher also die Interferenzreflexion stattfindet. Sie wird, ebenso wie ¢;, 
unendlich groB fiir (7}b;)=0 und ist um so kleiner, je gréBer die Feld- 
komponente senkrecht zur reflektierenden Netzebenenschar ist. 


Liegt die Zeitspanne von ¢; —|/m/2e, (6,)| bis t+ | m/2r69| (%b,)| 
ganz innerhalb der Beobachtungszeit, so wird der Reflexionskoeffistent 
fiir die durch 6; charakterisierte Netzebenenschar 


aM sont ()? m|U,)2 
Res (+22 by = i 
i Tan (O®  — Meg lO” i42) 


Ersetzt man hier noch den Fourter-Koeffizienten U; des Gitterpoten- 
tials durch die Breite AE;=2)|U,) der Energieliicke, welche mit der 
betrachteten Reflexion gekoppelt ist}, fiihrt man ferner den Winkel «; 
zwischen dem Feld und der Normale der betreffenden Netzebenenschar 
ein und ieee schlieBlich, daB 1/6; gleich dem zugehdérigen 


1 Veal. z.B. den Bericht von A. SoMMERFELD u. H. BETHE im Handbuch der 
Physik, Bd. XXIV/2, bes. S. 385 ff. 
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Netzebenenabstand d; ist, so kann man statt (12) auch schreiben: 


m (AE;)* d; 


Haba of aoe “t3) 


Dann wird aber wegen der vorausgesetzten Kleinheit der Reflexions- 
wahrscheinlichkeit die Durchtrittswahrscheinlichkeit durch die verbotene 
Zone D;=1—R; exp (—R,) gegeben durch 

{ m(AE;)?d; | 


D; ss exp } 40% ey F cosa; J” 


Diese Formel entspricht der bekannten ZENER-Formel fiir den eindimen- 
sionalen Grenzfall schwacher Gitterstérung, enthalt aber dariiber hinaus 
auch die von W. FRANZ! angegebene Abhangigkeit von der Feldrichtung 
relativ zum Kristall. 


(14) 


Koln, Institut fiir Theoretische Physik der Universitat. 


1 Vgl. den zusammenfassenden Bericht von W. FRANz in den Erg. exakt. 
Naturwiss. 27 (1953), insbes. S. 13—20. 
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Daimpfung von Ultraschallwellen in Fliissigkeiten. 
Von 
S. PARTHASARATHY, C. B. Trpnis und M. PANCHOLY. 
Mit 1 Figur im Text. 


(Eingegangen am 6. November 1954.) 


Zur Messung der Absorptionskonstanten von Ultraschall in Fliissigkeiten wurde 

eine verbesserte optische Methode angewandt. Das Verfahren ist schnell und ein- 

fach und liefert Resultate, deren Genauigkeit der anderer Prazisionsmethoden 
vergleichbar ist. Ergebnisse fiir 10 Fliissigkeiten werden angefihrt. . 


1. Einleitung. 


Bekanntlich werden Ultraschallwellen in Flissigkeiten starker ge- 
dampft, als man nach der Theorie von STOKEs [/) erwarten sollte. Um 
diese Anomalie zu erklaren, sind mehrere Hypothesen entwickelt worden. 
Wahrend indessen Ubereinstimmung dariiber besteht, daB die experi- 
mentellen Werte allgemein tiber den theoretischen liegen, weichen die 
Angaben tiber nach verschiedenen Methoden bestimmte Absorptions- 
konstanten dennoch stark voneinander ab. Im Hinblick auf diese Unter- 
schiede werden in diesem Laboratorium seit ein paar Jahren Messungen 
an Flissigkeiten nach verschiedenen Verfahren durchgefiihrt. Bisher 
wurden Resultate veréffentlicht, die nach mechanischen Methoden [2], 
[3], [4], nach kiirzlich entwickelten thermischen Methoden [5] bis [11) 
sowie schlieBlich nach einer verbesserten optischen Methode [12], die 
durch die gleichzeitige Beugung zweier Lichtbiindel gekennzeichnet ist, 
erhalten worden sind. Die nunmehr vorliegenden Resultate wurden 
nach einem Verfahren bestimmt, das eine Kombination der letzteren 
Methode mit einem optischen Verfahren von WILLARD [13] darstellt. 
Die Vorteile des Verfahrens liegen in Schnelligkeit und Sauberkeit der 
Ausfihrung. 


2. Experimentelles. 


Die experimentelle Anordnung stimmt im wesentlichen mit der in [2] 
beschriebenen iiberein und ist in Fig. 1 dargestellt. Zwei schmale mono- 
chromatische Lichtbiindel (A= 5461 A) mit parallelen Achsen durch- 
queren die Fliissigkeit in der Kiivette und kreuzen das Schallbiindel 
senkrecht zur Ausbreitungsrichtung. Die an den Schallwellen entstehen- 
den Beugungsbilder werden in der Brennebene zweier passend gewahlter 
Kameralinsen beobachtet. Als Schallquelle dient ein in X-Richtung 


——, 
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geschnittener Quarzkristall, der durch einen Hochfrequenz-Oszillator 
variabler Frequenz angeregt wird. Die beiden Elektroden am Kristall 
sind mit einem Rédhrenvoltmeter verbunden, bei dem das Anzeige- 
instrument durch ein Galvanometer ersetzt ist. Mit dicsem Galvano- 
meter kann man die an den Elektroden anliegende Wechselspannung 
ablesen. Die Spannung ist in der Endstufe des Generators willkiirlich 
regelbar. 


Die Oberflache der Fliissigkeit ist durch einen ausgehéhlten Schaum- 
glasblock abgeschlossen, um die Ausbildung stehender Schallwellen zu 
vermeiden. Die Messung geht in folgender Weise vor sich: 


Die Wechselspannung am Quarz wird, bei Null beginnend, stetig 
erhoht, bis die Lichtintensitat der Linie nullter Ordnung des unteren 


fenster £ 
Kollimator- Kuvette 


Kamera-Linse Auge 


Lampe ‘ wited | qh 
=—TH< | | 
pee AL so =F ( 
vad sa MN — ——_}  Kamera-Linse a fhefy- 
sate Fenster yrrctay — Flassigkeit meter 


Fig. 1. 


Beugungsbildes ein Minimum erreicht. Bei der Beobachtung durch 
ein Kathetometer werden die Linien héherer Ordnung durch einen 
Schirm abgedeckt, so daB nur die Linie nullter Ordnung im Gesichts- 
feld ist. Die diesem Minimum korrespondierende Spannung wird am 
Galvanometer abgelesen. Vorausgesetzt wird, daB die Schwingungs- 
amplitude des Quarzkristalles proportional der Amplitude der an- 
regenden Wechselspannung ist. Dann gilt fiir die Amplitude A, der 
Schallwelle am Beobachtungsort im Abstande x, von der Oberflache 
des Kristalles: 

hols ll ill 
oder 

a, = const Vi 2", (1) 


wenn J, die abgelesene Spannung ist, « die halbe Absorptionskonstante 
der Fliissigkeit und A, die Schallamplitude an der Oberflache des 
Kristalles. 

Daraufhin wird die Spannung derart gedndert, daB die Beugungs- 
linie nullter Ordnung des oberen Lichtbiindels (Abstand x, von der 
Oberflache des Kristalls) ihr Intensitatsminimum erreicht. Ist die ent- 
sprechende Spannung J so gilt 


A; = const V,*e~**. (2) 
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Da nach dem erlauterten Einstellverfahren die Amplituden der Schall- 
wellen A, und A, gleich sein miissen, folgt aus (1) und (2) 


_ MV,/V; 


(%_ — 44) 


av? = logig Ve — logiy V,/0,4343 «d+ 2”, (3) 


wobei d der Abstand der Achsen der beiden Lichtbiindel ist, und » die 
Frequenz des Schalles. Aus einer gréBeren Anzahl von Ablesungen 
wurden V, und V, bestimmt und «/y? als Mittelwert der jeweiligen Be- 
obachtungsreihen ermittelt. Die Frequenz des Hochfrequenzgenerators 
wurde mit Hilfe eines Prazisions- (Absorptions-) Wellenmessers fest- 
gestellt. Die Wechselspannung wurde nur fiir kurze Zeiten am Kristall 
aufrechterhalten, um eine unzuldssige Erwarmung der Fliissigkeit durch 
absorbierte Schallenergie zu vermeiden. 


3. Ergebnisse und Diskussion. 


Die Untersuchungsergebnisse sind in Tabelle 1 angegeben, die 
experimentellen Werte fiir «/y? in Spalte 5. Die Resultate sind in guter 
Ubereinstimmung mit solchen, die nach anderen Prazisionsmethoden er- 
halten und in Spalte 7 aufgefiihrt worden sind. Die Werte stimmen gut 
mit den Werten der optischen Methode | 72) iiberein, die der Vorlaufer des 
jetzt angewandten Verfahrens ist. Die Resultate bestaitigen anomale 
Absorption in den untersuchten Fliissigkeiten. Sie sind mit anderen in 
diesem Zusammenhange gewonnenen Daten im Hinblick auf eine még- 
liche Klarung der anomalen Absorption analysiert worden. 


Tabelle 1. 

ae en an SR as 2 ee Ce 
Nr. | Substanz on 5 a in — Friitherer Wert a 
1 ~~ | 3 4 5 6 | 7 
a a 
1 Bengols.... © 5 4,911 29,7 906,2 830,0 13 
Se OU Oc Sil x os 4,895 31,0 96,12 90,0 13 
3 | m-Xylol . .. 6,896 30,2 75,85 78,0 13 
4 Chloroform. 4,900 30,1 388,00 380,0 13 
5 | Nitrobenzol . 6,814 30,0 94,05 90,0 14 
6 | Hexylalkohol . 4,901 29,8 134,5 
7,| Butylazetat . 4,895 29,1 130,7 
8 Propylazetat . 4,885 29,7 144,7 
9 | Amylazetat. . 6,844 29,5 106,5 
10  Athylazetat a) 4,905 29,8 157.4 516 beiv=10% sec-1! 75 

| Athylazetat b) 2,961 29,8 252,3 37 beiv=69-108sec 76 


Athylazetat c) 6,825 20,8 110,4 
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Berichtigung und Erganzung 


66% 


zur linearen Feldtheorie“*. 
Von 
WILLY SCHERRER. 


(Eingegangen am 6. Dezember 1954.) 


Eine numerische Angabe iiber das Weltalter wird korrigiert. AnschlieBend werden 
die allgemeinen Folgerungen iiber das Weltalter gezogen, welche die vorliegende 
Theorie gestattet. 


I. Bei der numerischen Auswertung der periodischen Variante ist 
ein Versehen unterlaufen. Der in §6 (23c) angegebene Wert muB ersetzt 
werden durch 


%q = — 3,8° 10%’ cm. (1) 
Fiir das ,,Weltalter“ erhalten wir daher an Stelle von §6 (24) 
t= 7 4 *© — 0,67 - 10% sec = 2,1 - 10° Jahre. (2) 


Dieser Wert steht nicht im Einklang mit dem aus der Radioaktivitat 
erschlossenen Mindestalter von etwa 4+ 10° Jahren. 


Ganz allgemein schlieBt man aus §6 (133) auf 


2» 
Arctg (—, | 
t= mage Me 2,2+10° Jahre. 
3ke 2p 3h i hai (3) 
3R 


Wenn man also Anschlu8 an die radioaktiven Daten gewinnen will, 
mu man auch in der vorliegenden Theorie die periodische Variante 
ausschalten. Einen Vorbehalt zu dieser Folgerung werde ich unter IV. 
machen. 


Il. Als einzige Méglichkeit zur Verbesserung bleibt nur die erste 
Variante, also die aperiodische Variante mit singulaérem Anfang [§6 (6,) 
und (43,)]. 


Setzen wir 


K = —A, (4) 


* Siehe Z. Physik 139, 44 (1954). 


i 
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so kénnen wir das Weltalter fiir diese Variante durch folgendes Integral 


ausdriicken 
a 


ror Zu 
eae. a 
. | | 40K (5) 
le hr = 


co Oc? 


wobei « die Massendichte darstellt. 
Im Rahmen der NeEwronschen Gravitationstheorie findet man fiir 
die entsprechende Zeit das Integral 
“ 


i eee ree [ BP cue! (6) 
2}Gne 4 ll Vut+ Bui 


wobei C die NEwronsche Gravitationskonstante ist. 
Sollen die beiden Integrale am Anfang der Entwicklung, also fiir 
hohe Dichten uw, tibereinstimmen, so muB gelten 


oe (7) 


Gestiitzt auf (7) und (4) findet man weiter 


a 
av Sh [tie SES : (8) 
und die Auswertung von (5) liefert 
das i 14° hk. fat biel xu 
ie 3kce OW 1—W |’ tice |) 3? (9) 


Wenn man schlieBlich vermittels »—>7y zur periodischen Variante zu- 
riickkehrt, geht (9) uber in 


2 1 . | /xpe ie 
= “ty, W= 4 
t Pu Arctg W, | 3h 


Die beiden Varianten sind also voneinander durch die kritische Dichte 


(10) 


ut = 3” _4,65-10-*%gcm 3 (11) 
xc 
getrennt. 
Durch die Analyse von (10) und (9) gelangt man nun leicht zu 
folgendem 


Satz1. Variiert die Dichte monoton von -+- 0% bis 0, so variwert 


das Weltalter t monoton von 0 bis +0. . 
Nach den Angaben der Astronomen ist die Dichte 1m Intervall 


4078 <= py AOFM (12) 
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zu suchen. KIENLE?) schlieBt die obere Grenze aus und schatzt, ~ werde 
am echesten zwischen 10~2% und 10728 liegen. Unsere Formel (9) liefert 
folgende Einzelwerte: 


/ = 10° gcm~*; t; = 4,7-10° Jahre (13) 
fig 3°10" Gam, t, = 3,6-10% Jahre. (14) 


Wir kommen also auf Grund der empirischen Daten gerade auf diejenige 
GroéBenordnung, welche dem aus dem Atomzerfall erschlossenen Welt- 
alter entspricht. 

III. Die an Hand von (6) vorgenommene Eichung des aperiodischen 
Falles, welche auf (7) und (8) fiihrt, entspricht nicht genau der an Hand 
von §6 (14) vorgenommenen Eichung des periodischen Falles, die ihrer- 
seits auf §6 (18) und (24 a) gefiihrt hat. Ich habe hier die frithere Eichung 
fallen lassen, weil im aperiodischen Fall keine minimale Dichte jig 
existiert. Man stellt aber unschwer fest, da8 der Zusammenhang der 
beiden Eichungen geregelt wird durch 

Satz 2. Beim Ubergang von der aperiodischen Eichung zu der periodi- 
schen muB man in den Formeln (9) und (10) u durch u/3 und in (7) K 
durch K/3 ersetzen. 

Wie man aus (11), (13) und (14) ersieht bleiben wir damit immer noch 
im empirisch zulassigen Intervall. Welche Eichung den Vorzug verdient, 
kann vorderhand wohl kaum entschieden werden. 

IV. Das differentielle Zerfallsgesetz fiir eine groBe Zahl N radio- 
aktiver Teilchen lautet bekanntlich 


re (15) 


wo « die Zerfallskonstante bedeutet. Beziehen wir uns statt auf die 
Zeit ¢ auf die Lange x,=c?, so haben wir zu schreiben 

yoTS (16) 
Die unbeeinfluBbare Konstanz der GréBen « bildet die Grundlage aller 
radioaktiven Altersbestimmungen. Stellen wir nun neben (16) die Defi- 
nition der HupsBiE-Zahl in der Gestalt 


is 
r=*, (17) 


so kann man die Frage aufwerfen, inwiefern nun die rechte Seite von (16) 
weniger verdnderlich sein soll als die rechte Seite von (1 7). Ich kann 
mir nicht vorstellen, wie man dafiir empirische Griinde geltend machen 
will. Diese Griinde sind hypothetischer Natur und wurzeln letzten 


} KIENLE: Naturwiss. 38, 92 (1951). 


; 
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Endes in der Vorstellung, es handle sich eben um einen diskontinuier- 
lichen Elementarakt. 

Aber auch (17) bringt ein kosmisches Elementarphinomen zum Aus- 
druck, und sobald man einen Zusammenhang zwischen der singuliiren 
Anfangsphase und dem Atomzerfall herstellen will, mu8 man auch eine 
Proportionalitat von « und & in Betracht ziehen: 


Me GCE 
N =e (18) 


r ke OL 
Hierin wiirde nun & in Abweichung von (17) nur die momentan gemessene 
HvuBBLE-Zahl bedeuten. Die Integration von (48) ergibt 


N B= 
et Rae (19) 
Nach (49) waren dann alle Altersbestimmungen nur Aussagen relativ 
zum Weltalter. Sollte eine vollstandigere kosmologische Lésung wieder 
auf den periodischen Fall zuriickfiihren, so béte also der Ansatz (18) 
eine Méglichkeit, um einem Widerspruch mit der Endphase auszuweichen. 
Vorderhand kénnen wir also folgendes feststellen: Weil wir auf die 
aperiodische Lésung gefiihrt worden sind, erzeugt die Annahme der 
absoluten Konstanz von « keinen Widerspruch. 


Bern, Mathematisches Seminar der Universitat. 


Zeitschrift fiir Physik, Bd. 140, S. 164—180 (1955). 


Zur linearen Feldtheorie. II. 
(Schwache Felder.) 
Von 
WILLY SCHERRER. 


(Eingegangen am 25. November 1954.) 


Fiir schwache Felder gilt ein System von Differentialgleichungen, welches eine 

Verallgemeinerung der Maxwetischen Vakuumsgleichungen darstellt. Eine 

spezielle Klasse bilden diejenigen schwachen Stérungen, deren zugeordnete Gra- 

vitationspotentiale ungestért bleiben. Sie breiten sich mit Lichtgeschwindigkeit 
aus und haben genau die Struktur der elektromagnetischen Wellen. 


$1. Einletiung. 

Die vorliegende Abhandlung setzt die Kenntnis meiner ,,Grundlagen 
zu einer linearen Feldtheorte’} (zitiert unter dem Zeichen G), sowie der 
daran anschlieBenden Note ,,Zur linearen Feldtheorte I‘‘* (zitiert unter 
dem Zeichen L I) voraus. Ihr Zweck ist zu priifen, wie weit die MAXWELL- 
sche Vakuumselektrodynamik in der linearen Feldtheorie enthalten ist. 


Zur Erleichterung des Lesers will ich aber die in G und LI gewon- 
nenen Grundgr6éBen und Grundbeziehungen noch einmal kurz zusammen- 
stellen. 


Im Zentrum steht das Wirkungsprinzip 


df Mdx=0 (4) 
S 
mit 
M=My (2) 
und 
3 
M=A,+>)(A;M;+ A, M). (3) 


t=] 


Aus ihm folgen die Feldgleichungen 


Dabei gilt 


ran OM 
e YA, “ 


} SCHERRER, W.: Z. Physik 138, 16 (1954). 
* SCHERRER, W.: Z, Physik 139, 44 (1954). 
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-und — dank der Beschaffenheit von 9% — 


A,r ! OM ra] WwW 
t 
2 


é 


' 
Zur Vereinfachung der Schreibweise werde ich mich der normierten 
Matrizen 


1 O 0) 0) 
Tac ho DA 4 (@) (0) (8) 
6) Oo —1 @) 
) 0) Oo —4 
~ und 
(0) € 0) 0 
= = 0 O 0|} 
a;%\| = 0 : (9) 
O LS 397 0|| 


bedienen. In den allgemeinen Formeln L I, §1 (3) fiir die 6 Invarianten 
M,, M; werde ich also setzen 


Ayo = 1; (shail 
qo — 4: qi — q?2 — g33— _ | (Sa) 
a4,=a°*=0 fir A+p | 
sowie 
Ay) =— Uo = &; Ao3 = — Ago, gil | (9a) 
alle iibrigen a,,,=0. f 
Nun empfiehlt sich natiirlich die Abkiirzung 
Ay = Ago; hy == Ay; usw., (8b) 


wobei jetzt der Index / in a, oder a’ bei Summation doppelt zahlt. 
Die traditionelle Summationsvorschrift kann hierauf folgendermaBen 
gehandhabt werden: Uber einen Index 4 wird suminiert, wenn er mit 
geradem Gewicht auftritt. 
Einen entsprechenden Vorteil sichern wir uns fiir die Matrix (9), 
wenn wir statt @,,, schreiben 4,7, wobei nun 7 den nach (9) eindeutig 
bestimmten, zu A ,,komplementaren‘' Index darstellt, also 


0=1, 1=0, Zaegre ve 2 (9b) 
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SchlieBlich kiirzen wir noch ab gemab 
a, = 455. (9c) 


Hier ist also zu beachten, daB bei der Anwendung der Summations- 
vorschrift sowohl der angeschriebene Index A als auch der nicht ange- 
schriebene Index 4 einfach zahlt. 

Im Sinne dieser Verabredungen erhalten wir nun fiir die 6 Inva- 
rianten folgende Tafel: 


M, = a, a, a {28 f28., (10;) 
My, = et ad ae (103) 
M; = On fE PS: (10s) 
WM, =i, a, a? p29: f7., (10,) 
M, = 7G, @, a’ f° fPB. : (10.) 
= "a, Pt re. (103) 


SchlieBlich noch einige Bemerkungen zu den aus der Kosmologie 
erschlossenen Beziehungen zwischen den Kombinationszahlen A,, Ag, 
A,, A,, Ay, Ag. Dieselben lauten nach LI, §4 (7) 


2A, + A, + Ay + 27? (A, + 2A,) = 0, (113) 
2A, + A,+ 2A,=0, (119) 

e4,=0, (113) 

3y2—e +0; A, +2A,+0. (12) 


Die Bedingung (11,) erwies sich als unerlaBlich, um die Isotropie des 
Raumes zu gewahrleisten. Die Bedingungen (11,) und (41,) waren dar- 
tiber hinaus notwendig, um dem Raum Kugelsymmetrie zu verleihen. 
Die Bedingung (12) schlieBlich verhindert, da8 durch die Bedingungen 
(11) die Wirkungsfunktion M annulliert wird. 

Die-zur Verfiigung stehenden Konstanten A bilden also immer noch 
eine dreidimensionale lineare Schar. Doch wollen wir diese Freiheit 
nicht einschranken, bevor dafiir zwingende Griinde vorliegen. 


$2. Die Feldgleichungen. 


In der Einleitung haben wir alles zusammengestellt, was wir brau- 
chen, um die Feldgleichungen aufzustellen. Gestiitzt auf (2), (3) und (10) 
berechnen wir die Ausdriicke (5) und (6) und setzen dieselben in (4) 
ein. Es empfiehlt sich, alles entsprechend (3) in Komponenten aufzu- 
lésen. In diesem Sinne zerlegen wir also die Tensordichten t*#”” und 


| 
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yA Bon. die! K. Auy ZA 2 r : 
ae e Komponenten t;*“", t; und J", TA". Die Rechnung 
iefert: 


tH Day ag atyr yh” Phy, (11) 
the a, (YRH AY — ye fram (14) 
ther a, (yeh? — yr phn faery, (14) 
BAM" = 20,0, a? yet yh” Ey, (1,) 
ft = 20, dat yn yh Phy, (i) 
et E PO ey + 

EER, ey Py. 
Entsprechend ergibt sich fiir die Komponenten von 37” die Tabelle: 


Tne 7s M, pre = 4a, a3 a? ye" ph Posy 


(21) 
Th" = M, y>" — 2a, (y>" fae es foc po i) y, (25) 
Taek — May" — 2a, (ye PE PS — POV Ph) y, (2s) 
Jae es M, p*" sa 4a,a; a? yy # FB, ak (2;) 
x, Be M, y* | ee 44,4, a’ ye fP8: a y, (Za) 


Ei = My" — 24, y*" PE Ley + | 
+a, (PO4 PE — PR POM) yp. 
Fiir den Fall der bedingten Variation ware hier noch hinzuzufigen 
der dem %,—y entsprechende Term 


TyeX = Moy", (2o) 


wahrend offenbar das zugehdrige t)*”” Null ist. 
Gestiitzt auf die Tabellen (1) und (2) erhalten wir schlieBlich t*”” 
und &?“ als Kombinationen: 


3 
for es D (Ay te’ + 4,4"), (3) 
i=1 
.) I] ~ A Fi 
Dh = Ny tg PV te Ay"). (4) 
i=1 


Es ist aber wohl zu beachten, daB die hier und in (1) und (2) auftretenden 
Zeiger 1, 2,3, 0 Invarianten numerieren und demnach nicht mit Formen- 
oder Koordinatenzeigern verwechselt werden diirfen. 

Die Ausdriicke (3) und (4) waren nun also in §1 (4) einzusetzen. Zur 
expliziten Berechnung der Feldgleichungen miiBten schlieBlich noch die 
Ableitungen der t””” ermittelt werden. Damit kamen die zweiten Ab- 
leitungen, also auch die Kriimmungskomponenten ins Spiel. Fiir die 
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Zwecke der vorliegenden Arbeit bendtigen wir aber die Kriimmungs- 
komponenten noch nicht. Wir kénnen also ihre explizite Berechnung 
aufsparen bis zur Analyse der Gravitationsfelder. 


§3. Schwache Felder. 

Eine invariante Charakterisierung schwacher Felder erhalten wir am 
einfachsten durch folgende Forderung: Die FeldgréBen /;,; seien so 
klein, daB man ihre Quadrate vernachlassigen darf. 

In nullter Naherung gilt daher 


fi.z1 = 0 (1) 
und es folgt 3 
__ OV 
ea Sag (2) 


wobei y;, ein variabler Formenvektor ist, dessen Komponenten wir als 
unabhangige Funktionen voraussetzen miissen, weil fiir die zugehdérige 
Jakobiante gilt 

CY. 
Cx, 


y = Det |p; ,| = Det +0. (3) 
Jetzt fiihren wir ein ausgezeichnetes Koordinatensystem ein durch 
den Ansatz 


Xx; = yj. (4) 
In diesem System gilt nach G §2, (1) und (11) 
Vik = San Yia 
Oxy I: 
CX Via 


Gestiitzt auf (2) und (4) folgt also 


Min 8%, Om 
Wir sind somit zu folgendem Ergebnis gelangt : Zu einem schwachen Feld 
existiert immer ein ausgezeichnetes Koordinatensystem derart, daB in 
demselben die Basismatrix in nullter Naherung diagonal und folglich 
die zugehérige Metrik pseudceuklidisch ist. 
In diesem Koordinatensystem wollen wir jetzt arbeiten und der Be- 
quemlichkeit halber die Querstriche fallen lassen. 
In erster Naherung kénnen wir nun ansetzen 


Cx; OR; ay - 
aut = bine (5) 


d 


Vir = Ov» 4 Pik: (6) 
Dabei sind die @;,, so klein, daB man ihre Quadrate vernachlissigen darf. 
Also gilt 


y= 6.4 — Dir (7) 
und 


y=i1+%.. (8) 
wobei in (8) itber 7 summiert wird. 


~ ew 
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Also folgt 
; 203 OP. OQ; )\ 
| td oa eine oad 9) 


Gestiitzt auf §2 ist es nun leicht, die Feldgleichungen in erster Nahe- 
rung anzugeben. Dabei wollen wir, entsprechend der quasieuklidischen 
Metrik, nur den Fall A,=0 in Betracht ziehen. 

Ohne Schwierigkeit erhalten wir die Tabelle 


{, 44” = 2a* a, a, f¥>, (10,) 
ther = a, fh? — a, fra, (104) 
ther = a, OME — a MPE, (103 
i>" = 2a, @, f, (10;) 
t,he" = 24'4, a,” fed, (103) 
her = a, OR PO + 4a, (Oh Ae — ohm HE. TO, 


Die Energiekomponenten §2 (2) aber werden offenbar von zweiter Ord- 
nung klein. Wir haben also einfach zu setzen 


TAF = 0. (14) 


SchlieBlich ist es, gestiitzt auf (6) und (7) sowie auf §14 (8a), in erster 
Naherung erlaubt, alle oberen Indizes auf der rechten Seite der Gln. (10) 
herunterzuziehen. Fiir die korrekte Durchfiihrung dieser Indexverschie- 
bung hat man die im AnschluB an die Gln. (13) angegebene Vorschrift 
zu beachten. Wir erhalten so fiir die Teilkomponenten der Feldglei- 


chungen 
a t” uy 


Seale 0 (12) 
Ox, 
fiir schwache Felder die Ausdriicke 
a,AuyY at. 
oh = : Ofi, py 4 
Se = 24,4, 4, ax, ’ (131) 
dt,” Cf, ps OfuvA 4 
a : 
Ox, ~ ax, “aaah oe ead (132) 
atn* Che, uo Ofo,ve { 
— a 
Ox, #\ Ox 02, 4 ” OF (133) 
gia ryt (G,) 
a a 
Ox, lee Nae hae Cte eo 1 
Ap c of x : 
i ar 24,4, a, - sa , (4 2) 
Ox, " 0%, 
oa, A4,¢Y B 4 of cas 
a z Stoke 4 1g gy (Blunt _ 5 728 | 1 
ax, bag, oe ge A Ox, un Ox (133) 
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Wie schon gesagt beziehen sich die Gln. (6) bis (13) auf ein ausge- 
zeichnetes Koordinatensystem. Von dieser Beschrankung kann man 
sich aber frei machen fiir alle diejenigen affinen Koordinatentransforma- 
tionen, welche zusammen mit einer passend gewahlten Formentransfor- 
mation den Ansatz (4) invariant lassen. Der Ansatz (4) zeigt schon, 
daB man in diesem Falle jeden Formenzeiger gleich behandeln muB wie 
einen gegenstandigen Koordinatenzeiger. 

Eine elementare algebraische Analyse zeigt, daB man die genannten 
Erfordernisse am bequemsten erfiillt durch folgende 


Vorschrift: Man behandle jeden Formenzeiger wie einen gegenstdéndigen 
Koordinatenzeiger und beschranke die linearen Koordinatentransformatio- 
nen auf LORENTZ-Transformationen. 


Wir setzen daher 
Pig = 4,0; ,, (14) 


wo nun w;, einen kovarianten Koordinatentensor zweiter Stufe darstellt 
und erhalten aus (9) 


fint = 4: Nj, x1 (15) 
mit 
ps 23 CW; & Cw; 1 
hi = 2 Ge ae) (16) 


Die Relationen (13) gehen damit iiber in 


a in 
Oty vid oe SAE Ohy ue 
Ox, gle ae ee (17;) 
Aue 4 P 
Oty si J — @a Chy, wa é hy, wa 
e By uy éx, @ %, } ’ (1 7») 
Aur 2 a 
ot,” ae ch 
- —— aa (- 2. He — e@,"e@ 
ax, oe ie 93, 1 % ax, ), (17s) 
TA my A 
ot," — - Oh, ag 
ar th nee — 
Ox, edict TR (17) 
az Amy s 
Cty —9544 Oh, 08 — 
Ox, a ar (172) 
ye’: Cy ila oe i. ah 
= = a - , a a rt Shs Be. yaa we = 
Ox, 7M Ox, a EN Bes 93, 1% ox, | (175) 


Die endgiiltigen Feldgleichungen fiir schwache Felder lauten nun also 


; 


3 
ayAur me 
at  ’ v athe? ete” 
=)'(4 +4, 26") =o. (18) 
t=} 


Ay 
Ox, 


Ox, ral Xy 
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_ Zusammen mit den aus (16) folgenden Gleichungen 
/ F Oh) xo 
ater wa ag (19) 


_ liegt also eine Vervierfachung der MAXweELtschen Gleichungen fiir das 
~ Vakuum vor. 


$4. Lichtwellen. 


In G §7 haben wir die der Basismatrix ||; ,|| zugeordnete Metrik 
erklart und darauf hingewiesen, daB diese Zuordnung nicht von vorne- 
herein eindeutig feststeht. Wir wollen jetzt eine Entscheidung treffen, 
indem wir setzen 


Gay =O Wi ar, n- (1) 


Dieser Festsetzung geniigt als Spezialfall das in LI §4 fiir die kosmolo- 
gische Lésung gebrauchte Linienelement. 

Gestiitzt auf die Relationen (6) und (14) des vorausgehenden Para- 
graphen erhalten wir daher in erster Naherung fiir schwache Felder 


G, = 7 + DW) + Diy i- (2) 


Es scheint schwer, ja vielleicht unméglich, aus dem metrischen Feld 
in invarianter Weise Lichtwellen zu gewinnen. Wir richten daher unser 
Augenmerk auf diejenigen schwachen Felder, deren Gravitationspoten- 
tiale nicht nur in nullter, sondern auch in erster Ordnung konstant sind: 


O,,+@,,=0. (3) 


Die weitere Entwicklung verlauft nun ein Stiick weit nach klassi- 
schem Vorbild. Durch die Gln. §3 (19) sind die GréBen nur bis auf ein 
Gradientenfeld 67,/0 x, festgelegt. Uber letzteres kénnen wir daher noch 
so verfiigen, daB die w,,, fiir jeden Index 4 der LorENtz-Konvention 
geniigen. Dieselbe lautet in unserem Falle in erster Naherung 

ew," 


= gt—+# —0, 


Ox, Oxy 


und dafiir diirfen wir schreiben 


peed 0. (4) 
x 


SchlieBlich benétigen wir noch den D’ALEMBERTschen Operator LJ. In 
unserer Schreibweise lautet er 


eto. (5) 


” Ox, Oty 
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Fiihren wir nun die Ausdriicke §3 (16) in die Beziehungen §3 (17) 
ein und machen wir Gebrauch von den eben gewonnenen Gln. (3) und (4), 


so erhalten wir 
A,pyv 


Ot 1 
= as w ; 
s Oo, (6) 
54 ABP 
a ye oa (62) 
Oxy 
Ot,”"? 
x" =0, (6,) 
Ox, 
A,pey 4 OW, 5 f 
ot 2h @ (a BA) 
ney oye ae : 
i eg oe) (6,) 
athe” x, 2 (a, one ; 
== == Jeg Sy 
ax, 2 “ 6x, v ex, ( 2) 
7 Apr — tal Cw 
ex, 2" day \* oy (63) 


Aus diesen Ausdriicken erhalt man durch die vorgeschriebene lineare 
Kombination offenbar dann und nur dann Lichtwellen, wenn die Be- 
dingungen 


— CW) 5 
a, “Oa, =0 (7) 
erfiillt sind. Da neben diese Bedingungen noch die LORENTZ-Konventio- 
nen (4) treten, sind die Méglichkeiten starker eingeschrankt als in der 
Maxwe.tschen Theorie. 
Zusammenfassend stellen wir also fest: Um Lichtwellen zu erhalten, 
miissen wir die Gleichungen 


Oe, =0 (8) 
unter den Nebenbedingungen (3), (4) und (7) lésen, wobei natiirlich die 
Kombinationszahlen A; der Bedingung 
genugen miussen. 

Zur Aufstellung einer konkreten Lésung empfiehlt sich der Ansatz 
Wo,1 = — 9 =A, p(x 


( 
Wo,3 = — Wy 9 = Ag P(x 
(10) 
( 


\ 
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Dabei sei y(y) eine Schwingungsfunktion 
P(x) = — v(x) (11) 
und ¥ eine Linearform der Koordinaten 
X=Pr%, (12) 


wahrend die A; und B; Konstanten sind. Durch (10) ist (3) erfiillt. 
Infolge (12) wird (8) durch 


a, p, p, = 0. (13) 


Geniige geleistet, so daB nur noch die Diskussion der Bedingungen (4) 
und (7) ibrig bleibt. Durch (10) erhalten sie folgende Gestalt: 


Ayp, + Agpe+ Asp3 = 0 


Ay po + Bsp2.— Byp3 = 0 (4’) 
Aspo— Bshy + B,p,;=0 
Aspy+ Bopi— Bhp = 0, 
eA; po +17 Asp, — A, $3 =0 
eA, p, — 7 By p, — 4 Bg p3 = 0 ' (7’) 
— €Bs py + €Agp, + HB, py = . 
€B, py + €Az hy +7 By ps =0. 


Eine elementare Analyse zeigt, daB man die einfachste Lésung erhalt, 
wenn man die Bedingung ~;-+$;=0 zugrunde legt. Mit Riicksicht auf 
(13) heiBt dies, daB man auszugehen hat von den Bedingungen 


p2=h3 = 0; pi=—Cbo; C= 1. (14) 
Es folgt 
Azz 0; Beme o-t er) A rel Be. (15) 


Dank der Antisymmetrie (3) kann man nun die vollstandige Losung, 
d.h. alle w,; , entnehmen aus der Tafel 


o1 = 0; O23 = Byp 
Wo,2= —0 Bs; 3,1 = Bay (16) 
Wy3= 6 Beg; as nip 
wobei 
P= {ho(% —f Ma)5 (17) 


die durch (11) charakterisierte Schwingungsfunktion ist. 
Gestiitzt auf §3 (16) und (15) erhalt man schlieBlich folgende Tafel 
fiir die FeldgréBen f;,;, wobei die Abkiirzung 


P=}ho9' (18) 
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beniitzt wird: 


fo,o1 = 9; fo23 = 0 

fooz= CBP; fos1 = — B,P } (199) 
fo.o3 = =—(igk; fore ia Bb, P 

hor = 93 hes = 0 | 

hoz= BP; fiai=—¢B,P (194) 
hios= — Bak. haz _tp,p,| 

feo1 = 9; fe23 = 0 

feo2=9; feai= CBP (19g) 
feos = Bye; foi2=9, 

fs,o1 = 95 hirer =9 | 

fs,o2 =— B,P; fa.31 = 0 (195) 
fs.o3 = 93 fhsi2= -B, P.| 


Es ist lehrreich, auf Grund dieser Tafel die in G §13 (8) und (9) de- 
finierten phanomenologischen Feldstarke /,, zu berechnen. Nach den 


ad 


eben zitierten Gleichungen kénnen wir sie schreiben in der Gestalt 


F, =RC,yfi yr, 


ur 
wobei wir entsprechend den getroffenen Verabredungen alie Indizes tief 


gestellt haben. 
Wenn wir noch durch 


B= F,5; Eo = Fig, E,=F,, | (21) 
H,=Fy,, H,=K, H,=F, J 
die klassischen Vektoren einfiithren, erhalten wir 
By =U: H,=0 
E,=kP([(B,(C,o4 C,) — B,C,]; HH, =C&, (22) 


E, = — kP[B,(Coo + C,) — B, Cs); H,= —CE,. 


Der phanomenologische Strahl F,,, ist also bei jeder Wahl der Kon- 
stanten C,, transversal und erfiillt daneben automatisch die Bedingungen 


C= €EH=0; [(€$) =— 0,0). (23) 


Er stellt somit genau eine elektromagnetische Welle dar. 

Die Berechnung der iibrigen simultanen Lésungen zu (13), (4’) und 
(7') ergibt zu jeder von (¢,0,0) verschiedenen Raumrichtung genau 
einen Strahl, Diese weiteren Lésungen miissen sich also auch ergeben, 
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wenn man auf die geschilderte Lésung passende Lorentz-Transforma- 
_tionen ausiibt. 

SchlieBlich sei nocheinmal darauf hingewiesen, daf diese Folgerungen 
_ fiir jede Wahl der Kombinationszahlen A;, A; gelten, sofern nur die Be- 
_ dingung (9) erfiillt ist: 24, — Ag+ 0. Sie gelten also speziell auch, wenn 
~ sich auf den symmetrischen Fall beschriinkt, indem man ¢ = “= 0 setzt. 


$5. Die Energiedichte. 

Bekanntlich ist die Energiedichte 7°.° beim heutigen Stand der 
Dinge die wichtigste Komponente des Energietensors. Wir wollen sie 
_ daher berechnen und auf den Fall der Lichtwellen spezialisieren. 

Man kann immer, d-h. also auch fiir starke Felder, das Koordinaten- 
system in einem Weltpunkt so wahlen, daB in demselben gilt 


Dr = 63,1 . (1) 
_ Dadurch vereinfachen sich die Komponenten des in §2 (2) angegebenen 
Energietensors. Um die Vorzeichen sichtbar zu machen, sondern wir 
den Index 0 aus und verwenden fiir die iibrigen Indizes lateinische Buch- 


_ staben. Es empfiehlt sich, zuerst die Teilkomponenten M,, M; der Wir- 
kungsfunktion 


ee 


M = 2 (A; M; + 4;M)) (2) 
zu notieren. an 
M, = — 2fo,01 fo,o1 + fo,xi fort + 2hiorhior — farhier (31) 
M, = — foorto.or + 2foertior + hiorho: + heihyie- (32) 
M, = — fo,orfo,0: o< 2foorfe, et “+ hioiton =fezilear (33) 
M, = 24, (ef,,01 + 1f,,23)" — | 3,) 
1 
— 4a,€n (ft, orby,23 +hozly31 + toats,12) + j 

M, = 44, (e fyo1 + Nf, 23)" (39) 
My = 2 (€ fy,01 + 11,23) fe,00° (33) 


In Verbindung mit §2 (2) erhalten wir nun weiter die Teilkomponen- 
ten der Energiedichte T°.°. 


T° = 2fo,01f0,01 + fo,nifo,n1 — 2fi,o1fi,01 — Kinrti,ny (41) 
T°.” = fo or fo,o1 — fi,orht,08 + hinihsin> (4) 
T° = fo o1fo,01 — fi,oih,o1 — fini lev (45) 
T,°.° = 24, (€ fy, 01 + Nhy,28)” + 44 (e fr,o1 — 1 fy,23) > (41) 
Tie! cm AG, (€f/,01+ 1 ty,23)” — 84, € fy 01 (Efs,01+  ty,23)» (42) 
T39:° = 4 (€ fo,o1 + 1fo,28) foo + 2 (€fi,o1 + 14,28) foie: (4s) 
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Bei der Spezialisierung auf Lichtwellen findet nun eine erhebliche Ver- 
einfachung statt dank der Lorentz-Konvention §4 (4) und der Strah- 
lungsbedingung §4 (7). 

Aus der Lorentz-Konvention folgt namlich wegen der Antisymme- 
trie von @, , 


OW4 CWy 2 ja] 
BOF, ie ae Cx; 
= 24, h,, Ap 
— 2h eke a 0, 
also 
Pa =e 0. (5) 
Die Strahlungsbedingung aber kénnen wir wegen 4, = — 4; schreiben in 
der Gestalt 
_ 6035 1 /— 25 _ 9035 
a, ad (a, Z oar as <a -} 
Oxy 2 Cx, Ox, 
as, rae wba 
= 2 é Cx, Ox; 
— a, hy sy 
= a, Qa ss =0 
Somit folgt 
a, hh, a 0, 


was ausfithrlich geschrieben liefert 
Efior +N fhi23 = 0. (6) 
Aus (6) ersieht man unmittelbar, daB fiir Lichtwellen gilt 
T,°.° -— T,°.° = T,°:° =, 


und aus (5) tolgt nach einer kleinen Rechnung 7,°.° = 0. 
Wir haben also nur noch die Kombination 


Ta C... « ] 0.0 ) | 0.0 
: 1%.9 =: A, T,° + Ay ey 
zu berechnen. 


Gestiitzt auf die Tafel §4 (19) erhalt man als Resultat 
T°,° = — 4(2A, —A,) B? P? 
oder also, unter Beachtung von §4 (18) ausfiihrlich geschrieben 
T°.° = — (2A, — Ay) Bi pig’®. (7) 


Die Energie erscheint also als definite GréBe. Doch sind wir hier noch 
an unser ausgezeichnetes Koordinatensystem gebunden. Es stellt sich 
somit die wichtige Frage, ob das Vorzeichen der Energiedichte der 
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Strahlung invariant ist gegeniiber den fiir die schwachen Felder maB- 
gebenden Transformationen, d.h. also gegeniiber den im AnschluB an 
§2 (13) erlauterten LorENtTz-Transformationen. 

In unserer Bezeichnung ist eine lineare Transformation 


Sore 


4 gde &e (8) 
eine LORENTZ-Transformation, wenn die Spaltenrelationen 


a, al gi? — a (9) 


oo 
erfuillt sind. Diese Relationen miissen wir also beachten, wenn wir den 
Energietensor gemaB 
“ye 
t Pec eet al Bg (10) 


: * : = ; 
auf ein neues SystemS transformieren. Speziell lautet nun unsere 
Frage, ob dem definiten 7°.° ein gleichsinnig definites 


70,0 = g 99° TH? (11) 
entspricht. 
Zur Beantwortung der Frage muB also der ganze Energietensor be- 
rechnet werden. Gestiitzt auf (5) und (6) findet man leicht 


Taf = Tot — Tt 7t* = 0, a2) 


so daB nur 7,*“ und T,*“ explizite zu berechnen sind. 
Die Ergebnisse kénnen wir in folgender Weise zusammenfassen : 


Tt = 2 1,7t, (13) 

tf 9 o 

aN 5 ae (aaa eae 8) 
| 7,4#|| — 02 B?|| : 14 
bre east Fl 4.) Py pe! ol i 

aad cate 

mit 

Q=fPoy’. (15) 


Die Auswertung von (11) ergibt nun tatsachlich ein zu T°.° gleichsinnig 
kk 
definites 7°: °. 
Es sei noch darauf hingewiesen, daf man fiir diese Auswertung 


speziell die Gleichung 
Anh’ af * =A (16) 


braucht. Sie ergibt sich wie folgt: Man schlieBt zuerst aus (9) auf die 


Inversionsformel | 
aM = Ay A, hy 4 (17) 
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und gewinnt hierauf aus (9) und (17) die zu (9) gehorigen Zeilenrelationen 


Ayal ao = 6 Bs ye: (18) 


Setzt man hier 9=o0=0, so folgt (16). 

Dieses Verfahren ziehe ich dem imaginaren x,=7t vor, weil es auto- 
matisch die Vorzeichen liefert. 

Es ist klar, daB diese Energiebetrachtung nur dem heuristischen 
Zwecke dient, den Anschlu8 an die MAxwettsche Theorie sicherzu- 
stellen. Nach dem leitenden Gesichtspunkt unserer Approximation 
sollten ja alle quadratischen Terme vernachlassigt werden, insbesondere 
also auch die Komponenten des Energietensors. 


$6. Schlupbemerkungen. 


Zum Schlu8 will ich die zwischen den Kombinationszahlen A;, A; 
bestehenden Relationen [§1 (11), (12); §4 (9) und §5 (7)] nocheinmal 
zusammenstellen. Dabei ordne ich sie nach MaBgabe des Gewichtes, das 
ich ihnen jetzt beimesse und iiberschreibe jede Relation mit dem fiir 
ihre Aufstellung maBgebenden Grund. 


I. Positiv definite Lichtenergie: 


A, — 2A,> 0. (I) 
II. Isotropie des Raumes: 
2A, + A, + Ag + 27? (A, + 2A,) = 0. (II) 
III. Kugelsymmetrie des Raumes: 
2A,+A,+24,=0; ¢eA,=0. (Illa, b) 
(37? — s*) (A, + 2A,) +0. (IIIc) 


Ich halte diese Bedingungen nicht fiir unumstéBlich. Aber sie entspre- 
chen dem, was sich bis heute bewahrt hat so gut, da8 man keine von 
ihnen preisgeben kann ohne zwingenden Grund. Am starksten hypothe- 
tisch belastet scheint mir III. zu sein, die Kugelsymmetrie des Raumes. 

Die wichtigste Frage, die noch offen steht, ist die, ob die antisymme- 
trische Matrix unvermeidlich sei, ob also zum mindesten 7 + 0 zugelassen 
werden miisse. 

Im symmetrischen Falle gilt e—%=0 oder, was praktisch auf das- 
selbe hinauskommt, 4, = 4y= A; =0. Dann aber gehen (IT) und (IIIa, b) 


liber in 


2A, a A, + Ay = i. (II’) 
2; + Ay + 243=0, (III’) 
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_ wahrend nach LI, §4 (5) und (6) die Bedingung (IIIc) durch 
2A; + A, + 343+ 0 (L11'c) 


_ ersetzt werden sollte, was nach (II’) und (III’) nicht mehr méglich ist. 


Sicher ist es, daB man (III’) preisgeben kann, wenn man sich nicht 
scheut, einen Torusraum einzufiithren, d.h. also einen Kubus, in dem je 
zwei gegentiberliegende Wande identifiziert werden miissen. In einem 
derartigen Raum existieren namlich dieselben kosmologischen Lésungen, 
welche ich in LI angegeben habe. 


Die Vereinfachungen im symmetrischen Falle sind so erheblich, daB 
ich es fiir zweckmaBig erachte, denselben vorderhand als Arbeitshypo- 
these einzufiihren. 


Insbesondere werde ich in der nachsten Mitteilung zeigen, daB im 
symmetrischen Falle neben den hier entwickelten Gleichungen fiir 
schwache Felder die E1NsterNschen Gravitationsgleichungen fiir be- 
liebige Felder streng giiltig sind. 


Nachtrag zur ,,linearen Feldtheorie II“. 


Im symmetrischen Falle lassen sich die Differentialgleichungen fiir 
schwache Felder, also die Differentialgleichungen §3 (17,) bis (17) und 
(18) mit A,—A,—=A,—=0 ohne Schwierigkeit vollstandig lésen und 
liefern folgende Ergebnisse: 


1. Die in §4 beschriebenen elektromagnetischen Wellen gemaB 


04, + @,,= 0. (4) 


2. Reine Gravitationswellen gemaB 


1 p — Oya = 0. (2) 
Sie breiten sich ebenfalls mit Lichtgeschwindigkeit aus. Die ihnen zu- 
geordneten FeldgréBen f;,,, sind samtlich Null. Also ist auch der zu- 
gehérige Energietensor fiir jede Wahl der Kombinationszahlen Aj, vas Ay 
gleich Null. Diesen Gravitationswellen kann daher zum mindesten bei 
Beschrankung auf schwache Felder keine physikalische Realitat zuge- 


schrieben werden. 


3. Zwei Typen von Materiewellen gemab 


O34, — Oya + 0. (3) 
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Der eine Typus kann durch eine geeignete LoRENTZ-Transformation 
in eine stehende Welle verwandelt werden und stellt demnach eine 
DE BroGLiE-Welle dar. Der andere Typus geht aus einer DE BROGLIE- 
Welle dadurch hervor, daB man in der letzteren die Zeit x) mit einer 
Raumkoordinate vertauscht. 

Die diesen Materiewellen zugeordneten FeldgréBen sind nicht alle 
gleich Null. Trotzdem verschwinden alle Komponenten des zugehérigen 
Energietensors infolge der Isotropiebedingung 


2A,+A,+A,;=0. (4) 


Hier haben wir also ein erstes Symptom dafiir, daB man vielleicht auch 
die Isotopie des Raumes preisgeben muB. 


Bern, Mathematisches Seminar der Universitat. 


’ 
. 
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Zur Theorie der Leitfihigkeit und Viskositiit 
von Salzschmelzen. 


Von 
W. OLDEKOP. 
Mit 4+ Figuren im Text. 


(Eingegangen am 7. November 1954.) 


In Anlehnung an die molekular-statistische Theorie der Fliissigkeiten von GREEN 

wird ein einfaches Modell der Salzschmelzen, insbesondere der Alkalihalogenid- 

schmelzen entwickelt, welches gréBenordnungsmaBig die experimentellen Befunde 
richtig beschreibt. 


1. Einlettung und Problemstellung. 


Wahrend die an der Elektrizitatsleitung beteiligten physikalischen 
Vorgange bei den kristallinen Salzen und bei den waBrigen Salzlésungen 
weitgehend untersucht sind, besteht hinsichtlich der beim Stromtrans- 
port und der inneren Reibung sich abspielenden molekularen Vorgange 
in Salzschmelzen wohl noch keine endgiiltige Klarheit. Die bekannten 
Vorstellungen tiber Platzwechselvorgange im Rahmen der_,,Locher- 
theorie‘‘! haben zwar bei vielen Fliissigkeiten ein qualitatives Verstand- 
nis des Zahigkeitsverhaltens ermdéglicht, doch laBt die relativ schwache 
Temperaturabhangigkeit der Leitfahigkeit von Salzschmelzen es fraglich 
erscheinen, ob die zu einer exponentiellen Temperaturabhangigkeit 
fiihrenden Platzwechselvorstellungen einen geeigneten Ausgangspunkt 
fiir die theoretische Behandlung dieser Substanzen darstellen. 

In der vorliegenden Arbeit soll eine von molekularen Verteilungs- 
funktionen ausgehende statistische Beschreibung der Vorgange in den 
Salzschmelzen versucht werden, die sich in ihren Grundgedanken an 
die in den letzten Jahren insbesondere von BorN und GREEN? ent- 
wickelte allgemeine statistische Theorie der Fliissigkeiten anlehnt und 
in gewissem Sinn eine zur Platzwechseltheorie entgegengesetzte Aus- 
gangsnaherung darstellt. Wenn es sich hierbei auch nur um einen 
ersten Versuch handelt und die makroskopische Berechnungsweise 
der atomaren Vorgange stellenweise recht gewagt sein mag, so soll 
die qualitativ und gréBenordnungsmabig richtige Wiedergabe der ex- 
perimentellen Ergebnisse doch nahelegen, daB die im folgenden ent- 
wickelte Auffassung einen méglichen Ausgangspunkt fiir die Behand- 
lung der Salzschmelzen darstellt. Die vorliegende Arbeit beschrankt 


x Vel. etwa EuckEN, A.: Lehrbuch der chemischen Physik. 
2 Vgl. GREEN, H. S.: The molecular Theory of fluids. Amsterdam 1952. 
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sich auf die Alkalihalogenide, doch diirfte eine Ubertragung der Grund- 
gedanken auch auf kompliziertere Falle méglich sein. 

Wir fassen das geschmolzene Alkalihalogenid als eine Ansammlung 
einwertig positiv und negativ geladener starrer Ionen bestimmten Durch- 
messers auf, die durch Coutoms-Krafte miteinander in Wechselwirkung 
stehen. AuBer den Coutoms-Kraften wirken auf die Einzelionen noch 
elastische StoBkrafte bei ZusammenstéBen mit Nachbarionen. Derartige 
StoBkrafte spielen in Gasen, wo der mittlere Teilchenabstand groB gegen 
die Reichweite der intermolekularen Kriifte ist, eine wichtige Rolle. In 
Fliissigkeiten dagegen, wo sich jedes Teilchen dauernd im Kraftfeld 
seiner Nachbarn befindet, sind es vor allem die standig wirkenden inter- 
molekularen Krafte (in unserem Falle die CouLomB-Krafte), die das 
makroskopische Verhalten der Substanz bestimmen. Wir wollen im 
folgenden die StoBkrafte ganz vernachlassigen und nur die COULOMB- 
Krafte zwischen den Einzelionen beriicksichtigen. Die Ermittlung der 
Leitfahigkeit und Zahigkeit wird dann auf Grund folgender Uber- 
legungen durchgefiihrt. 

Jedes Ion verdrangt aus seiner Umgebung die gleichgeladenen lonen 
und zieht die entgegengesetzt geladenen Ionen an, so daB es im statisti- 
schen Mittel von einer Raumladungswolke umgeben ist, die (genau 
so wie bei der DEByYE-Htcketschen Theorie der waBrigen Elektrolyte) 
aus Poisson-Gleichung und BoLTZzMANN-Verteilung berechnet wird. 
Bei kugelsymmetrischer Raumladungswolke wirkt auf das betrachtete 
Ion keine resultierende Kraft. Wird dagegen das Ion durch ein auBeres 
Feld bewegt, so bleibt die Raumladungswolke wegen der zu ihrem Auf- 
bzw. Abbau erforderlichen Relaxationszeit hinter dem Ion zuriick und 
iibt eine zur Geschwindigkeit proportionale Bremskraft aus, die bei 
gegebenem Feld die Wanderungsgeschwindigkeit des Ions bestimmt. Aus 
den Wanderungsgeschwindigkeiten erhalt man die 
Leitfahigkeit. Der wesentliche Unterschied gegen- 
iiber ahnlichen Betrachtungen bei der DEBYE- 
Htcxketschen Theorie der waBrigen Elektrolyte 
’ wird in der hier anders vorzunehmenden Ermitt- 
lung der Relaxationszeiten liegen. 

Zur Bestimmung der Zahigkeit wird ein makro- 
ne 1. Deformation der  skopisches Geschwindigkeitsfeld v, = const - z (vgl. 
ugelsymmetrischen —La- - . 
dungsverteiling durch ein. Fig.) betrachtet. Die Raumladungswolke eines 

arene eee in der Ebene z=0 befindlichen Ions wird dann 

y in der skizzierten Weise deformiert werden. Die 

im oberen Halbraum (z> 0) befindlichen Ionen iiben wegen der asym- 
metrischen Ladungsverteilung um das Zentralion auf dieses eine Kraft 
in der x-Richtung aus. Summiert man die Kraftwirkungen, die von 
den oberhalb der z-Ebene befindlichen Ionen auf die unterhalb der 
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z-Ebene befindlichen Ionen ausgeiibt werden, so erhalt man die zum 
Geschwindigkeitsgefialle gehérige Reibungsspannung. Aus der Reibungs- 
spannung ergibt sich die Ziahigkeit. 

Fir die quantitative Durchfithrung dieser Uberlegungen miissen die 
GréBe der Raumladungswolke und ihre Deformation bei den entspre- 
chenden Geschwindigkeitsfeldern sowie die Relaxationszeiten ermittelt 
werden, was in den folgenden Abschnitten getan wird. Im letzten Ab- 
schnitt werden dann die Ergebnisse der Theorie mit den experimentellen 
Befunden verglichen. 


2. Die das Einzelion umgebende Raumladungswolke 
und die daherriihrenden Kraftwirkungen., 


Bezeichnen wir mit ;" die Teilchendichte der Ionenart j in der Um- 
gebung eines Zentralions der Sorte 7, so ist 


—qo 
ss « RT 
nl =H -¢ ; (4) 


wobei 7 die (bei Alkalihalogeniden fiir beide Ionenarten gleiche) mittlere 
Teilchendichte, g; die Ladung des Ions; und g" das elektrostatische 
Potential in der Umgebung des herausgegriffenen i-Ions ist. (k = Bo.tz- 
MANN-Konstante; JT = Abs. Temperatur.) Wegen des endlichen Ionen- 
durchmessers ist die Giiltigkeit der BOLTZMANN-Formel hier nicht selbst- 
verstandlich, weil starke Ionenanhaufungen nicht vorkommen kénnen; 
doch wird sich spater zeigen, daB Le nirgends wesentlich gréBer als 
4 ist, so daB keine schwerwiegenden Uberschreitungen der aus den 
Ionenradien sich ergebenden Maximaldichte zu befiirchten sind!. Wie 
bei DEBYE-H&cKEL soll auch hier die Exponentialfunktion entwickelt 
werden; man erhalt dann fiir die Ladungsdichte o in der Umgebung 
des 7-Ions: 


Gn, 
=) gn =—) are”, (2) 
ees ex 
j i 


da > 9;% offenbar = 0. 
j 


Aus der Porsson-Gleichung A g'? = — im o folgt somit im kugel- 
symmetrischen Fall: 
a2 (t) ; 
oe ee ) = x2 +7 gp) (3) 
mit 
gn 8 oan 
rt = 32 ae ee (4) 
€ J kT 2, ki 
j=1,2 


4 Far genauere Betrachtungen mu man die von WickE und ErGEN [Z. Elektro. 
chem. 56, 551 (1952)] hergeleitete Formel verwenden. 
2 
Zeitschrift fiir Physik. Bd. 140, 13 
| 
} 
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da hier |g,| =|q,| = Elementarladung ey. « ist die Dielektrizitats- 
konstante. 


Die Ionen der Art; kénnen sich dem Zentralion? nur bis auf die 
Entfernung a,;+ a; nahern, wenn a, die Ionenradien sind. Da die Ionen 
entgegengesetzten Vorzeichens aber den groBeren EinfluB ausiiben 
werden, diirfte es gerechtfertigt sein, fiir beide Ionenarten mit dem 
gleichen Mindestabstand a= a,-+ a, zu rechnen. Die Gl. (3) ist dann 
zu integrieren mit den Randbedingungen 


a ot) ; 
g@® +0 fir roo; “£ ) == fe, 
r=a 
Die Lésung lautet: 
gp 1 : en *(r—a) (5) 


e(1+ xa) Yr 


Die Ladungsdichten o9/") =g; |" erhalt man aus der fir r= a giiltigen 


Beziehung 


Fa 


(#) __ a j (2) 


die aus (4) nach Entwicklung der e-Funktion folgt. 

Fir geschmolzenes NaCl ist #+1,5- 10cm ~3, a=2,8- 10-8 cm, 
e=2. (Wegen der Zahlwerte vgl. Abschnitt 4.) Bei T= 1000° K, folgt 
dann!: 

1 : gy? (a) 

© 048-108 em; | 19. 
Letzteres zeigt, daB die Entwicklung der e-Funktion in (1) zwar genau 
genommen nicht gerechtfertigt war, doch diirften die Fehler wegen des 
raschen Abfalls von bei gréBeren r-Werten nicht sehr wesentlich sein, 
so daB wir auf genauere Betrachtungen verzichten wollen. Da mit 
gp ~e *'— nach (2) auch 9 ~e~*'"~® gilt und 1/x <a ist, sind die Raum- 
ladungen praktisch auf eine diinne Kugelschale der mittleren Dicke 1/x 
begrenzt, so daB wir statt (5) auch schreiben kénnen: 


(1) —_ i 1 p—*(r—a) 
iS Banaieeeod a! (Sa) 
Die im Gegensatz zu den waBrigen Elektrolytlésungen hier sehr geringe 
Ausdehnung des Raumladungsgebiets bedingt den wesentlichen Unter- 
schied in der spateren Ermittlung der Relaxationszeiten gegeniiber der 
DeByeE-HtckeEt-Theorie. 

Wir untersuchen jetzt die Deformation der Raumladungswolke bei 
Relativbewegungen von Zentralion und Umgebung. Es sei wieder o\” 


‘ Die hier auftretenden GréBen sind natiirlich nur als statistische Mittelwerte 


aufzufassen. Die Ladungsverteilung um ein spezielles Ion wird im Einzelfall ganz 
anders aussehen, 


ee a malin 
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die Ladungsdichte der j-Ionen in der Umgebung eines i-Ions. Da es 

nur auf Relativbewegungen ankommt, denken wir uns das Zentralion 

ruhend bei r=0 und die j-Ionen mit der Relativgeschwindigkeit v;(r) 

bewegt. Betrachten wir den festen Ort r, so wird sich dort im statio- 
Aol)’ 


naren Falle die Ladungsdichte nicht indern: | bi =0. Bewegen 


. . f 4 Sp / pamCObst 
wir uns dagegen mit den j-Ionen mit, so wird in jedem Augenblick 
die Ladungsdichte 9") bestrebt sein, den oben ermittelten kugelsym- 
metrischen ,,Sollwert‘‘ 0") anzunehmen, was aber nur mit der Relaxa- 


ag 


tionszeit |" vor sich geht. Fiir den totalen Differentialquotienten 


(mitbewegt) schreiben wir daher ahnlich wie bei den analogen Be- 
trachtungen in der Elektronentheorie der Metalle: 


() 2 y(t) a (i) ats). p(t) a 
dot eo! Co! < g o eo! 


] 
ah aot Da | i ri) 8 
(v; habe nur eine x-Komponente.) Daraus folgt: 


a6 
ath je ee ad 
— py Tj 7 


(t) 9) (1) 
1) Bi Oie 


(7) 


Q; 
wobei wir noch im Korrekturglied 9\ durch die ungestérte Lésung 0}’ 
ersetzt haben, was bei Beschrankung auf lineare Glieder in v; méglich ist. 
Zur Berechnung der Ionenbeweglichkeit denken wir uns ein Ion 
herausgegriffen und durch eine auBere Kraft A, mit der mittleren 
Geschwindigkeit v gegeniiber seiner Umgebung bewegt. Statt dessen 
k6nnen wir natiirlich auch, wie bei der Herleitung von (7), das Zentralion 
festhalten und seine Umgebung mit der Geschwindigkeit —v bewegen!. 
Auf das Zentralion wirkt dann die durch die CouLoMsB-Krafte bewirkte 
, Bremskraft“: 


of) ae ie acd) Go") 
4g eee Se BE spe sty gh pas de. 
B ey ey 2 ; ox # 
j=1,2 j=1,2 


wobei dt = dx dydz ist. (Das in (7) auftretende kugelsymmetrische 
Glied 9!" liefert offenbar keinen Beitrag zur Bremskraft.) Bei Ein- 
fiihrung von Polarkoordinaten laBt sich das Integral sofort auswerten 
und man erhalt wegen (6) und (5a): 
2 id AE le in aga 
i 2 ( ae! beth yal ; 
Ky = ep, ? fi. (0;" (29) — 05" (a)) i Pee, Ue | ex@RT 
q j 
1 Auf Komplikationen, die sich aus der verschiedenartigen Bewegung der beiden 


Ionenarten im elektrischen Feld ergeben, soll hier nicht eingegangen werden. 
fia 7 
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Im stationdren Falle muB die auBere Kraft K,=Kg, sein. Fur die 
Beweglichkeit 5; des i-Ions, definiert durch die Beziehung 0;=6,;- Kg 


folgt somit: 
—1 


(8) 


Die Leitfahigkeit ¢ errechnet sich dann aus den Beweglichkeiten zu 


a= > Gnb,, (9) 


i=1,2 
denn die Stromdichtej ist offenbar gegeben durch 


j= 4% = 17G;0G,E =9-E. 
Da x nach (4) durch die atomaren und makroskopischen Material- 
konstanten ausgedriickt ist, gehen in obige Gleichungen nur noch die 
Relaxationszeiten als zunachst unbekannte GréBen ein. 
Zur Bestimmung der Zahigkeit mu8 man bei einem 
Geschwindigkeitsfeld v, = const - z = ( os) -z (vgl. 
, C2 10 
Fig. 1 und 2) die Kraftwirkungen berechnen, die von 
den oberhalb der z=0-Ebene befindlichen Ionen auf 
die unterhalb der z = 0-Ebene befindlichen Ionen aus- 
getibt werden. Dazu greifen wir ein Ion der Sorte 7 
heraus, das sich gemaB Fig. 2 im Abstand € von der 
Fig. 2, Skizze zur Be- Ebene z=0 befindet. Die auf das Ion in der x- 
rechnung der Kraft- Richtung wirkende Kraft ist gegeben durch 


wirkungen auf ein Ion 
im Abstand ¢ von der 


= \? q; 0 
Ebene s=0. K;= —> fas xdt, 
— w 
j=1,2 2>0 


wobei die Integration iiber den Halbraum z> 0 zu erstrecken ist, da 
nur die dort befindlichen Ladungen fiir die Berechnung der Reibungs- 
spannung heranzuziehen sind. Fiir g{" folgt aus (7) mit der Relativ- 


. . . fal - 4 
geschwindigkeit v = (3 % (z+ ¢), wenn wir noch z+¢=2’ setzen und 


Polarkoordinaten x=r-+cos?; 2'= r+ sin? cos g einfihren: 


i) — 2) ty. OR (22). cin A: 
0" = 0) Tj * Gp OSB (5- ) r-sind:+cos@. 


Also wird 
ee 1 qi (Ov i eg” ‘ 
K,= >! (So) ,3 py 100s? Osin? i cospdrdddg. 


j=1,2 Y=r-sindcosg>c 
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Die g-Integration gibt den Faktor: 2- i a a ay Also erhalt man 


; y-sin® 
fiir obenstehendes Integral nach der Substitution cos # = u den Ausdruck 


eo” / ya 
2 f = wyr—e.l/4 sli” wdudr 
ut<1— 4 O($) 
5 ¥2\2 ; 2 
ay [ at (1— se) YA — Bards. 


Hier gibt die ¢+Integration den Faktor 2/16. Bei der tibriggebliebenen 

r-Integration, die wegen des Mindestabstandes a der Ionen von a bis co 

2 (#) 

zu erstrecken ist, kann man wegen oe ~~we-*' mit _<a die ubrigen 

14 

y-abhangigen GréBen an der Stelle y =a vor das Integral ziehen und er- 
halt dann wegen (6) und (5a) fiir K;(¢) den Wert 


fa 2 2h (i) 2 Vo 
Ke) = (2 15 95" Tj ( ree 
i(0) beans exakT 4 A a2) : 
j=1,2 


Aus diesen Einzelkraften erhalt man die Reibungsspannung, wenn man 
obigen Ausdruck iiber alle Ionen, die sich unterhalb der Flacheneinheit 
befinden, summiert bzw. integriert. Nach Division durch (@v/02), folgt 
dann fiir die Zahigkeit: 


ne eR Ds [KO at *n. 


1=1,2 
Die von 0 bis a zu erstreckende £-Integration gibt den Faktor ;4a, so 
daB man schlieBlich fiir die Zahigkeit 7 den Wert 
an 4G" 4s 


Las... auer 3 (10) 


erhalt. 
Damit ist die Berechnung von Leitfahigkeit und Zahigkeit auf die 
Ermittlung der Relaxationszeiten zuriickgefihrt. 


3. Berechnung der Relaxationszetten. 


Ein ruhendes Ion ist im statistisch-thermischen Gleichgewichts- 
zustand von der im vorigen Abschnitt berechneten Raumladungswolke 
umgeben, deren Einstellung (bzw. Verschwinden bei Entfernung des 
Zentralions) aber nicht momentan, sondern mit einer gewissen Rela- 
xationszeit erfolgt. Bei den waBrigen Lésungen ist die Relaxationszeit im 
wesentlichen durch Diffusionsvorgange bestimmt; denn denkt man sich 
das Zentralion plétzlich entladen, so werden sich die Raumladungen auf 
Grund der thermischen Diffusion der beteiligten Ionen etwa in derjenigen 
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Zeit verlaufen, welche die Ionen zur Zuriicklegung einer Strecke von — 
der GroBe des Raumladungsgebietes brauchen. Diffusionsvorgange 
kénnen aber nur wirksam werden, wenn die Reichweite des Raumladungs- 
gebiets gréBer ist, als die mittlere freie Weglange der Ionen. In unserem 
Falle erfolgt der Abfall der Raumladungen auf einer Strecke 1/x, fiir 
die bei NaCl ein Wert von 0,18-10-8cm gefunden wurde. Da diese 
GréBe natiirlich nur als statistischer Mittelwert aufzufassen ist, stellt 
ihre Kleinheit keine physikalische Unméglichkeit dar; von Diffusions- 
vorgingen kann auf solchen Strecken aber offenbar nicht mehr die Rede 
sein. An die Stelle der Diffusionsgeschwindigkeit tritt daher hier die 
mittlere thermische Geschwindigkeit der Einzelionen. Diese ist hier wie 
beim idealen Gas durch die BottzMANnN-Verteilung gegeben, da in der 
HamiILton-Funktion Raum- und Impulskoordinaten separiert werden 
kénnen. Das Raumladungsgebiet um ein Zentralion? kann man sich 
nach dem friiher Gesagten als eine diinne Kugelschale der Dicke 1/x 
vorstellen, die das Zentralion im Abstand a >> 1/x umgibt. Im Gleich- 
gewichtsfalle treten pro Flacheneinheit und sec #- u; lonen der Sorte} 
von auBen in die Raumladungszone ein, wenn u; die mittlere thermische 
Geschwindigkeit der 7-Ionen in der +7-Richtung bedeutet. Diese Be- 
ziehung gilt allerdings nur dann, wenn die Ionen beim Eintritt in das 
Raumladungsgebiet keine merkliche Potentialschwelle zu itiberwinden 


qj 


haben. Da wir aber stets - iP als klein gegen 1 vorausgesetzt haben, 


was im tiberwiegenden Teil des Raumladungsgebiets auch der Fall ist, 
kénnen wir vom EinfluB des Potentialverlaufs absehen.| Die gleiche 
Anzahl von Ionen verlaBt das Raumladungsgebiet in der umgekehrten 
Richtung. Denkt man sich jetzt die Dichte der j-lonen im Raumladungs- 
gebiet irgendwie um den Betrag 6»; iiber den Gleichgewichtswert hinaus 
erhéht, dann verlassen pro sec und Flacheneinheit iu; + dn, u; Ionen 
das Raumladungsgebiet, so daB ein resultierender Ionenstrom 6 n; > U; 


itbrig bleibt. Um diesen Betrag nimmt der Ioneniiberschu8 6n,- ; 

(pro Flacheneinheit der Kugelschale) ab, so da8 aus ‘ 
a adn, _ 
x et 

eine Relaxationszeit { 


LT Be 


u;ONn,, bn; ~ ee"! 


xu; 
fiir die j-Ionen folgt. Von der Art des Zentralions hangt die Relaxations- 


zeit nicht ab, Setzt man fiir «; den aus der BOLTZMANN-Verteilung sich 
ergebenden Wert ein, so erhalt man: 


iota 22 ym, = Ve |m 
j * VAT j 


; (11) 


wobei m; die Masse des j-Ions ist. 


AP ae @& © 
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Einsetzen dieses Ausdrucks in (8), (9) und (10) gibt schlieBlich bei 
Berticksichtigung von (4) mit |g,!=|g,| = Elementarladung e,: 
Leitfthigkeit o = 12° 492__/*7 __ (12) 
\2z )m + my 
Zahigkeit =» = 12% gg & Vt Vm (13) 


30 ea | Fag 
a—=Yr,+ 7%, = Summe der beiden Ionenradien. 
Damit sind o und » durch atomistische und makroskopische Kon- 
stanten ausgedriickt, die anderweitig bekannt sind, so daB ein Vergleich 
mit den experimentellen Ergebnissen ausgefiihrt werden kann. 


4. Der Vergleich mit den E xperimenten. 


Wir beschranken uns im folgenden auf die Alkalichloride, da hieriiber 
die meisten Messungen vorliegen und sich bei den anderen Alkalihalo- 
geniden ahnliche Verhaltnisse ergeben. In der Tabelle 1 stehen in der 
ersten Spalte die aus dem Tabellenwerk von LANDOLT-BORNSTEIN ent- 
nommenen Dichten s der Alkalichloride kurz oberhalb ihres Schmelz- 
punktes. (Die Temperaturabhangigkeit der Dichte wird in Anbetracht 
des relativ groben Modells im folgenden nicht beriicksichtigt.) In der 
zweiten Spalte stehen die aus der Dichte s berechneten mittleren Ionen- 
dichten 7,—7,—7. Die dritte Spalte gibt die Summe a der Ionen- 
radien nach GOLDSCHMIDT?. 


Tabelle 1. 
s n a %, +O, eC 
in g/cm? in cm”? in A in cem-* 
LiCl . 1,5 2,1 +107 2,59 3,6 - 10724 1,9 
NaCl 4.5 455-4077 2,79 37,0074 157 
KCl ‘5 1,2,- 10% 3,14 4,4-107%4 1,7 
RbCl 4 1,0 - 10? 3,30 5,1 -°10-%4 1,6 
CsCl 2,8 1,0 + 107? 3,46 6,4.- 10774 Ai7 


Die Dielektrizitatskonstante e wurde aus den Polarisierbarkeiten « 
der Einzelionen nach der Beziehung ¢=1 + 4207 (4, +») berechnet. In 
der vierten Spalte wurden die ~-Werte von FayAns und Joos verwendet”. 
In der letzten Spalte steht die daraus berechnete DK. 

In Fig. 3 sind die nach Gl. (12) mit obigen Zahlwerten berechneten 
Leitfahigkeiten der Alkalichloride aufgetragen. Die gestrichelten experl- 
mentellen Kurven sind dem Taschenbuch fiir Chemiker und Physiker 


1 Nach Morr u. GuRNEY: Electronic Processes in Ionic crystals. Oxford 
University Press. 
2 Fayans, K., u. G. Joos; Z. Physik 23, 1 (1924). 


190 W. OLDEKOP: 


von D’Ans-Lax entnommen?. In Fig. 4 ist die nach (13) berechnete 
Zahigkeit yon NaCl aufgetragen; die Kreuze geben Me8punkte nach 
D’Ans-Lax wieder. MeBwerte der tibrigen Alkalichloride waren dem 
Verfasser nicht bekannt. 

Die Darstellungen zeigen, daB die experimentellen Daten gréBen- 
ordnungsmaBig und in ihrer Abhangigkeit von der Zusammensetzung 
der Salze durch die Theorie befriedigend beschrieben werden, was in 
aihnlicher Weise auch fiir die iibrigen Alkalihalogenide gilt. Beriicksich- 
tigt man die Temperaturabhangigkeit der Dichte, was hier nicht getan 
wurde, so fallt bei den Leitfahigkeiten der Anstieg mit der Temperatur 


0 400 00 1200 °C 1600 0 400 800 1200°C 
a — 
Fig. 3. Leitfabigkeit der Alkalichloride, -------- Ex Fig.4. Zahigkeit von NaCl. x X xX MeBpunkte 
perimentelle Werte; — berechnet nach Gl. (12). nach D’Ans-Lax; berechnet nach Gl. (13). 


allerdings zu schwach aus, wogegen bei der Zahigkeit ein starkerer Abfall 
mit der Temperatur erhalten wird. Nach provisorischen Untersuchungen 
wird aber in beiden Fallen eine starkere Temperaturabhangigkeit im 
richtigen Sinne erzielt, wenn man die hier vorgenommene und nicht 
—% 
ganz gerechtfertigte Entwicklung der Exponentialfunktion ¢ AT bei der 
Berechnung der Raumladungswolke unterlaBt, und die BoLTzMANN- 
Formel in (1) durch die WicKE-EIGEN-Beziehung ersetzt. Dadurch wird 
eine bessere Beriicksichtigung des Potentialanstiegs im Raumladungs- 
gebiet etreicht, so daB durch die in die Leitfahigkeit und Zahigkeit dann 
eingehenden Exponentialfunktionen auch eine gewisse Annaherung an 
die Platzwechseltheorie erzielt wird ®. 

Da die zugehérigen Rechnungen aber nur noch numerisch zu be- 
waltigen sind, soll auf ihre Wiedergabe hier verzichtet werden. Bei dem 
relativ groben Modell ist mehr als eine gréBenordnungsmaBige Uberein- 
stimmung mit dem Experiment ohnehin nicht zu erwarten. Die Ver- 
mutung, daB die hier zugrundegelegten Vorstellungen tiber reine COULOMB- 


* D'Ans-Lax: Taschenbuch fiir Chemiker und Physiker. Berlin: Springer 1949. 


® Fiir wertvolle Hinweise in dieser Richtung méchte ich Herrn Professor Dr. 
EE. Wicke danken. 
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Wechselwirkungen zwischen vollig dissoziierten Ionen der Wirklichkeit 
recht nahe kommen, scheint aber berechtigt zu sein, so daB man aus 
einem eventuellen krassen Versagen der Theorie in anderen Fallen auf 
wesentlich andere Verhdltnisse bei diesen Substanzen wird schlieBen 
diirfen. In diesem Zusammenhang sei erwahnt, daB eine Anwendung 
obiger Formeln auf geschmolzenes SiO, zu vdéllig falschen Resultaten 
fiihrt. Homeopolare Bindungen und Assoziationen kénnten hierfiir ver- 
antwortlich sein. 


Herrn Professor Dr. F. SAUTER und Herrn Dr. W. HANLEIN danke 
ich fiir wertvolle Ratschlage sowie Herrn Direktor Dr. A. SIEMENs fiir 
die Méglichkeit, diese Untersuchungen durchfiihren zu kénnen. 


Rathsberg iiber Erlangen, Siemens-Schuckertwerke A.G., Laborato- 
tium ZW/LS. 


Zeitschrift fiir Physik, Bd. 140, S. 192—214 (1955). 


Wechselwirkung eines nichtrelativistischen 
Teilchens mit einem skalaren Feld 
fiir mittlere Kopplung. I’. 
Von 
G. HOHLER. 
Mit 2 Figuren im Text. 


(Eingegangen am 16. November 1954.) 


Die Energie des Grundzustandes eines nichtrelativistischen Teilchens, das mit 
beliebiger Starke an ein skalares Feld gekoppelt ist, wird mit einem neuen Variations- 
ansatz berechnet, der die Vorteile der Ansatze von PEKaR einerseits und von 
GuRaRI und LEE, Low, PINEs andererseits vereint. Er stellt eine Verallgemeinerung 
der ,,intermediate coupling method‘‘ von TomonaGa dar, bei der spezielle Korrela- 
tionen zwischen den Quanten beriicksichtigt werden. Die Anwendung auf die 
Theorie der Leitungselektronen in polaren Kristallen wird kurz diskutiert. 


1. Einleitung. 


Eine zufriedenstellende Theorie der Bewegung von Elektronen in 
polaren Kristallen hegt gegenwartig nicht vor. Die Schwierigkeiten 
riihren daher, daB eine Ubertragung der Methoden der Metalltheorie 
nicht zulassig ist, weil sie die Voraussetzung enthalten, daB die Elek- 
tron-Gitter-Wechselwirkung als kleine St6rung behandelt werden kann. 
Das mag fiir homéopolare Kristalle zutreffen, gilt aber sicher nicht fiir 
polare Kristalle, bei denen die starke Wechselwirkung des Elektrons mit 
den Polarisationswellen hinzukommt. Es miissen daher neue Methoden 
entwickelt werden, die auch bei staérkerer Kopplung anwendbar sind. 
Von den bisherigen Ansaitzen in dieser Richtung nennen wir vor allem 
die von PEKAR?, FROHLICH’, PrnEs*.5 und ihren Mitarbeitern. Weitere 
Literaturangaben finden sich in den zitierten Berichten. 

Nach klassischen Vorstellungen ruft das ruhende Elektron in seiner 
Umgebung eine Polarisation des Gitters hervor, die auch nahezu un- 
verandert mitgefiihrt wird, wenn es sich langsam bewegt. Da dieses 


' Vorlaufige Mitteilung: HéuLER, G.: Z. Naturforsch. 9a, 801 (1954). 

* PeKaR, S.1.: Fortschr. Physik 1, 367 (1954) (Bericht). — Untersuchungen 
iiber die Elektronentheorie der Kristalle. Berlin: Akademie-Verlag 1954. 

8 Frouticn, H.: Adv. Physics 3, 325 (1954). In diesem Bericht findet man auch 
eine klare Darstellung der PEKARschen Theorie sowie der Methoden zur Berechnung 
der Beweglichkeit des Elektrons in polaren Kristallen. GuRARt, M.: Phil. Mag. 44, 
329 (1953). 

4 Les, T., F. Low und D. Prngs: Phys. Rev. 90, 297 (1953), als ,, LLP“ zitiert. 
TOMONAGA, S.: Progr. theor. Phys. 2, 6 (1947). 

’ Leg, T., u. D. Pings: Phys. Rev. 92, 883 (1953). 
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Gebilde, ein Elektron und die von ihm hervorgerufene Gitterdeformation, 
auch in der quantenmechanischen Leitfahigkeitstheorie PEKARs die 
Eigenschaften eines Teilchens hat, gab er ihm den Namen ,,Polaron‘, 
der von den anderen Autoren iibernommen worden ist. 

In den obengenannten Arbeiten werden vorwiegend zwei Teil- 
probleme untersucht: 


a) Wie hangt die Energie des Grundzustandes des Polarons von der 
Starke der Kopplung ab? 

b) Wie groB ist die Tragheit eines Polarons in einem auBeren elek- 
trischen Feld ? 

Wir werden zu a) ein Variationsverfahren angeben, das bessere Er- 
gebnisse liefert als die bisherigen Ansaitze. Die Anwendung auf b) soll 
in einer spateren Arbeit erfolgen. 

Da die Starke der Elektron-Gitter-Wechselwirkung erhebliche 
Schwierigkeiten mit sich bringt, muB das Modell soweit wie méglich 
vereinfacht werden. a) Reduziert sich dann auf ein fundamentales 
Problem der nichtrelativistischen Feldtheorie, namlich auf die Bestim- 
mung der Selbstenergie eines Teilchens, das mit einem skalaren Feld in 
Wechselwirkung steht, als Funktion der Kopplungsstarke. Fiir die Leit- 
fahigkeitstheorie braucht man dann weitere Eigenschaften des ,,ange- 
zogenen Teilchens**, das hier als Polaron bezeichnet worden ist, insbe- 
sondere seine Tragheit. AuSerdem treten durch die Einfiihrung des 
Begriffs der Temperatur Fragestellungen auf, die bisher in der Feld- 
theorie nicht behandelt worden sind. 


2. Der HAMILTON-Oferator. 


Da vor kurzem zwei Berichte iiber Elektronen in polaren Kristallen 
erschienen sind?., verzichten wir an dieser Stelle auf eine Begriindung 
des HamILTon-Operators. Wir wahlen die tibliche Form 


ee ¥ Be Dts b,— iho 4 7) 8D | (b, e't® — be eit"), (2.1) 
m 


rt, p sind Koordinate und Impuls des Elektrons, (0¢ + br) und (by — be) /i 
verallgemeinerte Koordinaten und Impulse des Gitters bei Zerlegung 
der longitudinalen optischen Gitterschwingungen nach laufenden 
Wellen, f ist der Ausbreitungsvektor, w die als konstant angenommene 
longitudinale aah ores Pre era V ist das Normierungsvolumen, 
u-t=(2mw/h)~* eine Lange der GréBenordnung 1077 cm. Die atomisti- 
sche Struktur des Gitters wird lediglich dadurch beriicksichtigt, daB die 
Summen im Raum nur iiber eine Kugel vom Radius ky erstreckt 
werden. Die Kopplungskonstante g? (FROHLICHs «) kann berechnet 
werden, wenn neben Materialkonstanten des Kristalls auch noch die 
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effektive Masse m des Elektrons in dem periodischen Potential des — 
Gitters bekannt ist (g? ~|\m m). Von dieser kann man aber gegenwartig 
mit Sicherheit nicht viel mehr als die GréBenordnung angeben, so daB 
tes Zahlenwerte fiir g? nicht vorliegen. FROHLICH schatzt fiir NaCl 
~6, PEKAR fithrt ein (nicht ganz iiberzeugendes) Argument fir 
28,8 an, indem er aus Messungen an F-Zentren auf m= 2,78 Elek- 
tronenmassen schlieBt (vgl. S. 181 seines Buches’). 
Wir gehen zu dimensionslosen GréBen iiber®: 


t>t/u, prhup, tout, V>V/u, H->hiod. 


Dann wird aus (2.1) 
H= p+ 5) by by — tx > ; (bp e§** — db, e=***). (2.2) 


Den haufig vorkommenden Ausdruck (42 g?/V)! haben wir mit x ab- 
gekiirzt. Fiir das Elektron wahlen wir die Ortsdarstellung; p? wirkt 
dann als Differentialoperator —A. b; und >, sind Erzeugungs- bzw. 
Vernichtungsoperatoren fiir Schallquanten der Wellenzahl f. ®, sei der 
, Vakuumzustand™, also b,®,—0 fiir alle —. Aus ihm kénnen wir in 
bekannter Weise den allgemeinsten Zustand aufbauen. 

Wir wollen versuchen, eine Naherung fiir den Grundzustand des 
HAMILTON-Operators (2.2) zu finden, die fiir alle g* brauchbar ist. Dazu 
benutzen wir das Variationsverfahren. Wir machen also fiir den Grund- 
zustand einen Ansatz, der noch unbestimmte Funktionen oder Parameter 
enthalt, berechnen den Erwartungswert von H und wahlen dann bei 
jedem g? diejenigen Funktionen oder Parameter, die ihn zum Minimum 
machen. — In 3. besprechen wir den fiir schwache Kopplung guten An- 
satz von LEE, Low und PIngs, in 5. den fiir starkere Kopplung bisher 
besten Ansatz von PEKAR, weil wir spater zeigen wollen, da8 unser Ver- 
fahren zu beiden in enger Beziehung steht. 


3. Der Ansatz von LEE, Low und Pines 4, 
Der Operator der Gesanitwellenzahl 

W= p+ Voy (3.4) 
ist mit 7 vertauschbar®. Daher existiert ein vollstandiges System simul- 
taner Eigenzustinde der beiden Operatoren. Geben wir einen Eigen- 
wert WW’ von YW vor, so ist damit ein Unterraum des HILBERT-Raumes 
festgelegt, der noch viele Eigenzustande von H enthalt, darunter hat 
einer den tiefsten Eigenwert. Ziel der Rechnung ist es, eine gute Nahe- 


rung fiir diesen Zustand und seinen Eigenwert E als Funktion von QW’ 
und g? zu gewinnen, 


* Der Unterschied gegeniiber dem Erhaltungssatz fiir den paige stg 2 wird 
von FROHLIcH® (§3) eingehend diskutiert. 


Ane ti, St 


> o- 
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Fur kleine Kopplung gehért der Grundzustand des Systems zu W’ = 0 
genauer (0, 0,0). Wir nehmen an, ohne daB wir einen Beweis dafiir 
haben, dies gelte fiir alle Kopplungsstairken. Wenn der Grundzustand 
zu einem YW’ =: 0 gehérte’, wire er stark entartet, denn E hangt nur von 
(W’)? ab, weil in H keine Richtung ausgezeichnet ist. 

Die Tatsache, daB wir den tiefsten Zustand von H jeweils nur in dem 
durch einen Eigenwert von W festgelegten Unterraum suchen, ermég- 
licht es, die Operatoren des Elektrons zu eliminieren. Dazu transfor- 
mieren LLP unitar mit 


oe 


S = ef WB Leoy bh) (3.2) 


Dann sind die tiefsten Eigenwerte von SHS in dem durch p’=0 
(wegen S71 WS=W’+p) festgelegten Unterraum zu ermitteln. Fiir 
Zustande aus diesem ist aber S1HS gegeben durch den Operator 


+ + 7. 
H = (B'— Yb) + Vo te, (3.3) 


der nur noch von den 6, und 5; abhangt. Der Variationsansatz von LLP 
lautet 
Dip = ZO POW) g,. (3.4) 


Wenn wir uns nun auf Y’=0 beschranken, wird der Erwartungswert 
von H mit dem normierten Zustand (3.4) ein Minimum, wenn 


() =i (3.5) 


gewahlt wird. Man erhalt fiir die Energie des Grundzustandes 
Evrp = — 8°. (3.6) 


7 In einer soeben erschienenen Arbeit [Progr. theor. Phys. Il, 449 (1954)] 
kommt E. Haca zu dem Ergebnis, daB fiir g? > 6 die tiefste Energie nicht mehr 
zu YW’ =0 gehért. Wir vermuten, daB er den Giiltigkeitsbereich seiner Naherung 
iiberschatzt hat, denn die nach unserer Meinung zuverlassigere Storungstheorie 
von LLP gibt bei g?=6 nur eine Korrektur von etwa 20%. — Im Anschlu8 an 
Yoxora versucht Haca (S. 450, 455), Resultate von Gross [Phys. Rev. 84, 818 
(1951)] fiir das System 1 Elektron + 1 Oszillator unmittelbar auf unser System 
zu iibertragen. Das ist nicht zulassig, denn die beiden Systeme unterscheiden sich 
wesentlich voneinander. Von vielen Argumenten hierfiir bringen wir nur zwei: die 
Variationsverfahren?,4,® kénnen nur fiir unser System, nicht jedoch fiir das von 
Gross, formuliert werden, weil dazu viele Oszillatoren mit verschiedenen f-Werten 
notig sind; bei uns gibt der tiefste Eigenwert von H, (in 4.) eine obere Schranke, 
im Spezialfall eines Oszillators gibt er jedoch eine untere Schranke fiir die Energie 
des Grundzustandes bei ¥8’=0. — Im Falle der von Haca ebenfalls erwahnten 
Supraleitung hangt die Starke der Kopplung mit der Anwesenheit sehr vieler 
Elektronen zusammen, so da®S hier ein unmittelbarer Vergleich ebenfalls frag- 


wiirdig ist. 
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Eine andere Form des Ansatzes folgt mit der Beziehung 


gfbt—f*b _ g—hif* gfb* e—f*> (3.7) 


Man erhalt 
Opp =e ZN JT POND, 
t 


=e EI A) OF EDT Ht) F(t) 85 8%, + -f Do 
3 7 


f,, €. 


(3.8) 


Wenn wir dies mit der Entwicklung eines beliebigen Zustandes nach Zu- 


standen mit 0, 1, 2,... Teilchen vergleichen 
B= {ey + Delt) + Dearth t) +f 3.9) 
t vfs 


kénnen wir die Mannigfaltigkeit der von LLP beriicksichtigten Zustande 
so charakterisieren: es sind alle Funktionen g, (f,, ..., £,) zugelassen, die 
als Produkt in der Form f(f,) /(f,) ...f(£,) darstellbar sind, wobei die 
Funktion f fiir alle » die gleiche ist und von vornherein beim m-ten Term 
ein Koeffizient (7!)-! steht. Man kann von einem HARTREEschen Ansatz 
(im f-Raum) sprechen, bei dem sich alle Schallquanten in demselben 
Zustand /(f) befinden. Diese Formulierung zeigt die Verbindung zu 
TomonaGas ,,intermediate coupling method‘. 


4, Der Riickstopterm. 


Wir wollen nun den RiickstoBterm in (3.3) [das ist () f by })2, denn 
vorlaufig betrachten wir nur W’ = 0) zu einer oberen und einer unteren 
Schranke® fiir den tiefsten Eigenwert von H in Beziehung setzen. Der 
tiefste Eigenwert von H liegt sicher héher als der von 


é 71 + q ae b om bf é 
H ~ (S015 b)P = VKH + HD I ; — (4.4) 
mit 
Cob +1 = (4.2) 
entsteht aus (4.1) 
Tike =, | 
2 — 8D (4.3) 


Der tiefste Eigenwert dieses Operators hangt wesentlich von ky ab. 


> 
to 
> 
—) 
3 
n 


gt. (4.4) 


Die Existenz dieser unteren Schranke ist also eine Folge des Abschnei- 
dens (s. u. 6.). 
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Nach dem gleichen Verfahren, jedoch mit 


= db +4 * 
kann man eine exakte Lésung auch noch erhalten, wenn man nur einen 
Teil des RiickstoBterms wegliBt. Wir ordnen die Operatoren um nach 
dem Gesichtspunkt: Vernichtungsoperatoren nach rechts: 


(Std; br) = ZAG be +S (tas te) BF BE, Oy be (4.6) 
und zerlegen ‘H in ‘HH, und H’ 
Ho =D, (1 +b) Fe +x DI a (4.7) 
HH’ = >” (ty, fe) Bg, Bf, be, by. 
Der tiefste Eigenwert von H, ist§ 
E=— 2. (4.8) 


Das ist eine obere Schranke, denn es stimmt genau mit dem Variations- 
ergebnis von LLP iiberein (3.6). 

Die exakte Kurve beginnt bei g?=0 an der oberen Schranke. Um 
dies einzusehen, gehen wir von dem Operator ‘H, aus und behandeln ‘H’ 
als kleine St6rung. In der zweiten Naherung finden wir eine negative 
Korrektur zu (4.8), die mit g* wachst und bei g?=5 erst etwa 6% be- 
tragt*. — Die gewéhnliche St6rungsrechnung mit dem Kopplungsterm 
als St6rung kann, da die erste Naherung keinen Beitrag gibt, nur auf 
(4.8) fiihren, denn das Ergebnis muB proportional zu g? sein und wir 
haben gesehen, daB (4.8) fiir kleine Kopplung giiltig ist. Die Zerlegung 
H =H, +H’ ist jedoch besser, weil bei einer Entwicklung des Grund- 
zustandes in der Form (3.9) im Falle schwacher Kopplung die beiden 
ersten Terme dominieren. H’ wirkt auf diese aber wie ein Faktor Null. 
AuBerdem bewirkt der in H’ enthaltene Faktor (f,, f,) eine Verkleinerung 
des Matrixelements (bei g?=5 ist die mittlere Zahl der Schallquanten* 
bereits 2,5). — Mit ansteigender Kopplung nimmt der Einflu8 von H’ 
dann aber schnell zu. 

Eine Verbesserung gegeniiber LLP ist nur dadurch méglich, dab He 
Variationsansatz so gewahlt wird, dab 


a) sein Erwartungswert mit H’ negativ ist, 
b) sein Erwartungswert mit H, ,der sicher hoher als — g? liegt, die 
Verbesserung a) nicht wieder aufhebt. 


8 Wir haben bei der Berechnung von (4.8) wie LLP die f{-Integration bis co 
erstreckt. Beriicksichtigung des Abschneidens gabe den Faktor (2/71) arc ty hy 0,9. 
9 Siehe auch Haaa, l.c. ’. 
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Der Erwartungswert von H’ mit ®,;p ist Null. Da fiir kleine Kopp- 
lung ®,;» schon recht gut ist, wird man nicht einen ganz neuen Ansatz 
suchen, sondern eine Verallgemeinerung, die auf negative Erwartungs- 
werte von H’ fiihrt. Es bieten sich zwei einfache Méglichkeiten: 

a) Wir bleiben bei der HARTREE-Naherung im fRaum, lassen aber 
neben f(f), das sich fiir Y@’—=0 als nur von |f| abhangig erwiesen hat, 
noch weitere von der Richtung von f abhangige Zustande zu, damit der 
Faktor (f,, £,) in H’ zur Wirkung kommen kann. Diesen Weg haben LEE 
und Pings eingeschlagen®, indem sie neben dem ,,s-Zustand“ /(&) noch 
die drei ,,p-Zustande“ g(k) k,/k, g(k)k,/k, g(k) k,/k beriicksichtigen. 

8) Wir verlassen den Rahmen der HArtTREE-Naherung, indem wir 
jedem der n-Teilchen-Terme in (3.8) einen Faktor anfiigen, der nicht 
mehr als Produkt von Funktionen mit je einer Variablen f; geschrieben 
werden kann. Das entspricht der Einfiithrung von Korrelationen in den 
Variationsansatzen fiir das H,-Molekiil und das He-Atom, jedoch auf 
den f-Raum iibertragen. Ein Beispiel ware (wir beschranken uns auf gy) 


82(f1, f) = x f (Ry) f(Re) [1 + 0? (f + €)7)*. (4.9) 


Man sieht leicht, daB die Paare von Vektoren mit negativem (f,, f,) 
bevorzugt werden, denn sie haben ein kleineres (f, + f,)*. Mit dem Wert 0 
des Variationsparameters « kommen wir auf ®;;p zuriick, so daB das 
Ergebnis nicht schlechter sein kann als das von LLP. 

Der Hinweis auf diese naheliegende Verallgemeinerung ist natiirlich 
erst interessant, wenn es gelingt, einen Weg zur Berechnung des Er- 
wartungswertes von H zu finden. Durch Ubergang zu einer FouRIER- 
Integraldarstellung fiir den hinzugefiigten Faktor kommt man zu einem 
Produktansatz zuriick, kann analog zu (3.8) einen geschlossenen Aus- 
druck fiir den neuen Ansatz finden und auch den Erwartungswert 
von ‘H berechnen. Wir fiihren das nicht naher aus, weil sich in 10. der 
Zusammenhang mit einer anderen Methode zeigen wird, der wir uns 
nun zuwenden wollen. 

Es gibt einen einfachen Ansatz von PEKar, der fiir g? > 10 bald 
erheblich besser wird als der von LEE und Pings. Die vorliegende Arbeit 
ist aus der Vermutung entstanden, daB diese Verbesserung auf Korrela- 
tionen in obigem Sinne zuriickzufiihren ist, die in PEKARs Ansatz ent- 
halten sind. Wir werden zunachst den Ansatz besprechen (5.), dann 
eine Modifikation angeben, durch die es méglich wird, die Form (3.9) 
herzustellen (9.) und zeigen, daB die Vermutung zutrifft (10.) Wahrend 
bei PeKar die Starke der Korrelation vorgegeben ist, regeln wir sie durch 
einen Variationsparameter und geben ihr damit bei jeder Kopplung den 
giinstigsten Wert. In 10. werden wir schlieBlich feststellen, daB der 
Ansatz verallgemeinert werden kann, ohne daB neue Schwierigkeiten 
bei der Berechnung des Erwartungswertes von H auftreten, 
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3. Der Ansatz von PEKAR", 2,3, 


PEKAR kommt, von klassischen Uberlegungen ausgehend, zu der 
Vorstellung, daB man in guter Naherung den Grundzustand des Polarons 
beschreiben kann, indem man fiir das Elektron und fiir das Gitter je 
- einen ,,besten“ Zustand y(t) bzw. ®, angibt. Diese beiden Zustiande 
sollen in folgender Beziehung stehen: 

a) Die aus dem Elektronenzustand y(r) zu berechnende Ladungs- 
_ dichte ey*y erzeugt durch Polarisation den Gitterzustand @,. 

b) Im Zustand @, sind im Gitter Ladungen aus ihren Gleichgewichts- 
lagen verschoben. Man kann in der Kontinuumsnaherung auf einfache 
Weise das als Folge der Verzerrung auftretende elektrische Potential 
ausrechnen. y(r) soll der tiefste Elektronenzustand in diesem Potential- 
topf sein. 

Die Forderungen a) und b) sind gleichbedeutend damit, daB man zur 
Bestimmung des Grundzustandes von H den Ansatz 


YO = p(t) D (5.1) 


macht. Es ist giinstig, y(r) so zu wahlen, daB die Ladung des Elektrons 
praktisch auf einen kleinen Bereich (von der GréBe einiger Gitterkon- 
stanten) lokalisiert ist. Obwohl physikalisch alle Punkte des Kristalls 
gleichberechtigt sind, zwingt der Ansatz zur Auszeichnung eines Ortes, 
ohne daB das Ergebnis von seiner Lage abhangt. Wir wahlen den Ur- 
sprung des Koordinatensystems, weisen durch den oberen Index an ¥;” 
darauf hin und setzen y(r) kugelsymmetrisch an, weil in H keine Rich- 
tung bevorzugt ist. — Die Variation von (y®,, Hy@,) nach @, bei 
festem y und Beachtung der Normierung fiihrt auf eine exakt lésbare 
Eigenwertgleichung (Feld mit vorgegebenen Quellen, ,,BLocH-NorRD- 
stECK-Problem‘) 
(y, Hy) Pg = Ely) ®. (5.2) 
Mit der Abkiirzung 
oft) =fe|p2dr=o((t), da y)=y(r) 6.3) 


entsteht ein HaAmitton-Operator fiir Oszillatoren mit verschobener 
Ruhelage 
Beebe 
+ o(h) “t t 
(y, Hy) =—(y,4y) + oa +e (5.4) 
Wir bringen den in 4, b* linearen Term wie in (4.2) zum Verschwinden 
durch 
. k 
oy = by tin OO) (5,5) 
10 Wir betrachten hier nur den Produktansatz von PeKaR. Zur adiabatischen 
Naherung s. !”. 
Zeitschrift fiir Physik. Bd. 140. 14 
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und erhalten 


(y, Hy) =a Por ye — (yp, Ay). (5-6) 


Der tiefste Eigenwert dieses Operators ist leicht abzulesen, weil die cy 
denselben Vertauschungsrelationen geniigen wie die b,. AuBerdem gehen 
wir von der Summe zum Integral iiber 
ky oo 
Z F : 4 cy \2 
Eply]= —¢2 | o?(k) dk + 4a || a) de (5.7) 
0 0 

Bei der Bestimmung des Grundzustandes von (5.4) beachten wir, daB 
die verschiedenen f in (5.4) separiert sind. Wir machen daher den 


Produktansatz »3) 
a ire ') 

Zur Bestimmung der C,,(f) benutzen wir die Relation ¢®,=0, durch 
die der gesuchte Zustand definiert ist. Mit (5.5) und (5.8) folgt das 
Gleichungssystem 


Jn +1C,..() = e C,, (8) (5.9) 
mit der Lésung 
Ca(f) = nt (x 2) Sol (5.10) 
(n!)? 


Cy ist aus der Normierung zu bestimmen. Zu jedem y(r) lautet also der 
bestmégliche normierte Gitterzustand 


_i,.y? th 5 .23> so S @ p+ 


b —ix by +-by) 
t (os + 
@, =e 2° 4k» ak @,=e hat BP 


®,. (5.41) 


Bemerkenswert ist die Ahnlichkeit zu (3.4), auf die wir spater zuriick- 
kommen. 

Die Bestimmung des giinstigsten y(r) wird zweckmaBigerweise so 
ausgefiihrt, daB man eine geeignete normierte Funktion mit einigen 
Variationsparametern einsetzt und diese so bestimmt, daB der Er- 
wartungswert von H ein Minimum wird. Wir beschranken uns auf zwei 


einfache Ansitze mit je einem Parameter und vernachlassigen zuniichst 
das Abschneiden bei kp. 


voy=teay tens; oy alr 4 (4) )* 
(5.42) 12 
Ep(g*) = — H¥ng* = — 0,008 gt; a= 18g 2, 
r B? k* 
v0) =(2iape *; o(k)=e 8 | 
E,(e2) = a Ge hi a , (5-43) 
ga Bede dy = = — 0,106 g*; B=3\xg | 


‘Dieser Ansatz wird auch von FrRértLtcr diskutiert, 
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PEKAR kam mit einigem Aufwand an numerischer Rechnung noch etwas 


- weiter? 


y (r) =0,1229 a3 (1 + v/a + 0,4516 r2/a2) e> « 
( [% = 516 r?/a?) e (5.14) 


E p(g?) = — 0,109 g*. 


Dieser Wert liegt fiir g?> 9,2 tiefer als das Ergebnis von LLP und fiir 
g° > 10,2 tiefer als das von LEE und Pines (Fig. 2). 


6. Der Grenzfall starker Kopplung. 


LEE und Pings haben fiir ihren An- 
satz den Grenzfall starker Kopplung 
diskutiert. Sie finden, wenn sie bei 
ky abschneiden, fiir das Variations- 
ergebnis E(g?) eine Annaherung an 
die untere Schranke und damit an 
die exakte Lésung. PEKarR erhielt ,- ¥ 
ohne abzuschneiden E ~g* (5.14). Da | # 
die untere Schranke wie g? geht (4.4), ih 
liegt die Vermutung nahe, daB mit 
Beriicksichtigung. des Abschneidens 
ein ahnliches Verhalten eintritt. Die 
Rechnung macht mit den Ansadtzen 
(5.12) und (5.13) keine Schwierig- 
keiten! und bestatigt die Erwartung. 4 
Obwohl LEE und PrNes nicht den ™ 


200 


r Fig. 1. £(g*) in doppelt-logarithmischer Dar- 
penzen V erlauf i sondern nur das stellung. ——_— Mit Beriicksichtigung des Ab- 
asymptotische Verhalten angegeben _ schneidens bei ky =7,5; ------- mit ky = 00. Ler 


und Pines *® haben nur einige Punkte (x) mit 
haben, darf man wohl annehmen, daB k,=oo gerechnet und gezeigt, daB sich ihre 


jhre Kurve hoher liegt als unsere Kurve fiir g*> co der unteren Schranke (4.4) 
(Fig. 1). asymptotisch nahert. 

Die in (5.13) vorkommende GrdéBe f ist ein MaB fiir die Ausdehnung 
des Potentialtopfes und damit des Gebietes, in dem das Gitter merklich 
polarisiert ist. Es nimmt bei kleinen und mittleren Kopplungen wie 
eg ab (5.13). Bei g?~15 erreicht es die GréBenordnung der Gitter- 
konstanten. Nun beginnt das Abschneiden merkbar zu werden, indem 
es zu einer langsameren Abnahme fiihrt. Die kinetische Energie des 
Elektrons, die ja unmittelbar mit der GréBe des Potentialtopfes zu- 
sammenhangt, verhalt sich entsprechend. Die Tatsache, daB sie bei 
starker Kopplung gegeniiber den anderen Energien relativ klein wird, 
daB also ein einfacher ,,Grenzfall starker Kopplung“ existiert, ist somit 
ein Abschneideeffekt und hat mit unserem eigentlichen Problem nichts 


VAvan abbr 
14* 
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Die Frage, ob es auch ohne Abschneiden einen einfachen Grenz- 
fall starker Kopplung gibt, behandeln wir im zweiten Teil der Arbeit. 


Wir werden im folgenden alle £-Integrationen bis nach oo fihren, 
weil sonst erhebliche rechnerische Komplikationen auftreten. Die Er- 
gebnisse sollten bis etwa g? 15 und damit fiir alle bisher untersuchten 
Kristalle anwendbar sein, denn sowohl bei PEKARs Lésung als auch bei 
der von LEE und Pines (° S. 888, abgesehen von §) wirkt sich das Ab- 
schneiden erst von dieser Kopplung an aus. 


7. Der Translationsoperator des Polarons. 


Wir haben oben gesehen, daB in PEKARs Ansatz ein Ort ausgezeichnet 
ist: der Mittelpunkt der Ladungsverteilung und zugleich auch der Gitter- 
deformation. Fiir diesen ,,Ort des Polarons“ hatten wir den Ursprung 
des Koordinatensystems gewahlt. Mit gegebenem Zustand des Elektrons 
und dem zugehérigen besten Gitterzustand ergab sich fiir den Zustand 
des Systems 

lay wi ~ix) a”) at 


P=ypr)e 7H Ke * Dy. (7.4) 


Darin ist @ aus py zu berechnen (5.3). 

Wir interessieren uns nun fiir den Fall, da8 sich das Polaron an einem 
anderen Ort befindet, etwa in rg. Die Rechnung verlauft wie in §., 
an Stelle von p(r) haben wir die gleiche Funktion von |r— ry|, an Stelle 
von o(h) 


Or, (f) = f ew" | (|e — to) Pde =e"! o(k). (7.2) 


Bei der Bildung von (5.4) erhalt } jetzt noch den Faktor e‘** und 6* 
den Faktor e~''", Der Zustand des Polarons am Ort t, lautet schlieBlich 


yr el’) 


p= kt 


1 ~ 
“i ain yer gate 
(t) : ; ¥. % . x? ix bp é 
PE = yp(|r— tl) O™; Gee 2 i 2 ®,. (7.3) 


Der. obere Index an ®) bezeichnet den Mittelpunkt der Gitter- 
deformation, den unteren Index ,,G‘‘ lassen wir ven jetzt ab fort. 


Nun konstruieren wir einen unitaéren Operator T(tg), der die Verschie- 
bung bewirkt 


Pp = T (tq) PE. (7.4) 
Der Elektronenanteil in T muB offensichtlich der Verschiebungsoperator 
ent (toP) — g—toarad (7.5) 

sein. Der Gitteranteil 7; ist 


Ti; (tp) BO = eZ Mtt) Goo) _ Hie) (7.6) 
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Zum Beweis dieser Formel benétigt man zwei Zwischenergebnisse der 
zu (3.3) fihrenden Rechnung 


TG (to) be Te" (tq) = by e~ tthe tia) 
und 


TG (to) Py = Dy. (7.8) 


Der ,,Translationsoperator des Polarons‘‘ ist also! 
T(t») = en ite (P+ Eee) (7.9) 


Zur infinitesimalen Translation gehért der schon friiher eingefiihrte 
Wellenzahloperator (3.1) 


W= p+ DEH by. (7.10) 


Die zu der unitaéren Transformation T gehérige kanonische Transfor- 
mation lautet 


tT>t—t, pp, Desabeett, bf > df eth, CAG) 


Die Invarianz von H gegeniiber (7.11) entspricht der Vertauschbarkeit 
mit T bzw. W. In einem richtigen Gitter, also ohne die Kontinuums- 
naherung, ist die Situation anders, weil der obige Ubergang zur infini- 
tesimalen Transformation nicht mehr méglich ist. 


8. Erweiterung des PEKARschen Ansatzes. 


Man kénnte denken, daB es bei mittlerer und erst recht bei starker 
Kopplung nicht mehr sinnvoll ist, von den Teilsystemen ,,Elektron‘ 
und ,,Gitter“ sowie von ihren Zustanden zu sprechen. PEKAR tut dies 
jedoch in demselben Sinne, in dem HARTREE Zustande einzelner Elek- 
tronen behandelt. Das Elektron stellt seinen Zustand nach dem M7ttel- 
wert des Gitterzustandes ein und das Gitter richtet sich nach der muittleren 
Ladungsverteilung des Elektrons. Es ist bemerkenswert, daB PEKARs 
Ansatz gerade bei kleiner Kopplung versagt, aber bei starker Kopplung 
bis jetzt giinstiger war als jeder andere. 


Wenn man ihn verbessern will, muB man iiber den Produktansatz 
(5.1) hinausgehen, d.h. Elektronen- und Gitterbewegung diirfen nicht 
mehr unabhangig voneinander sein. Die folgende einfache Abhangigkeit 
kann ohne besondere Schwierigkeit beriicksichtigt werden!: der Zustand 
des Elektrons sei weiterhin y(r), aber der Mittelpunkt des Potential- 
topfes liege jetzt nicht mehr im Nullpunkt sondern bei 7r. Mit dem 


12 Die Bezeichnung stammt von PEKAR, er meint die Translation des lokalisierten 
Zustandes Y) und nicht etwa die Anwendung von T auf exakte Eigenzustande. 
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Gitteranteil (7.6) des Translationsoperators folgt’® 


> e—inte 


. 59) hes 
= Te (nr) =yre PHBE “HR  @,  (8A) 


y ist ein Variationsparameter, der die Starke dieser ,,Korrelation‘‘!* 
regelt. Mit dem Wert Null kommen wir auf den alten Ansatz zuriick. 
Wir berechnen den Erwartungswert von H 


(HY) =(YP,AYP); H=Te"(yt)HTe(yr). (8.2) 
H ergibt sich mit (7.7) zu 
H = — (grad — in >it by by)? + 


; Poe er 8. 
+ or bb — ix) ; (b, FtEC—D _ By eFFE—m) (8.3) 
Insgesamt folgt fiir den Erwartungswert der noch von der Funktion yp (7) 
abhangige Ausdruck 


0 


x (8.4) 


E(;¥) =4n { (2) rar— | 
af 2fe (k)o(R(1—m)) ak — fo dk— m[ornanl. | 
0 


0 0 


In dieser Formel gilt zwischen 0 und y die Beziehung (5.3). Andererseits 
kann man natiirlich einen Ansatz (8.1) machen, in dem o und w unab- 
hangig voneinander sind. Dann geht zwar die anschauliche Inter- 
pretation verloren, aber (8.4) bleibt giiltig, wir miissen nur die Funktion g 
mit dem Argument (1 — 7) durch einen Index von den anderen unter- 
scheiden, denn sie ist nach wie vor gem&B (5.3) aus y zu berechnen, — 
Nun kommen wir auf die nach (5.11) bemerkte Ahnlichkeit zu (3.8) 
zuriick. (8.1) geht in S®,;p [(3.2), (3.8), W’ = 0] tiber, wenn p(r)—14, 
—txo(hk)/k->f(k), n—>1. Ein geeigneter Variationsansatz (8.1) mit 
voneinander unabhangigen 0 und y ist also sicher nicht schlechter als 
der von LLP und der von Pekar. Da jedoch eine vorlaufige Auswertung 
nicht so gute Ergebnisse brachte wie das nun zu schildernde Verfahren, 
das tiberdies den Vorzug hat, Eigenzustande zu YW zu liefern, gehen wir 
hier auf spezielle Ansatze (8.1) nicht naher ein. 


18 Herr Professor FROHLICH hat mir freundlicherweise die in $6.2 seines Berichtes 
erwahnte Arbeit von BucKINGHAM (E.R.A. Report L/T 312) iibersandt und mich 
darauf aufmerksam gemacht, daB BucKINGHAMs Ansatz mit (8.1) eng verwandt ist, 

M Der Begriff ,,Korrelation’’ wird hier in anderem Sinne verwendet als in 4. 


— 


; ‘ : ae nae 
Wechselwirkung eines Teilchens mit einem skalaren Feld. I. 205 


9. Herstellung von Eigenzustdénden von YW aus PEKARs Produktlisung. 
In einer Bemerkung am Ende von 4. haben wir darauf hingewiesen, 


daB es interessant ist, den PEKARschen Ansatz mit dem von LLP zu 
vergleichen. Unmittelbar ist das nicht moéglich, denn es fehlt ihm eine 


_ Eigenschaft, die sowohl die exakte Lésung als auch die von LLP und 


Ler und Pines haben: er ist kein Eigenzustand von %, Man kann 
PEKARs Ansatz aber nach Eigenzustinden von 8 zerlegen und den zu 
W'= 0 gehdrigen Anteil % nehmen. Dann ist der Vergleich méglich 
und er fiihrt auch zum Ziel, denn wir werden sehen, daB Y, noch besser 
ist als Y. Im Hinblick auf die spatere Berechnung der Masse des Pola- 
rons beschranken wir uns vorlaufig noch nicht auf YW’ = 0 sondern be- 
handeln gleich den allgemeinen Fall. 

Zur Zerlegung von ¥, benutzen wir einen Projektionsoperator, der 
auf den zu W’ gehdrigen Unterraum projiziert 


Poy = (277)9 6 (W’ — W) YO = f iWWRIH YO, (9.1) 


Yay ist nicht normiert. — Um zu zeigen, daB es Eigenzustand von W 
ist, wenden wir zunachst den Operator T (7.9) an 


T (t) Pou = e 2 B% f c(B-WRARWO — — WH yy, (9,2 


Das gilt fiir alle ry, also ist Yq Eigenzustand von WW zum Eigenwert YW’. 
(9.1) kann auch anders gelesen werden 


Fe [ATR PP AR = Lar AP OR, (9.3) 


~~ 


Wir haben also Zustande des PEKARschen Typs mit verschiedenen Orten 
des Polarons so linear kombiniert, daB das Resultat Eigenzustand von YW 
ist. Diese Formulierung erinnert an die Metalltheorie. Bei seiner Nahe- 
rung fiir tiefe Energien ging BLocH von den Atomfunktionen, also 
lokalisierten Zustanden, aus und bildete eine analoge Linearkombination, 
um die Symmetrie des HAmiLton-Operators zu beriicksichtigen. Im 
Sinne der Stérungstheorie fiir entartete Systeme ging er zur ,,richtigen”’ 
Linearkombination nullter Naherung iiber. Der Erfolg war ein ver- 
besserter Wert fiir die Energie des tiefsten Zustandes. Auch in unserm 
Fall ist der Erwartungswert der Energie mit Y%, tiefer als das PEKARsche 
Ergebnis (Fig. 2). 
Wir setzen den Erwartungswert der Energie in unserer Naherung 
(9) — 2 Pw) (9.4) 
(Pay, Pay) 
15 FROHLICH weist in seinem Bericht* ausdriicklich darauf hin, dab es vorteil- 
haft ware, die PeKarsche Lésung so zu modifizieren, daB Eigenzustande von Xx 


entstehen (S. 352). Diese Méglichkeit wird auch von H., HAKEN diskutiert (Halb- 
leiterprobleme I, herausgegeben von W. SCHOTTKY, Vieweg 1954, S. 72, FuBnote 8). 
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in Beziehung zu PEKARs Resultat Ep 


f (Poy, H oy) dB’ = fe A-™) (WHY) dw’ dRdW 


(9.5) 
= (22)? V(W”, H Pp”) = (22)? VEp. 
Es folgt 
vy (‘Pay Pay’) Sit 
Ep = [£(B ea a B’. (9.6) 


PrKARs Energie des Grundzustandes ist also ein Mittelwert tber die 
tiefsten Energien bei gegebenem 2’ im Rahmen unseres Ansatzes (9.1). 
Sie liegt damit sicher Adher als unser tiefstes E (Y’). Wir haben ebenso- 
wenig wie bei der exakten Lésung einen allgemeinen Beweis dafiir, daB 
dieses tiefste E zu YW’ =O gehdért. Die notwendige Bedingung fir ein 
Minimum gradgy E (Y8’) =0 ist jedoch erfiillt, weil E nur von WW’? ab- 
hangt. 

Wenn wir durch Integration iiber YW’ (9.3) nach ¥" auflésen und 
beriicksichtigen, daB Y4y, nicht normiert ist 


* (Pay, Pay)* Poy 
yo — 1 [Yon d= | (tay Ww’) oe Ww dW’, (9.7) 
(22)8 : (2z)§ (Pay, Pay)! 
erscheint PEKARs Zustand als ,,Wellenpaket‘‘ aus Zustanden Yq. Die 
klassische Vorstellung eines bewegten lokalisierten Polarons kann nur 
eine brauchbare Naherung sein, solange das Wellenpaket nicht merklich 

auseinander flieBt. 
Fiir die weitere Rechnung ist es zweckmabig, (9.1) so umzuformen, 
daB es die Gestalt 
Voy = Sx (9.8) 
annimmt, wo S der in (3.2) definierte Operator ist. 
Voy = S+ f ef B(R-2) e-ipR o- i Dede be(R—xt) aR. (9.9) 


Die Wirkung der beiden letzten Exponentialfaktoren auf ¥{° kennen 
wir aus 7. Mit %—vr-—> und der in (7.3) eingefiihrten Bezeichnung 


fiir die-Gitterdeformation mit dem Mittelpunkt R entsteht 
4 = Se PX yp(R) OMAR. (9.10) 


Nun muBten wir mit (9.8), (9.10) den Erwartungswert von H ausrechnen. 


Zuvor wollen wir unser Ergebnis jedoch noch unter einem anderen Ge- 
sichtswinkel betrachten. 


10. Allgemeine Form unseres V ariationsansatzes. 


Die Vorgabe eines Variationsansatzes fiir y (9.8) kann dadurch er- 
folgen, da® man von einem Operator U ausgeht, der eine Funktion aller 
be und 0 ist und ihn auf den Vakuumzustand ®, wirken laBt. U enthalt 
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noch Funktionen und Parameter, die variiert werden kénnen, so 
daB U®, eine ganze Klasse von Zustinden umfabt. Zu einem brauch- 
baren Ansatz gehért nicht nur, daB die Mannigfaltigkeit der beriick- 
sichtigten Zustande gut ausgewahlt ist, sondern es muB auch moglich 
sein, den Erwartungswert des HAMILTON-Operators mit ertraglichem 
Aufwand auszurechnen. Bei der Behandlung des BLocu-NorpsIECK- 
Problems ist uns ein unitarer Operator begegnet (5.11), der als Ver- 
schiebungsoperator fiir die b; und }, wirkt, so da8 der Erwartungswert 
(U®,, HU®,) =(®,, U! HU®,) leicht ermittelt werden kann, Mit 
U=exp {/(f) by — /*(E) &} gilt [Formel (12) von LLP] 


U4, U=b + /*(8); U1b,U=h)h+f7(f). (10.1) 
Hierauf beruhen die Methoden von Tomonaca und LLP. Mit (3.7) folgt 
eZ (HOE HOE) GeO By, (10.2) 


Die allgemeinste Form unseres Ansatzes besteht nun darin, daB wir fiir 
die Operatorfunktion U eine Integraldarstellung 


f G(R) cZFENE ARD, (10.3) 


verwenden. Durch diesen Kunstgriff ist die Berechnung des Erwartungs- 
wertes von H, was die Operatoren betrifft, fast ebenso wie einfach bei LLP. 
Es kommen allerdings sechs 3i-Integrationen hinzu und man muB bei der 
Wahl von F und G darauf achten, daB fiinf davon und die f-Integrationen 
analytisch ausfiihrbar sind. Die letzte Integration kann dann numerisch 
gemacht werden. (10.3) umfaBt natiirlich eine viel gréBere Mannig- 
faltigkeit von Zustanden als (10.2). 

In 9. haben wir durch Anwendung eines Projektionsoperators auf 
die PEKARsche Lésung einen Spezialfall von (10.3) erhalten 


G(R) =y(R)e™*; F(R) =—in PM et (10.4) 


Eventuelle Bedenken gegen die mathematisch nicht strengen Methoden 
in 9. haben keine Folgen fiir die Brauchbarkeit des Ansatzes (10.4), denn 
ein Variationsansatz bedarf keiner Begriindung. Seine Giite zeigt sich 
bei der Berechnung der Energie. Wir werden in 12a ein Beispiel zu 
(10.4) bringen. 

In (10.4) ist 9 die Fourrer-Transformierte (5.3) von |p|’. Diese 
Beziehung wird jedoch bei der Berechnung des Erwartungswertes von H 
nicht benutzt. Sie ist lediglich von Bedeutung, wenn wir den Zusammen- 
hang mit der PEkKaRschen Lésung untersuchen. Ahnlich wie in 8. be- 
trachten wir daher kiinftig 9 und y als unabhangig voneinander und 
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fiigen auch wieder einen Variationsparameter 7 ei. 


Voy = 3%; g= Je" PUR) eae | 


ey pt & -—intR | 
g —tx) bee 
D®) — ¢ pF k 0° 


(10.5) 
Der Erfolg ist wiederum, daB bei geeigneter Wahl von 9, p, 7 als Funk- 
tion von g? ein Ansatz resultiert, der nirgends schlechter ist als der von 
LLP und iiberall besser als der von PEKAR. Im Gegensatz zu 8. haben 
wir jetzt aber Eigenzustande von YW. Einige Rechnungen mit (10.5) 
bringen wir in 12b. 

Nun wollen wir noch die Verbindung zu 4. herstellen. Dazu ent- 
wickeln wir (10.5) in der Form (3.9) und beschranken uns auf W’= 0 
und auf g, 


galt, te) ~ SE? 2 (ty + b)I- (10.6) 
1 2 
Dabei ist 
es —_ 1 ' —iRR 7S Y 
p(k) = ane | y(R)e am. (10.7) 


Speziell werden wir in 12b 9 =(1+?)1, p~e-** wihlen. Dann ist 
p(K) ~ [1 +o? K?)|~? und wir kommen mit (3.5) auf den in 4. angegebenen 
Ausdruck zurtick (7 = 1). 


11, Berechnung des Erwartungswertes von H. 


Der Zustand (10.5) ist bereits so geschrieben, daB wir leicht auf H 
(3.5) tbergehen kénnen!® 
(Yay, H Pay) (x, Hy) 


E (%") = = 
Me | (Voy, Pon’) (x. x) 


(41.1) 

Zuerst rechnen wir (DP), OP") aus. Dazu ersetzen wir die Summe im 

Exponenten durch ein Produkt von Exponentialfunktionen, schreiben 

fiir diese die Reihenentwicklung und finden nach einigen Umformungen 
. wr Ve 
mit 8 = R’— *R 

(Dr®) Din") = el(nZ)- L(y Z) = g? 3 0? (k) Po int8 dt (14 2) 

’ , / —_ 2 k2 — 

(OB), HOP") fiihren wir auf (14.2) zuriick, indem wir mit der Ver- 
tauschungsrelation 

§ 


[op eX PFE) FE ) eA ALR) (11.3) 


die by nach rechts durchziehen und mit den df dasselbe machen, nach- 
dem wir sie als 6, auf die andere Seite des Produkts gebracht haben. 


‘© Die Klammern beziehen sich im mittleren Term auch auf die Variablen des 
Elektrons, im rechten Term jedoch nicht mehr. Ein Faktor V wurde gekiirzt. 
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Der vom Gitter herriihrende Anteil ist dann 


(Da), >} be by Pink") ss, L(y Z) el (nZ) | (11.4) 
Beim Riicksto8term nehmen wir die Zerlegung (4.6) vor und finden 
(2B) 2 e! (2), (11.5) 
g* 2 re 0? (Rk) —intR q 
£, [Ge — ate) 20) cinta apetna, (11.6) 
Z If ry OP(R) QP (R’) ,—ing(tse) 1(Z) 1(nz) | 42 (Z) |? 
Ge a 78(t+ r o(3z). ol (7 Lt \)) ) |? = 
4a J, (Ef) Tire rh ASA aT ai TA ae 


Der Wechselwirkungsterm bringt wie (11.7) einen negativen Beitrag 
gf elk) (intr pint® (1 
os f ate ne beh Mt) dE gi (02). (11.8) 


Uber die nun auszufiihrenden f-Integrationen ist allgemein nichts zu 
sagen. Bei den ‘t-Integrationen erhalten wir Integrale des Typs (wp ist 
reell) 

T=S p(B’) y(R") &23 (3) aw aK”. (11.9) 


Fiir die drei Funktionen im Integranden setzen wir Integraldarstellun- 
gen ein* [vgl. (40.7) | 

h(3) =f &8Fh(g) de. (11.10) 
Die Integrale iiber 3’ und *%t” ergeben 6-Funktionen, so daB auch die 
RX’ und die &’’-Integrationen sofort ausfiihrbar sind. Es bleibt 


J = (22)? f &®3h(3) f e #3 [G(K)PdRdB. (4.114) 


Da das letzte Integral nur von 3 abhangt und in 12. nur /(Z) vorkommt, 
kann die Integration tiber die Winkel noch allgemein ausgefiihrt werden 
Jim g2zct f BOP LZ" bi2) fo=t#2 [p(K)24@dZ. (14.42) 
Diese Formel erlaubt die Berechnung der Normierung sowie der Gitter- 
und Elektronenanteile der Energie. Die Wechselwirkungsenergie gehért 
nicht zum Typ (11.9); man kann sie entweder direkt ausrechnen oder 
auch Integraldarstellungen (11.10) verwenden. Es geniigt, sich auf den 
ersten Term in (11.8) zu beschranken, denn der zweite bringt noch 
einmal dasselbe. 
12. Numerische Durchfiihrung. 
a) Wir haben in 5. gesehen, daB die einfache Gauss-Funktion beinahe 
ebenso gut ist wie der beste PEKARsche Ansatz. Daher wollen wir mit 
ihr die Giite des in 9. geschilderten Verfahrens erproben. y und ¢@ sind 


* Diese Methode zur Auswertung von J verdanke ich Herrn Dr. LEIBFRIED. 
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in (5.13) angegeben. Mit y=1Z/6, §=g?/P erhalten wir fiir die Gitter- 
energie 

ey is Se 
a - fe vig 2 2 ey dy (12.1) 


0 


Eg= 


@(y) bedeutet hier das Fehlerintegral. Der Erwartungswert der kine- 
tischen Energie des Elektrons besteht, wie in 4. naher erlautert, aus einem 
positiven und einem negativen Anteil 


= 22%) ee .> 

Ee = : s fe y e 27° PY er ae 
Nuf \\a 

0 


[ 


(2) May. (12.2) 


Die Wechselwirkungsenergie betragt 
" PY i. MPS o % ? 
Rae gr a ag , 42. 
Ey 4Efe ye Fy (ae) ) (12.3) 


0 


N ist bis auf einen Faktor das Normierungsintegral 


© oy)  ¥ 
= | e vy e€ tw ytdy, (12.4) 


. 


0 


Im Grenzfall kleiner Kopplung gilt 
= ~ {8 —1 —-*} gt. (12.5) 


Mit dem giinstigsten Wert von = |7°+-1/2 ergibt dies E = — 0,77? 
Da PEKARs Kurve wie g* einmiindet, ist die Verbesserung betrachtlich. 
E = — g? wird nicht erreicht, weil o von der fiir kleine Kopplung giinstig- 
sten Funktion 1/(1 + ®) verschieden ist. Numerische Rechnungen sind 
nur fiir 71 ausgefiihrt worden. 

FE ist eine Funktion von f und g*. Eigentlich miiBte man zu jedem g?® 
das giinstigste 8 bestimmen. Es ist aber einfacher, E als Funktion von 
&=g?/8 und g* zu betrachten und zu jedem g® das beste & zu suchen. 
Ee-+Ey hangt tiberhaupt nur von & ab und zwar beinahe linear 


7 ara ees “1,06&; ge<1 
« - Tv a iE; 
Eg + Ew =A(é) = ; ‘ (12.6) 
— ——€; ge > 141. 
Ja 


A(&) und das durch Ey, = B(&) g~4 definierte B(é) werden durch nume- 
rische Integration bestimmt. Wenn wir E iiber g? auftragen, gibt es 
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zu jedem & eine Hyperbel 
E (g®, €) = A(é) + 2 : (12.7) 


Die Einhiillende ist unser gesuchtes E (g?), Fig. 2 zeigt, daB wir eine 
merkliche Verbesserung der PEKARschen Kurve erreicht haben. Der 
Schnitt mit E——g* liegt jetzt bei g2=7,2 gegeniiber dem fritheren 
Wert 9,212". 
b) In 12a wollten wir feststellen, 0 
welche Verbesserungen die in 9. vor- 
genommene Modifikation des PEKAR- L 
schen Ansatzes fiir E(g*) bringt. Un- 
befriedigend bleibt vor allem das Ver- [ 
halten fiir kleine Kopplung. Man _ ,| 


kénnte daran denken, fiir 9 = 1/(1 + ?) 
und fiir py die mit diesem o aus (5.3) £& 
zu berechnende Funktion einzusetzen. | 
Dann kommt man aber erstens auf 
schwierige Integrale und zweitens 
divergiert in der zugehérigen Lésung 
vom Produkttyp (5.1) die kinetische 
Energie des Elektrons. Daher gehen 
wir zu (10.5) tiber, wahlen 0 und wy 
unabhangig voneinander und setzen *% 
o(k) =[1 + (BR/2)*)-"; etc Ligeti Hein 
poe 


(12.8) in 12a. und 12 b. Die Werte von Ler und Pines® 
sind mit * bezeichnet. 


o und y sind so beschaffen, daB fiir die Werte n=1, B=2, 7=0 der 
Variationsparameter der Ansatz von LLP resultiert (%%’=0) und fiir 
die Werte n=2, «=f, 7=1 ein Ansatz des in 9. diskutierten Typs 
vorliegt. Da wir uns vorwiegend fiir mittlere Kopplung interessieren, 
haben wir zunadchst nur 1 gerechnet. 


Es sei w=2nZ/B, u=2g7/B, A=2an/B 


mht Lee, ayy mteli 04 Po 
[ w? 2 12: 
etree Zp 
m,(w) = 4w (1 — e-”) m,(w) = = 


17 Der Unterschied zwischen dem Spezialfall (9.1) unseres Ansatzes (10.3) und 
dem PeKarschen Produktansatz betragt etwa 3/2hw (Fig. 2). PEKaRs adiabatische 
Naherung? liegt ebenfalls um 3/2hq tiefer als der Produktansatz. Eine einfache 
Beziehung zwischen beiden Theorien ist nicht ersichtlich, es lohnt sich aber, den 
Zusammenhang genauer zu untersuchen, weil dabei zugleich endgiiltig geklart 
werden kann, worauf PEKaRs groBer Wert fiir die Masse des Polarons zuriickzu- 
fiihren ist (mp=17m,, vgl. FuBnote auf S.416 in PEKARs Bericht). 
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Dann wird 
* 1 r Lm, (w)—w/A ws w 1 ag ”"0) “9 ‘ 
com tnfennt( +P atlmiadde, 240 
0 
7 et A Sf gum(w)—wlaf A) @ 4 ay2\|_* ow )| dew, 12.44 
En=ya |? (at t+@)[|5 6% me) (12.14) 
0 
1 um, (w)—wA . 

Ey =— ow fe sedi | 
pes (12.12) 
fmt), 20+4 30 [2 a3 | w)} dw. 

“tenga te © Las ee 
N ist bis auf einen Faktor das Normierungsintegral 
N= cae We + w*) dw. (12.13) 
? = i 
0 
Die Energie E =E,+ Ey, + Ey hat wieder die Form 
E (g2; 4, uw) =A(A, pw) + ag ; (12.14) 


Die drei Variationsparameter «, $,7 kommen also nur in der Kombi- 
nation A, vor. Hieraus entsteht die gesuchte Funktion £ (g*), wenn 
man fiir jedes g? dasjenige Wertepaar A, « einsetzt, das E zum Minimum 
macht. Da die Berechnung von A und B fiir ein Wertepaar bereits 
mehrere numerische Integrationen erfordert, beschranken wir uns zu- 
nachst auf A=1. A und B wurden fiir ~=5, 10, 20, 30 berechnet und 
konnten dann fiir Zwischenwerte gut interpoliert werden, weil A() 
beinahe linear ist und B(u) sich beinahe wie uv’ verhalt. Dann wurden 
die Hyperbeln (12.14) gezeichnet. Ihre Einhiillende ist eine Naherung 
fiir das zu unserem Ansatz gehérige E (g*), die sicher zu hoch liegt, weil 
wir die Werte fiir » und A fest gewahlt haben, statt sie zu variieren. Die 
Kurve ist in Fig. 2 dargestellt. 

Da A= 0 auf ®,;p zuriickfiihrt, erwarten wir fiir kleine g? kleine Werte 
von A. (12.14) gibt bei einer Entwicklung nach Potenzen von 4 zunachst 


in nullter Ordnung 
3 


ef D - 2 1 
E (g** 4) = at oa (42.15) 
Das beste wv ist gleich g? und es wird, wie erwartet E = — g®. — Die erste 


Ordnung hangt nur von |A| ab!® und muB einen positiven Koeffizienten 
haben, weil nach 4. jede Verbesserung von dem dort definierten H’ 
bewirkt wird und dieser Anteil der kinetischen Energie des Elektrons 


#8 Wenn wir auch negative 7, 4 betrachten, muB im Ergebnis der f-Integrationen 
in 11. || statt 7 stehen. «, B sind stets positiv. 
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erst bei 2? beginnt. Fiir sehr kleine 4 gibt also (12.14) stets eine héhere 
Energie als mit A=0. 

Wenn wir zu jeder Kopplung das beste Wertepaar A, als Punkt 
in der A, u-Ebene markieren, entsteht eine Kurve, deren Punkten Werte 
von g* zugeordnet sind. Nach obiger Uberlegung beginnt sie bei g=0 
im Ursprung und hat zuniachst die Richtung der v-Achse. Bei einem noch 
nicht bekannten g?, das nach Fig. 2 jedoch unter 5,5 liegen muB, macht 
sie einen Sprung, weil sehr kleine A-Werte, d.h. kleine Korrelationen, bei 
keiner Kopplung die giinstigsten sind. 


Schlupbemerkung. 

Das am Ende von 12. angedeutete Verhalten kénnte ein grundsatz- 
licher Mangel des Ansatzes (10.3) sein. In diesem Falle miiBte man ihn 
so erweitern, daB die Abweichungen des exakten Grundzustandes von 
®,;p besser approximiert werden. Einmal bietet sich die Méglichkeit, 
den Ansatz (9) der Arbeit von LEE und PinEs® durch Ubergang zu einer 
Integraldarstellung ahnlich zu verallgemeinern, wie wir es mit dem Ansatz 
von LLP getan haben. Andererseits kann man Anregungen fiir einen 
Variationsansatz bei mittlerer Kopplung erhalten, wenn man die Eigen- 
zustande erster Naherung derjenigen Stérungsrechnungen, die in nullter 
Naherung E = — g? liefern (4. sowie LLP), in der Form (3.9) darstellt. Wir 
folgen damit Methoden von FROHLICH, PELZER, ZIENAU!® und GURARI?. 

Die Frage nach der Tragheit des Polarons kann beantwortet werden, 
sobald F (8?) fiir kleine YW’? bekannt ist!*. In unserer Naherung macht 
die Berechnung dieser Funktion keine besonderen Schwierigkeiten, die 
numerische Auswertung ist jedoch noch nicht abgeschlossen. @, das 
jetzt auch von der Richtung von f abhangen kann, ist so zu wahlen, daB 
E (%8'?) méglichst tief liegt, die Masse also méglichst groB wird. — Es ist 
keinesfalls sicher, daB ein Verfahren, daB eine tiefere Energie fiir den 
Grundzustand liefert, nun auch einen genaueren Wert fiir die Masse 
ergibt. Da die Masse so empfindlich von der Art der Naherung abhangt, 
ist ihre Berechnung eigentlich erst dann von Interesse, wenn zugleich 
die Abweichung vom exakten Wert abgeschatzt werden kann (Stérungs- 
theorie von LLP, s. auch, S. 353). Auf die Kritik! an PEKARs Methode 
werden wir zuriickkommen, wenn die in !7 genannte Untersuchung ab- 
geschlossen ist. 

Zusammenfassung. 

Es wird ein Variationsansatz fiir die Berechnung des Grundzustandes 
des Polarons angegeben”, der folgende Eigenschaften hat: 

4. Er ist eine Verallgemeinerung der Ansitze von Gurari und LLP 
durch Beriicksichtigung von Korrelationen zwischen den Schallquanten. 

9 Préuzicu, H. Pewzer u. S. Zrenau: Phil. Mag. 41, 221 (1950). 

20 Nach einem Hinweis von FROHLICH? (S. 355) kénnte das Verfahren auch fiir 

die Supraleitung von Interesse sein. 


914 G. Héurer: Wechselwirkung eines Teilchens mit einem skalaren Feld. I. 


2. Er ist zugleich eine Verallgemeinerung des Produktansatzes von 
Pekar, die stets Eigenzustande von ¥ liefert. Damit ist die Beziehung 
zwischen diesen beiden Ansatzen geklart. 

3. Die Energie liegt fiir alle Kopplungsstarken nicht hoher als bei 
den in 1. und 2. genannten Autoren. Vorlaufige Rechnungen, bei denen 
die verfiigbaren Variationsparameter nur teilweise ausgenutzt wurden, 
ergaben das in Fig. 2 dargestellte Resultat. Unsere Eigenzustande von YW 
liegen also fiir alle Kopplungsstarken tiefer als die lokalisierten Produkt- 
zustande PEKARs. 

4. Es ist wiinschenswert, das Verfahren auf andere Typen der Kor- 
relation zu erweitern, damit das noch nicht zufriedenstellende Verhalten 
bei kleiner Kopplung verbessert wird. 


Ich danke Herrn Professor BECKER fiir sein freundliches Interesse 
an der Arbeit, Herrn Dr. H. LEHMANN fiir zahlreiche Diskussionen und 
Herrn A. MULLENSIEFEN fiir Unterstiitzung bei den Rechnungen. 


Anmerkung bet der Korrektur: Fir starke Kopplung kénnen die 
Integrale in 12a ohne numerische Rechnung ausgewertet werden. Man 
erhalt eine Entwicklung nach fallenden Potenzen von g*: 

tear kid Were 7 


- 
5 
32 = a 


Die Differenz zur Energie des Grundzustandes in PEKARs adiabatischer 
Naherung !” verschwindet also wie g~‘. 


Gottingen, Institut fiir Theoretische Physik der Universitat. 


+ 


Zeitschrift fiir Physik, Bd. 140, S. 215—220 (1955). 


Zerfall und Erzeugung von A-Teilchen. 
Von 
GERNOT EDER. 
Mit 1 Figur im Text. 


(Eingegangen am 19. November 1954.) 


Auf Grund einfacher feldtheoretischer Wechselwirkungsansitze wird das Ver- 

haltnis der Lebensdauern von freien und gebundenen A-Teilchen berechnet; die 

Konsequenzen dieser Resultate fiir verschiedene A-Erzeugungsmechanismen werden 
diskutiert. 


CroK und Mitarbeiter! haben auf Grund der vorliegenden Messungen 
gezeigt, daB im Kern gebundene A-Teilchen (A) eher direkt in ein 
Nukleon (N) iibergehen, als entsprechend dem Schema A-+N +2 zu 
zerfallen. Die beiden Prozesse sollen hier feldtheoretisch miteinander 
verglichen werden. 

Das A-Feld kann man durch einen Spinor oder einen Pseudospinor 
beschreiben. Wird es durch einen Spinor beschrieben, dann nimmt die 
Wechselwirkung zwischen den Nukleonen und A-Teilchen die Form 


Hy =18 (Py Ys T; Wa) ®; (1) 


an. Die reziproke Lebensdauer eines freien, ruhenden A-Teilchens ergibt 
sich daraus zu 
w, = (g?/420) [(M? + m? — 2)? (2M)-2 — m?}! x | 0) 
% (2M)-2 [(M — m)? — p2]t, j 
wo , m und M die Pion-, Nukleon- und A-Masse bedeuten. Numerisch 
bedeutet dies: w(A° +p +27) =3,6 - 10 sec? (g?/4); fiir den Zerfall 
A°->n-+2° ergibt sich der doppelte Wert. Mit t=3-107 sec fir 
A°+p+a” folgt 
(g?/42) = 0,93 -107"". (3) 


Der Zerfall eines im Atomkern gebundenen A-Teilchens kann auch 
ohne z-Emission vor sich gehen und zwar dadurch, daB der UberschuB 
an Energie und Impuls an ein Nukleon abgegeben wird. CHESTON und 
PRiMAKoFF? haben den nichtmesonischen Zerfall von A-Teilchen nach 
der pseudovektoriellen Theorie in nichtrelativistischer Naherung be- 
rechnet. Hier wird der kovariante Formalismus zur Berechnung der 


1 eyo, Pp, M. Danysz u. J. GrERULA: Nuovo Cim. I1, 436 (1954). 
th=c=1. 

2 CHESTON, W., u. H. PRIMAKOFF: Phys. Rev. 92, 1537 (1953). 
Zeitschrift fiir Physik. Bd. 140. 15 
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Ubergangswahrscheinlichkeit verwendet. Nehmen wir fiir die Nukleon— 
Nukleon-Wechselwirkung einen Term der Form 


H' = if (Py ys t; Py) Di (4) 


an, so erhalten wir in zweiter stérungstheoretischer Naherung entspre- 
chend dem Ubergang Fig.1 fiir die Ubergangswahrscheinlichkeit pro 
Mp) Mg) = Np) Ag) Nukleon im Kern w/A nach Mittelung tiber Spin 
und isotopen Spin (wir nehmen das Anfangsnukleon 
und das A-Teilchen als ruhend an): 


ath) alk) i, 
Ww 1 
a — >| eg (T ° Tx)? x 
Np) Aig) Njfp) Nolq) | 1 ae (5) 
Fig.1. 2-freier A-Zerfall _, (M — m)? |) (M + m)?/4 — m* 
bee STE sk lag Riel ha: ri ye + m(M — m)}?(M +m) ~ 
N—N-StoB. - 


Setzen wir fiir das Kernvolumen V = (42/3) A und vergleichen wir die 
Wahrscheinlichkeiten fiir den Zerfall ohne bzw. mit a-Emission, so 
ergibt sich aus (2) und (5): 
re 
=a 3 F (t4° ty) X sucfles 
: (\Z — m)? M? 3 (M + m)? — 4m? 


(M+ m) [u? + m(M— m)}? [((M — m)? ane pe) (M2+ m? — y2)2/ M2 — 4m? > 


Numerisch bedeutet dies: 
w/w, = 1,634 (f/42) (t4* Ty)®. (6) 
Fiir (v..° ty) =1 wird daher mit (f*/42) = 10: 
w/w,= 16,31. (7) 


Ein im Atomkern gebundenes :1-Teilchen wird also eher ohne Emission 
eines Pions zerfallen als nach dem Schema A—»+N +a. 

Es sollen nun einige Konsequenzen dieser Ansatze zum A-Teilchen- 
zerfall fiir die Erzeugungsprozesse diskutiert werden. Der zum zweiten 
Zerfallsschema inverse ProzeB (vgl. Fig.1) wird wegen der kleinen 
Konstante g einen sehr kleinen Wirkungsquerschnitt haben. Um ein 
Zahlenbeispiel anzugeben, nehmen wir fiir das einfallende Nukleon 
eine kinetische Energie von T= 2,2 GeV an; fiir die Erzeugung in 
der Einfallrichtung ergibt sich mit (3) numerisch: do/dQ(#=0) = 
5,2- 10°! mb (ty, Ty)? im Lab.system; dabei wurde (/2/427) = 10 gesetzt. 
Die Schwellenenergie fiir diesen ProzeB betragt im Schwerpunktssystem 
Fcnw = 4 (M +m), dh. Exay =m +373 MeV im L.S. Obwohl also die 
Schwellenenergie sehr niedrig liegt, ist der Wirkungsquerschnitt fiir 
A-Teilchenerzeugung sehr klein, wenn er durch die Wechselwirkung H’(1) 
vermittelt wird. 
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Ein anderer Erzeugungsmechanismus fiir .1-Teilchen liegt im Falle 
der Paarerzeugung vor: Beschreiben wir die a — A-Wechselwirkung 
analog zu (4) durch einen Term 


Hs =?f' (Paysti pa) ®, i) 


so vermittelt (8) unter anderem den Ubergang eines (virtuellen) Pions 
in ein A-Teilchen und ein entsprechendes Antiteilchen (1). Der ein- 
fachste ProzeB, der eine 1-Paarerzeugung ermoglicht, ist die Nukleonen- 
paarvernichtung, deren Bedeutung allerdings vom Nachweis eines Anti- 
protons oder Antineutrons abhangt. Das Matrixelement fiir einen solchen 
Ubergang ergibt sich aus H’ = (4) + (8) zu 


<A (p')Alg’) |Se| N(6) N(q)> =i ff fax d*y A. (x—y) (P_gystiyp(2)) x 
x (Wy: V5 Ti Y-q (¥)) : 


Daraus folgt ein totaler Erzeugungsquerschnitt von 


(9) 


O4 = (ff'/4a)? (ty TH)? (* + 1) (* + 1 — 2M?/m?)? x 


1 39 
‘N > 


* (x + 1 — u?/2m?)-? (x — 1)7 22 


wo x~=(1—v?)~4 (v = Geschwindigkeit des einfallenden Nukleons im 
L.S.) und 4, die Compton-Wellenlange des Nukleons bedeutet. Die 
Schwellenenergie betragt fiir diesen ProzeBim L.S. : E,4,, =m + 779 MeV. 

Uber die GréBenordnung von f’ geben die Messungen der A°-Bin- 
dungsenergie in 3H! und in Be? Auskunft: Die A°-Bindungsenergie in 
3H ergab sich zu (0,9 +1,1) MeV im Vergleich zu einer Neutronbindungs- 
energie von —6 MeV. Nehmen wir eine kinetische Energie pro Nukleon 
in der GréBe der FERMI-Energie T = 14,7 MeV an, so wird die potentielle 
Energie des A°-Teilchens bzw. Neutrons: (— 16+ 1) MeV bzw. — 21 MeV. 
Wenn der zu (4) analoge Ansatz von (8) berechtigt ist, so folgt daraus, 
daB f und /’ in der gleichen GroéBenordnung liegen. Fiir (ty + ty) =1 
und /?=/’?=40za ergibt (9): 


G4 = 139 mb (x + 1) [(x +1 — 2M?/m*)/(x — 1)}§ (# +1 —p?/2m*)*. 


Etwa fiir x =2 (E =2m) erhalten wir 0, =19,2 mb. 

Im Normalfall wird sich ein Nukleon-Antinukleonpaar in ein Pion- 
paar umwandeln durch Vermittlung des Paarterms, welcher der Wechsel- 
wirkung (4) entspricht: 


Hy’ = (f°/2m) (py py) ®- (10) 


Das Matrixelement fiir die Vernichtung eines Nukleons mit dem Energie- 
impulsvektor # und eines Antinukleons (q) unter Aussendung zweier 


1 Bonetti, A. u. a.: Nuovo Cim. I1, 210, 230 (1954). 
2 Fry, W. F., u. G. R. WHITE: Nuovo Cim. 11, 551 (1954). 
15™ 


218 GERNOT EDER: 


Pionen (k und k’) lautet nach (10): 
Cox* (it) 2 (k’) |S,| N (6) N(q)> = — i (PP/2m) f d*x (Pg ¥p) De De (x) - 
Der totale Wirkungsquerschnitt ergibt sich daraus zu 
= (0/4) (f2/420)? [(% — 1) (% +1 — 2p2/m?)}* (x +1) 14%. (14) 


Das Verhiltnis der Wirkungsquerschnitte fiir die Erzeugung eines A- 
bzw. a-Paares betragt nach (9) und (11): 
a iis (x + 1)? x +1 —2M*/m? 
pied. ( ’] (ty tH)" (% — 1) («+ 1 —p2/2 m2)? »| x+1—2y2/m? * (12) 
Fir /’ =f, (ty * ty) =1 und x = 2 bedeutet dies numerisch: o,,/o, = 1,04. 
Es wiirde also ungefahr jedes zweite Nukleonenpaar bei seiner Vernich- 
tung ein A-Paar erzeugen. 

Betrachten wir schlieBlich den Paarterm, welcher der Wechsel- 
wirkung (8) entspricht: 


Hs = (f'?/2M) (yi ys) @. (13) 


(13) vermittelt unter anderem die Erzeugung eines A — A-Paares durch 
z-Paarvernichtung. Das Matrixelement niedrigster Ordnung lautet fiir 
diesen ProzeB: 


<A (p') A(q’) |Sy| 0+ (6) -e (Q)> = — i (f'2/2M) Jf dx (Py p_y) Bf Dy (x). 
Daraus folgt ein totaler Wirkungsquerschnitt : 
= (21/2) (f'2/4a0)? (x + 1 — 2M®'?)? (xe + 1) 72 (xe — 1) ~4 (m/M)22%. (44) 


Die Schwellenenergie liegt sehr hoch: Ey = 2.M2/u —u =u + 17,43 GeV. 
Mit /’? = 40a” wird aus (14): 


o = 49,1 mb (x + 1 — 3M2/p?)8 (x + 1) (x — 1)- 


Das bedeutet z. B. fiir x = 150 (E = 21,09 GeV): a= 3,33 mb. 

Die drei fiir A-Erzeugung angefiihrten Mechanismen sind entweder 
wegen der kleinen Kopplungskonstante oder wegen der hohen Schwellen- 
energie unwahrscheinlich. ONEDA! hat versucht, einen gréBeren Er- 
zeugungsquerschnitt zu erhalten, indem er eine hinreichend starke 
Wechselwirkung der t-Mesonen mit dem N- und A-Feld annimmt. Es 
fihrt aber auch diese Annahme zu Schwierigkeiten: Beschreiben wir 
das t-Meson durch eine pseudoskalare FeldgréBe ®7, so kénnen wir fiir 
die A—N—t-Wechselwirkung folgenden Ansatz machen: 


Hy = 1G (Pays ti Py) DF + iG (Wy ys T py) PF. (15) 


1 OnepA, S.: Progr. Theor. Phys. 6, 633 (1951). 
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(Da ein Zerfall des A-Teilchens in ein Nukleon und ein t-Meson aus 
Energiegriinden ausgeschlossen ist, so bleibt die Lebensdauer der freien 
A-Teilchen durch (15) unverandert.} Durch (15) wird bei Nukleon— 
Nukleon-St6Ben entsprechend Fig.1 die Erzeugung von A-Teilchen ver- 
mittelt. (Das ausgetauschte Meson ist nun ein t-Meson.) Die Schwellen- 
energie ist natiirlich die gleiche wie beim ersten ErzeugungsprozeB. Bei 
einer kinetischen Energie von 2,2 GeV (L.S.) des einfallenden Nukleons 
betragt dann der differentielle Wirkungsquerschnitt in der Einfalls- 
richtung in niedrigster Naherung: 


do/dQ (9 = 0) = 0,046 mb (GG"/42)? (ty, Ty,)?. (16) 


(16) gabe plausible Werte, wenn man (G*/427) in der GréBenordnung von 1 
annehmen kénnte. Da aber (15) zum z-freien A-Zerfall beitragt, so wird 
die von den t-Mesonen vermittelte Zerfallswahrscheinlichkeit der 
A-Teilchen im Kern um den Faktor 


(GG w+m (AM — m) \2 G\2/G 
Fr Paar Sayre eae) a2) 
groBer als (5). Fir G’/f0,1 muB (G?/42) © 10-8 sein, damit die A-Le- 
bensdauer in der GréBenordnung von 107" sec bleibt. Das bedeutet aber, 
daB der Beitrag der t-Mesonen zur A-Erzeugung sehr klein wird. 
Nimmt man an, daB sich die A-Feldgr6éBe wie ein Pseudospinor trans- 
formiert, so wird 


Ay = 8 (py T pa); (1’) 
und die Zerfallswahrscheinlichkeit fiir ein freies, ruahendes A°-Teilchen: 
w(A° > + 2°) (2 
= (g?/4z) [((M? + m? — ?)? (2M)-? — m?}} [(M + m)? —p?] (2M) ?; 


dies bedeutet, daB 
(g?/422) = 0,26: 10° (3’) 


wird. Fiir das Verhaltnis der Zerfallswahrscheinlichkeiten gebundener 
und freier A-Teilchen ergibt sich 


ae pee Se 

Wa 4 42 % 7 
M—m (M+ 3m) M? 13 ee +m)? —4m?* 
M+m - (M+ m)?— pi?) [u2+ m(M— m)}? Y (M? + m? — 40)/M? —4 me 


Numerisch bedeutet dies: 
w/w, = 0,305 (f?/42) (ta * Ty)” (6) 
(6) unterscheidet sich von (6’) um den Faktor 


(M — m) [(M + m)® — 2] 3 (M + 3m) [(M — m)? —j2)]}4 = 5,35. 
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Der schnellere Zerfall des gebundenen A-Teilchens nach der Spinor- 
theorie im Vergleich zur Pseudospinortheorie kénnte vielleicht bei einer 
groBeren Zerfallsstatistik zwischen beiden Theorien entscheiden. Der 
A-Erzeugungsquerschnitt beim StoB zweier Nukleonen (vgl. Fig. 1) wird 
in der Vorwirtsrichtung fiir eine kinetische Energie von 2,2 GeV 


do|dQ (8 = 0) = 0,79 - 10 mb (/?/42) (ty, Tx)? 
Der Wirkungsquerschnitt fiir die A-Erzeugung durch Vermittlung der 
t-Mesonen wird 
do/dQ (9 = 0) = 0,216 mb (GG'/42)? (ty, * Ty). (16’) 
Der Faktor (17) geht iiber in 
0,77 - 10!3 (G’/7)? (G*/4z), (17’) 
so daB auch nach der Pseudospinortheorie G sehr klein sein miBte: 


(G?/42t) ~ 1074. 
Fiir die pseudovektorielle Kopplung 


Hy =g (/(M + m)) (pyo-t, ps) V®, (1”") 
erhalten wir 
Ww, = (g?/4a0) 3? (2M)-? [4 — ue? (M + m)-?), (2”) 


a 


eat = (M? + m? — ?)? (2M)~?— m*. Daraus folgt mit w=! = 3 - 107 sec: 
(e2/42) = 1,11 - 1071, (3°) 


also ein Wert, der im wesentlichen mit (3) iibereinstimmt, wie auf Grund 
des Aquivalenztheorems zu erwarten ist, wahrend sich in der Naherung 
von CHESTON und PRIMAKOFF (g?/42) =0,31 - 10°" ergibt. 


Gottingen, Max-Planck-Institut fiir Physik. 


Zeitschrift fiir Physik, Bd. 140, S. 221—240 (1955). 


Zur Statistik der Elektronenlawinen. 


Von 
WERNER LEGLER. 
Mit 11 Figuren im Text. 


(Eingegangen am 19. November 1954.) 


Die statistischen Probleme der Elektronenlawinen und der raumladungsfreien 
TowNSEND-Entladung werden untersucht. Ausgehend von der statistischen Ver- 
teilung der LawinengrGBe, die auf einfache Weise abgeleitet wird, wird die Struktur 
der einzelnen Lawinen und die Verteilung der Nachlieferungselektronen (durch 
positive Ionen und durch Photoeffekt) berechnet. In der Umgebung von 4 =y- 
(e~4_ 1) —1 treten mit merklicher Wahrscheinlichkeit Lawinenketten auf. Die 
zeitliche Folge der Generationen, die Verteilung der Generationenzahl und die 
Struktur dieser Lawinenketten werden untersucht. Dabei werden besonders die 
durch ein einzelnes Primarelektron ausgelésten Ketten behandelt. Ist w> 1, so 
kann ein einzelnes Elektron mit der Wahrscheinlichkeit 1— 1/u die Ziindung der 
Entladung herbeifiihren. Das erste Anwachsen dieser ziindenden Entladung wird 
naher untersucht und daraus die statistische Verteilung der Aufbauzeit abgeleitet. 


I. Einleitung. 


Beim Aufbau einer Gasentladung aus einem einzelnen an der Kathode 
startenden Elektron ist es nicht zweckmaBig, nur mit den statistischen 
Mittelwerten der Gasverstarkung und Elektronennachlieferung zu rech- 
nen. DaB diese Rechnungen dann insbesondere in der Umgebung der 
statischen Ziindspannung der Gasentladung zu falschen Ergebnissen 
fiihren miissen, zeigt folgendes Beispiel: An eine Entladungsstrecke wird 
eine solche Spannung gelegt, daB sich ein Riickwirkungsfaktor «4 = 1,08 
ergibt. w=y- (e*4—1)=1,08 bedeutet, daB jede Elektronenlawine im 
Mittel 1,08 Elektronen an der Kathode auslést, die wieder zu neuen 
Lawinen anwachsen. Erfolgt der Aufbau der Entladung in Ionisierungs- 
spielen (TOWNSEND-Aufbau), so wiirde die primitive Rechnung mit den 
Mittelwerten zu folgenden Aussagen fiihren: Jedes von der Kathode 
startende Elektron fiihrt zur Ziindung der Entladung. In der zehnten 
Generation laufen im Mittel “® = 2 Lawinen, in der hundertsten Gene- 
ration j°* = 2000 Lawinen. Wie im folgenden gezeigt wird, ergibt die 
Beriicksichtigung der statistischen Verteilungen dagegen die Aussagen: 
Von den an der Kathode startenden Elektronen fiihren 7,4 % zur Ziindung 
der Entladung. Im Falle der Ziindung laufen in der zehnten Generation 
im Mittel 21, in der hundertsten Generation im Mittel 27500 Lawinen. 

Die Frage, mit welcher Wahrscheinlichkeit ein einzelnes Elektron 
die Ziindung einer Entladungsstrecke bewirkt, wurde von BRAUNBEK [1], 
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Hertz [2] und WiysMAN [3] behandelt. Dabei rechneten BRAUNBEK 
und Hertz noch mit der mittleren LawinengréBe e*4 und berticksich- 
tigten nur die Streuung im NachlieferungsprozeB. WIJSMAN stellte auch 
die nicht unerhebliche Streuung der LawinengréBe in Rechnung. Fir 
die Ziindwahrscheinlichkeit W der Entladung durch ein einzelnes von der 
Kathode startendes Elektron erhalt Hertz die Beziehung 1—W=e~"", 


wihrend W1yJsMAN die Darstellung W= 1 — a angibt. Fig.1 zeigt deut- 


lich den EinfluB der Lawinenstreuung, die 
hier zu kleineren Ziindwahrscheinlichkeiten 
fiihrt. 

In der vorliegenden Arbeit werden auf 
Grund der statistischen Verteilung der 
LawinengréBe und der Anzahl der Nach- 
Fig.1. Wahrscheinlichkeit W fir  lieferungselektronen einer Lawine folgende 
peccpena eee,  Worgange genauer statistisch beschrieben : 
als Funktion von « = y (e*4—1). Der Aufbau der Elektronenlawine, die Struk- 
Kurve @ (nach [2]) ohne, Kurve? tur von Lawinenketten und die Entwicklung 


(nach [3]) mit Berticksichtigung der 
Lawinenstreuung. der TOWNSEND-Entiadung. 


II. Verteilung der Lawinengr6Be. 


Startet in einer Gasentladungsstrecke von der Kathode ein Elektron, 
so wird bei geniigend hoher elektrischer Feldstarke durch fortgesetzte 
StoBionisation im Gas eine Lawine von insgesamt 7 Elektronen an der 
Anode ankommen. Die Zahl wird als ,,LawinengréBe'* bezeichnet. 
Unter gewissen Voraussetzungen kann die Wahrscheinlichkeit v(m, d) 
berechnet werden, daB ein Elektron eine Lawine von » Elektronen er- 
zeugt, wenn es die Strecke d in Feldrichtung zuriicklegt. Diese Voraus- 
setzungen sind: 

1. Das elektrische Feld ist eine reine Ortsfunktion und wird nicht 
durch die Raumladung der Lawine verzerrt. Dann ist namlich die 
lonisierungswahrscheinlichkeit eines Elektrons unabhangig von den 
anderen gleichzeitig laufenden Elektronen. Diese Forderung beschrankt 
die maximal zulaissige LawinengréBe (je nach den Entladungsbedin- 
gungen ist m,,. etwa 10® bis 107). 

2. Die Ionisierungswahrscheinlichkeit der Elektronen soll nur eine 
Funktion der am Orte der Elektronen herrschenden Feldstarke sein und 
nicht von dem seit der letzten Ionisation zuriickgelegten Weg abhangen. 

Um diese Bedingung zu erfiillen, muB das Feld relativ homogen sein, 
da die Geschwindigkeitsverteilung der Elektronen und damit ihre 
lonisierungswahrscheinlichkeit wesentlich von der Feldstarke abhangt, 
die kurz vorher auf die Elektronen eingewirkt hat. Nach Messungen von 
Morton |4) darf die Feldinhomogenitat im Zylinderfeld dann nicht mehr 
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vernachlassigt werden, wenn die Anderung der Feldstarke langs einer 
freien Weglange mehr als 1 bis 2% ausmacht. Dies kann bei Zaihlrohren 
der Fal sein. 

Da ein Elektron nach einer Ionisation im allgemeinen nur noch eine 
Energie besitzt, die klein gegen die Ionisationsenergie e- U, ist, so muB 
es erst etwa einen Weg s = U;/E in Feldrichtung zuriicklegen, um erneut 
ionisieren zu kénnen. Soll es fiir die Rechnung berechtigt sein, die 
Ionisierungswahrscheinlichkeit als unabhangig von dem seit der letzten 
Ionisation zuriickgelegten Weg anzusehen, so muB jedenfalls s klein 
gegeniiber der mittleren Ionisationsweg- 1000 
lange 1/x sein. Es muB also E/x>> JU, sein. 
In Fig.2 ist E/x als Funktion von E/p 
fiir Luft und Argon dargestellt. In Luft 
wird z.B. E/p kleiner als 100 V/cm: Torr 
sein miissen, damit E/« groB gegen U, wird. 
Das bedeutet aber, da8 fiir den Aufbau bei 
kleinen #-d-Werten (f-d<10 cm: Torr) 
diese Rechnungen nur mit Vorsicht an- 
wendbar sind, da fiir kleinere p - d-Werte 
die zur Ziindung notwendigen E/p-Werte 
sehr rasch ansteigen. 

“ ‘i Fig. 2. Mittlere Energieaufnahme E/a 

Unter diesen \ oraussetzungen (4 und eines Elektrons zwischen zwei Ionisatio- 
2) erzeugt ein Elektron auf dem Weg- 22 inAbhangigkeit von £/p fir Luft [0] 

a und Argon [6]. JIonisierungsspannungen 
differential dx in Feldrichtung mit der zum. Vergleich. 
Wahrscheinlichkeit «-dx ein weiteres 
Elektron durch StoBionisation und « ist nur vom Orte abhangig. 
Die LawinengréBe m unterliegt dann folgender Verteilung: 


Argon \luff 


Sonsrerungssponnungen: 
—hAr —N, 


0 20 50 100 200 nA 


a at — 
P cm Torr 


E xakte Verteilung Ndaherung fiir ny > 1 
li 4-\s—1 tay Thea 
} = 4 d) = - e€ 0 
v(n, ad) =i (1 el v(n, d) . 
Mittelwert: 
‘ (1) 
ni =n, = exp fs a(x)-dx —— 
0 
Streuung: 
o= n° V1— 1/No C= MN. 


Kurze Ableitung der V erteilung. 


u(n, x) ist die Wahrscheinlichkeit, daB ein Elektron beim Durchlaufen der 
Strecke x in Feldrichtung eine Lawine von insgesamt » Elektronen erzeugt. Be- 
rechnung von uv(n, x +d): Fir die Erzeugung einer Lawine von 1 Elektronen auf 
der Strecke x-+dx gibt es genau zwei einander ausschlieBende Moglichkeiten, 
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Einmal durchlauft das Elektron mit der Wahrscheinlichkeit (1 — « - d*) die Strecke 
dx obne zu ionisieren und fiihrt dann mit der Wahrscheinlichkeit v(m, x) auf der 
Strecke x zu einer Lawine von » Elektronen. Im anderen Falle erzeugt das Elektron 
zunichst mit der Wahrscheinlichkeit «- dx in dx ein zweites Elektron. Auf der 
anschlieBenden Strecke x bildet das eine Elektron eine Lawine von k Elektronen, 
das andere eine von n —k Elektronen. Die Wahrscheinlichkeit, daB beide Elek- 
tronen auf der Strecke x zwei Lawinen von insgesamt  Elektronen hervorbringen, 

n—1 
ist also >) v(k, x): u(n—k, x). Somit wird 

k=1 is 

u(n,* + dx) = (1—adx)-v(n, x) +adx- > u(k, x)-u(n — k, 2) 


oder 


d n—1 
—v(n, x)= —a-u(n,x) +a d v(hk, x) - v(m — k, x). (2) 
ax k=1 


Die Lésung dieses Differentialgleichungssystems gelingt durch Einfiihrung der 


erzeugenden Funktion 
V(E, zF <= Du (m, 2) -%. 
nn 


Fiir V(t, x) erhalt man nach (2) die Differentialgleichung 


eV hae 
—— =a V(¥ — 4)- 
rs 
Mit den Anfangsbedingungen v(1, 0) =1, v(m, 0) =O fiir » +1, also V(#, 0) =¢ er- 
halt man die Lésung 
~ 
: 1 f as 1 \"-1 
VG, 2) = —-——— , “ea [1 — Te 
My 1—(1— 1/m)-t Los my Ng 
n=1 
. 
Dabei wurde zur Abkiirzung exp J x(x) -dx—=n, gesetzt. Durch Koeffizienten- 
0 


vergleich erhalt man sofort die oben angegebene Verteilung (1). 

Fiir die weiteren Rechnungen wird die Naherungsdarstellung (1) fiir 
Ny-> 1 benutzt, die im folgenden als ,,ideale’‘ Verteilung bezeichnet wird. 

Eine ahnliche Verteilung findet man bei den Ausgangsimpulsen von 
Sekundarelektronen-Vervielfachern. Der StoBionisation im Gas ent- 
spricht hier die Sekundirelektronenauslésung an den Dynoden des Ver- 
vielfachers | 77}. 

In manchen Fallen ist es zweckmaBig, die Streuung der Lawinen- 
groBe als Streuung eines fiktiven Ionisationskoeffizienten «* aufzufassen, 
der durch n= e*"“ definiert ist. Daneben sei y= e**. (Im homogenen 
Feld ist % =a.) Sei AW die Wahrscheinlichkeit, mit der «* in ein Inter- 
vall Aa* fallt, so ist mit An=d- e*4-Aa* 

AW = > re": An = d-Aa*- exp {(a* — a) d — ee 4, 
Fir die Wahrscheinlichkeitsdichte w= AW/d -Aa* erhalt man (s. auch 
Fig. 3) 
w = exp {(a* — a9) d — of" od}, 3) 
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Die Wahrscheinlichkeit, daB® eine Lawine mit ihrem «*d den Wert yd 
© 

um mehr als A iiberschreitet, ist durch das Integral f w-d(a*d) ge- 
: x : X d-+-A 

geben. So findet man z.B. fiir Ad =2 bzw. 3 die Wahrscheinlichkeiten 


0,6 - 107 bzw. 2- 1079. 
Der fiktive Ionisationskoeffizient «* darf aber nicht so interpretiert 

werden, daB sich im ebenen Feld eine Lawine der GréBe n—e**4 nach 

dem Gesetz n(x) =e*"* aus einem Elek- a 

tron aufbaut. Diese Frage des Aufbaus 


. : Q3 
behandelt der nachste Abschnitt. 4 
| ge 
w 
Die Struktur der Lawinen. Qtr 
Die bekannte Formel fiir die Trager- a 
. . . “a my | y + * a 
vermehrung in einer Elektronenlawine (c¢*-cy)d —e 
Fig. 3. Verteilung des fiktiven Ionisations- 
dn=—n-a-dx n = e** koeffizienten «* einer Lawine. 


ist fiir kleine » sicher nicht exakt giiltig, sondern stellt nur eine Mittel- 
wertsaussage dar. Aus dem vorhergehenden ist bekannt, daB die La- 
winengr6Be einer starken Streuung um den Mittelwert 7) = e*@ unter- 
worfen ist. Welches Aufbaugesetz gilt fiir eine Lawine, die aus insgesamt 
n Elektronen bestehen wird? Die genaue Antwort auf diese Frage kann 
natiirlich nur in Form einer Wahrscheinlichkeitsverteilung gegeben 
werden. 

Es sei W,,(m,, x) die Wahrscheinlichkeit, daB eine Lawine, die aus 
insgesamt » Elektronen bestehen wird, an der Stelle x schon 1, Elek- 
tronen umfaBt. Bei der Ableitung der Verteilung wird vorausgesetzt, 
daB m groB gegen n, ist. 

Bezeichnen wir mit v(,, ¥) die Wahrscheinlichkeit, daB sich ein Elektron auf 
dem Wege x auf , Elektronen vermehrt hat und mit v(n,,,d—-) die Wahr- 
scheinlichkeit, daB sich », Elektronen auf der Strecke (d— x) auf m Elektronen 
vermehren. Die Wahrscheinlichkeit, daB eine Lawine von 1 Elektronen an der 
Stelle x aus n, Elektronen bestand, ist dann 


Waly, =) = esl 2m Ea) (4) 
>) v(m, +) u(ny, 2, d — +) 
n=1 


Fiir v(n,, x) legen wir die exakte Lawinenverteilung [s. Gl. (1)| zugrunde 


v(m, #) = e~** (1 —e7**)™—?. (5) 
Fiir die Darstellung von v(n,,, @— x) benutzen wir die Naherung fiir 7) >1 in 
Gl. (1). Dies ist wegen n;<n berechtigt. Die Lawine von » Elektronen kann als 
Summe von n, an der Stelle x startender Teillawinen aufgefaBbt werden, wobei jede 


einzelne der Verteilung 


u(n’,d—*#)=e 


rade) (6) 


—ald—4) , exp (—n’-e 
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unterworfen ist. Fiir die Verteilung der Summe n gilt dann? 
—2 (4-2) a4 
(mn, — 1)! 


Setzt man nun die Verteilungen (5) und (7) in die Gl. (4) ein, so erhalt man 


u(n,,n, d — #) OF wei Ai (ne -exp(—n-e *@-)). (7) 


, err ees 8 
W(n,, 4 = oan Sol (8) 

Dabei wurde zur Abkiirzung 
A = = (e** —1) (9) 


gesetzt. Fiir die Darstellung der Lawinenstruktur wahlen wir den 
Mittelwert von m,. Dieser ergibt sich zu 


7000 AP ee Pe ee Le kk 
Nn, (x) =A+1= = e*4 (1 a (10) 
a In Fig. 4 ist dieser mittlere Entwick- 
is lungsgang einer Lawine fir verschie- 
{ dene Werte von n/n, dargestellt. 
Bw Schreibt man die Gl. (10) in der 
rS Form 
nd ji,(x) —1 = ™ (&*—1), (10a) 
No 
3 


so erkennt man sofort die Naherungs- 


darstellung 
n 


aa . - x 
gike My (x) © ——-e*, (41) 
0 
Fig. 4. Der mittlere Entwicklungsgang einer 


Lawine fiir verschiedene Werte von m/m,. die dann giltig ist. wenn sowohl Ny (x) 
als auch e** groB gegen 1 ist. Die hier 


fiir das ebene Feld gegebene Darstellung iibertragt sich auf das in- 


homogene Feld, wenn man «- x durch [ «(&)-d& ersetzt. 
0 
Gemessene Lawinenverteilung. 

Eine Messung der Verteilung der LawinengréBe liegt vor von 
S. C. CurRAN, A. L. Cockrort und J. ANcus (7). Es wurde die Impuls- 
verteilung eines Proportionalzihlrohres bei Lichteinstrahlung gemessen. 
(Daten des Zahlrohres: Drahtdurchmesser etwa 0,15 mm, Gasfiillung 
250 Torr Argon + 250 Torr Methan, Spannung 3550 V.) Die maximale 
Feldinhomogenitait (Anderung der Feldstirke auf eine freie Elektronen- 
weglinge) betragt etwa 1%. Die mittlere Gasverstarkung wird zu 
1,46 - 10° el” angegeben. Bei einer Anderung der Spannung um 105 V 
wurde eine Anderung der Gasverstirkung um den Faktor 2 festgestellt. 


y Siehe Anhang 1. 
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Niramt man fiir « eine Feldstarkenabhingigkeit der Form «= A « e~8/F 
an, so kann aus den obigen Angaben der Wert von E/« naiherungsweise 
berechnet werden. Fiir die Umgebung des Drahtes, in der der wesent- 
liche Teil der Verstarkung stattfindet, ergibt sich E/x = 60 bis 100 V. 
Die Bedingung E/x > U, ist also nicht gut erfiillt. Fig. 5 zeigt die 
gemessene und die entsprechende ideale Verteilung. Die letztere 
wurde unter Zugrundelegung 00 
des experimentellen Mittel- $a 
wertes %=1,46-10° gezeich- = y59 
net. Wesentliche Abweichungen 
der beiden Kurven voneinander 
findet man nur fiir Lawinen, 
die kleiner als 1/3 sind. Ob 
diese Abweichungen auch auf- 
treten, wenn die Voraussetzun- 
gen 1 und 2 der theoretischen 
Verteilung besser erfiillt sind, 


mu8 noch experimentell ge- Fig. 5. Vergleich deridealen Lawinenverteilung (- - - - - - ) 
klairt werden. mit Messungen am Proportionalzahlrohr. MeBwerte 
nach CuRRAN u.a. [7]. 


€ 
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III. Verteilung der Nachlieferungselektronen. 


Fiir diese Rechnungen werden nur solche Nachlieferungsprozesse 
betrachtet, bei denen die Elektronen an der Kathode erzeugt werden 
(gasionisierende Strahlung z.B. ist also ausgeschlossen). Ob nun die 
Nachlieferung durch positive Ionen, metastabile Atome bzw. Molekiile 
oder durch Photoeffekt an der Kathode erfolgt, stets handelt es sich um 
eine groBe Anzahl von Elementarprozessen, wobei jeder mit der (sehr 
kleinen) Wahrscheinlichkeit 6 ein Elektron an der Kathode auslésen 
kann. Die Anzahl # der Elementarprozesse sei der Elektronenzahl n 
der Lawine proportional. Die Wahrscheinlichkeit, daB durch eine 
Lawine k Nachlieferungselektronen ausgelést werden, ist also 


w(n, k) = (7) #0 — 6)P-*, (12) 

Da p>1 und 6<1, ist (12) in guter Naherung eine Potsson-Verteilung 
. = (5:p)* 9-5 

w(n, k) = Pe 8P. (13) 


Wegen der Proportionalitat zwischen » und # gilt 


w(n, k) = Ys a pe (44) 
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Dabei ist y die Anzahl der im Mittel pro Lawinenelektron ausgeldsten 
Nachlieferungselektronen. Diese Verteilung gilt sowohl fiir Rick- 
wirkung nur einer Art (z.B. nur durch Ionen), als auch fir gemischte 
Riickwirkung (z.B. Ionen und Photonen). Dann ist in Formel (14) ein- 
zusetzen: ¥=7tonen + YPhotonen: Dies ist sofort durch Berechnung von 


k 
w (1, k) 3 bs Wionen (1, v) * Wphotonen (n, ea v) 
v=0 
zu verifizieren. 
Ist nun die statistische Verteilung der Lawinengr6Be n durch Angabe 
der Wahrscheinlichkeiten v() gegeben, so ist die Wahrscheinlichkeit @,, 
daB eine belichige Lawine k 


10 
Nachheferungselektronen aus- 

08 . e 
: lést, zu ermitteln gemaB 
96 = 
a4 w, = > v(n)-w(n, k) 
. pr (15) 
; ~ fv(n)-w(n,k)-dn. 
“Toe a5. 1. 2 apes 

as Fiir die ideale Verteilung (1) 


Fig. 6. Verteilung der Nachlieferungselektronen einer La- hal : : d 2 
wine in Abhangigkeit von “= yrend, Die Wahrscheinlich- erhalt man mit (14) und mut 


keiten w, fiir k Nachlieferungselektronen sind ubereinander der Abkiirzung uN y die 

aa auch in Fig.6 dargestellte 
Verteilung der Nachlieferungselektronen einer Lawine, wie sie bereits 
fiir Nachlieferung durch positive Ionen angegeben wurde [3). 


ing osha 6 
hae cae Oe (1 ) 
(1 + 4) 

Wy hat bei 4 =0 den Wert 1 und geht mit wachsendem ~ monoton gegen 
Null. Die tbrigen w, haben ein flaches Maximum bei uw =f. Fiir groBe u 
ist die Verteilung sehr breit, denn @,., ist nur sehr wenig von @, ver- 
schieden. 

Wie man leicht zeigt, ist R=y, d.h. eine Lawine lést im Mittel u 
Nachlieferungselektronen aus. 


IV. Lawinenketten. 


Kommt die mittlere Anzahl der Nachlieferungselektronen einer 
Lawine in die GréBenordnung von 1, so treten mit merklicher Wahr- 
scheinlichkeit Lawinenketten auf. Dabei lésen die primar startenden 
Elektronenlawinen eine Schar von Nachlieferungselektronen aus, welche 
ihrerseits wieder Nachfolger haben usw. Es bildet sich eine Lawinen- 
kette, die aus einer ganzen Reihe von ,,Generationen‘’ bestehen kann. 
Ist «> 1, so wird die Kette mit einer von Null verschiedenen Wahr- 
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scheinlichkeit nicht abreiBen und somit zum Durchschlag der Entladungs- 
strecke fithren. Nach der Anzahl N, der Primiarelektronen unterscheiden 
wir Makroketten (Nj >1) und Mikroketten. Bei letzteren werden wir 
speziell den Fall Nj=1 betrachten. Die Entwicklung der Makroketten 
_ bildet kein besonderes statistisches Problem, da hier sehr gut mit Mittel- 
werten gerechnet werden kann. Die Lawinenzahl der n-ten Generation 
ist dann einfach Ng-y"~'. Derartige Ketten wurden yon LAvER [8] 
beobachtet. Eine Frage, die sowohl Makroketten als auch Mikroketten 
betrifft, ist die zeitliche Folge der Generationen und deren Streuung. 


A. Die zeitliche Folge der Generationen. 

Zur Zeit f=0 mégen N, Elektronen von der Kathode der ebenen 
Entladungsstrecke starten. Mit welcher Wahrscheinlichkeit dW wird 
in dem Zeitraum ¢...¢+d¢ ein Nachlieferungselektron an der Kathode 
ausgelést ? 

a) Nachlieferung durch Photoeffekt. Der Einfachheit halber nehmen 
wir an, da8 nur ein Energiezustand durch ElektronenstoB angeregt 
werden kann. Die mittlere Lebensdauer des Zustandes sei t,. Werden 
in der Zeitspanne ?’...¢’ +d?’ in der Entladungsstrecke  (t’) - dt’ Molekiile 
angeregt, so werden von diesen Molekiilen in der Zeit ¢...¢+dt 

n(t') dt’ -e—¢-tin Ut 
TY 
ein Photon emittieren und in den Grundzustand zurtickkehren. Ins- 


gesamt werden also im Intervall dt 
t 
wT eit [ n(t') edt’ (18) 
Ty Fs 
Photonen emittiert. Dieser Photonenzahl ist die gesuchte Wahrschein- 
lichkeit dW proportional, wenn die Absorption im Gas zu vernachlassigen 
ist. Da die Zahl der anregenden StdBe m (¢’) - dt’ der Anzahl der zur Zeit ¢’ 
laufenden Elektronen — das sind N,- e**°"", wo t_ die Elektronen- 
laufzeit ist — proportional ist, gilt 
Min (t, 7) 
dW~edt- f [exp (ad/r_ + 1/r)¢'] dt’. (19) 


0 


(17) 


Insgesamt werden N, - ,, Nachlieferungselektronen ausgelést. Die auf 1 
normierte Wahrscheinlichkeitsdichte /(¢) wird definiert durch 


ay = ft), Hp" dt (20) 
und es ist wegen (19) 
peice fete Wr, — ef) tice oh 4a) 7. 
iT (e%4 =>) 
f()=4 ° : (21) 
/ 2a (cx4 Aga? a) e-t—H)im fiir t>t_. 
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Dabei wurde t_/xd+1,—T, gesetzt. Da im allgemeinen gt >41 ists 
kann (21) noch vereinfacht werden: 


t-t 
1 2 tied fr Off = '5_ 
f(i) = 4 ree (22) 
ey Wi. FS hae 


To 
Diese Verteilung ist in Fig. 7 schematisch dargestellt. Der Mittelwert t 
der Verteilung liegt bei 
t=t_(1 —1/ad) + 17,. (23) 
Dies ist der mittlere zeitliche Abstand zweier aufeinanderfolgender 
Lawinengenerationen. Er kann unter Umstanden erheblich von der 
fy) Elektronenlaufzeit t_ abweichen, da 
= T, etwa 10S sec betragt und damit 


1) , ; - 
besonders bei kleinem d mit T_ ver- 
gleichbar ist. 

T, - - . 
eel Fiir die mittlere Streuungo der 
0 r 2 7 Verteilung findet man 
Ti te ee 
ad *7 
oe 2 A 
Fig. 7. Zeitliche Verteilung der Wahrscheinlichkeit a= (T_ /xd) r Ty- (24) 
fiir die Ausl6sung von Nachlieferungselektronen Ss ’ 
durch Photoeffekt. Fir die Streuung o, der n-ten 


Generation gilt o,=n-o%. Diese 
wachsende Streuung bedeutet, daB die Lawinen einer Generation zeitlich 
immer mehr auseinanderlaufen. Wird diese mittlere Streuung gréBer als 
t, so kénnen die einzelnen Generationen nicht mehr voneinander ge- 
trennt werden. Die Anzahl der etwa noch zu trennenden Generationen ist 

7 (ad —1+ f)® : Ty 
.. -— —— t =_ 5 
a 1+ # “e t_jad (25) 
Ist P< 1, so gilt m= (ad —1)*, d.h. m kann von der GréBenordnung 100 
werden. Bei B=1 ist immerhin noch n= 4 (xd)* um erst bei B=ad 
auf etwa n= 4 abzusinken. 

b) Nachlieferung durch Ionen. Die Wahrscheinlichkeit fiir die Aus- 
lésung eines Nachlieferungselektrons ist proportional zu der jeweils an 
der Kathode ankommenden Ionenzahl. Unter Vernachlassigung der 
Diffusion kommen in der Zeit dé 

ap t's ; <— dt 
iIN=ad:Nve %+e ._< 
l od ne a (26) 
Tonen an. Dabei ist t, die Ionenlaufzeit. Mit dW = f(t) - Nott, di~dN 
folgt fiir die auf 1 normierte Wahrscheinlichkeitsdichte f(t) 


f(t) = ‘ie. fOr: Oatictetelees (27) 


0 tir £> ‘ty -+- 4... 
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Diese Verteilung ist der fiir Photonennachlieferung gefundenen dhnlich. 
Es ist dort nur t_ und t, durch t_ + t, bzw. 0 zu ersetzen. Fiir die ent- 
sprechenden Mittelwerte erhalt man also 


t= (t, + T_)(1— 1/ad) wt, (1 — 1/ad) 
tT, + Tt Ty 


Co = ] 


ad ad ° 


Die Anzahl der noch zu trennenden Generationen wird = (ad —1)?. 
Aus diesen Rechnungen folgt, daB bei einem Gasentladungsaufbau 
liber viele Generationen hinweg (TowNSEND-Aufbau oder Mehrlawinen- 
Kanalaufbau [9]) die Struktur dieser Generationen zumindest itiber 
10 bis 100 Generationen merklich bleibt und dann langsam verschmiert 
wird. 
B. Mikroketten. 

Wie bereits erwahnt sollen hier nur Mikroketten mit N,=1 be- 
trachtet werden, d.h. es wird immer nur ein Elektron durch duBere 
Einwirkung aus der Kathode ausgelést. Wegen des kleinen N, kann 
jetzt nicht mehr mit dem Mittelwert Nj ~*~} fiir die Lawinenzahl in 
der k-ten Generation gerechnet werden. Im folgenden Abschnitt soll 
untersucht werden, wie lange nach dem Start dieses Elektrons noch das 
Spiel der Elektronenlawinen in der Entladungsstrecke lauft. Daraus 
ergibt sich zwangslaufig eine Angabe iiber die Wahrscheinlichkeit, mit 
der ein Elektron zur Ziindung der Entladungsstrecke fiihrt. Anschlie- 
Bend wenden wir uns der Struktur der Mikroketten zu. Dabei wird die 
Frage behandelt, wie eine Kette von z.B. 10 Generationen Dauer ab- 
lauft, wieviel Lawinen also in jeder einzelnen Generation zu erwarten 
sind. 

Die Verteilung der Generationenzahl von Lawinenketten. 

Es werde mit K,, die Wahrscheinlichkeit bezeichnet, daB eine Kette 
aus mindestens n Generationen besteht. Aus genau n Generationen be- 
steht eine Kette dann mit der Wahrscheinlichkeit K,—K,.,. Die 
Numerierung der Generationen beginnt mit 1 bei der von anBen ausge- 
lésten Lawine. Mit der Wahrscheinlichkeit w, treten dann in der zweiten 
Generation k Folgelawinen auf [vgl. Gl. (15) und (16)]. Soll die Ge- 
samtkette mindestens »-+1 Generationen umfassen, so muB jedenfalls 
eine der k Nachfolgelawinen Ausgangspunkt einer Kette von mindestens 
n Generationen sein. Mit der Wahrscheinlichkeit (1 —,,) bricht aber 
eine Kette vor der n-ten Generation ab. Also wird mit der Wahr- 
scheinlichkeit 1 — (1 —K,,)* mindestens eine Kette der # Nachfolger die 
n-te Generation erreichen. Die Gesamtkette umfaBt dann mindestens 
n +41 Generationen. Fir K,, gilt somit die Rekursionsformel 


Pere) =So, (4 — (14 — K,)*). (28) 


Zeitschrift fiir Physik. Bd. 140. 16 
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Entnimmt man @, der Gl. (16), so erhalt man 
1 i 
Kyat 7 Ky” 
Da jede Kette mindestens aus einer Generation besteht, ist Ay =1 und 
man erhalt 


Ke ee fiir ent (29) 


K,, =1/n fir pw=1. 


Die Fig. 8 zeigt K,, als Funktion von yp. _ 
Bildet man K,,= lim K,,, so findet man die ,,Ziindwahrscheinlich- 

keit‘‘? ae 
Kg = 0 fir wai 


(30) 
K,=1—1/p fir p> 1. 


Es besteht also fiir 44> 1 eine von Null verschiedene Wahrscheinlichkeit, 
daB ein einzelnes Elektron eine Lawinenkette auslést, die nicht wieder 
abreiBt. Damit ist die Ent- 
ladung selbstandig geworden. 
Wie hieraus deutlich hervor- 
geht, bedeutet die ,,Ziindbe- 
dingung «> 1 nicht, daB auch 
jedes aus der Kathode ausgeléste 
Elektron zur Ziindung fiihrt. 


OF a we whee . . ” Die Struktur der Lawinenketten. 


Fig. 8. Die Wahrscheinlichkeit Ky, da8 ein Elektron eine Wie Gl. (29) zeigt, treten 
Lawinenkette von mindestens n Generationen ausliést. langere Lawinenketten mit be- 
achtlicher Wahrscheinlichkeit 
auf, sobald uw in die Nahe des Wertes 1 kommt. So fihren z. B. bei 
/4=0,9 bereits 6% aller Primarelektronen zu Ketten von iiber 10 Gene- 
rationen Dauer. Unter der Struktur der Lawinenketten soll nun 
folgendes verstanden werden: Wir betrachten bei festem yz alle Lawinen- 
ketten von genau 7 Generationen Dauer und fragen, wieviel Lawinen 
werden im Mittel in der k-ten Generation gleichzeitig laufen. 
An den Anfang unserer Rechnung stellen wir die Angabe der Wahr- 


scheinlichkeit W(k, m), daB iiberhaupt in der k-ten Generation einer 
Kette m Lawinen laufen. 


k—1 
W(h, 0) = 2 
ane (31) 
W(k, m) = 1.2 at. a ee ee. 0 
Pee Ai f we—1) Poa ro 


1 Diese Ziindwahrscheinlichkeit wird bei W1jSMAN [3] auf etwas anderem Wege 
abgeleitet. 
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Fiir ~=1 gehen diese Formeln iiber in 
W(k, 0) =1 —4/k 


m—1 
Wk, m) = @=1) tar mo. | (2) 


pmrd ‘. 


Der Beweis der Formeln (31) wird durch volistandige Induktion gefiihrt. In 
der (& — 1)-ten Generation laufen mit der Wahrscheinlichkeit W(k —1, v) »y Lawinen, 
Diese erzeugen mit der Wahrscheinlichkeit U(v, m) genau m Folgelawinen in der 
k-ten Generation der Kette. U(v, m) wird aus Gl. (16) berechnet und zwar ist! 


U(0,0)=1, Ul, m) =0 eer 
U(v, m) = ‘ 7; ai ™) aes (vy + 0). (33) 


Fiir W(a, m) erhalt man durch Summation iiber alle » die Rekursionsformel 
co 
W(k, m) = } W(k — 1, v) - U(v, m). 
v=0 


Da die Lawinenketten stets mit einer Primarlawine beginnen sollen, ist W(1, 1) =1 
und W(1, m) =0 fiir m +1. 

Um die Wahrscheinlichkeit W(r, k, m) zu bestimmen, daB in einer 
y-gliedrigen Kette in der k-ten Generation m Lawinen laufen, ist noch 
zu ermitteln, ob die in &-ter Generation auftretenden m Lawinen auch 
zu einer r-gliedrigen Kette gehéren. Dazu ist erforderlich, daB sich die 
m Lawinen noch weiter iiber genau r—k Generationen fortpflanzen. 
Dies geschieht in der Weise, daB mindestens einer der m Lawinen eine 
Kette von genau s=7-+ 1—£# Generationen bildet, wahrend die anderen 
Lawinen zu héchstens s-gliedrigen Ketten fiihren. Die Wahrscheinlich- 
keit dafiir ergibt sich zu 


w (m, s) = (1 — K,,,)” — (1— K,)”. G4) 


Es ist namlich (1— K,) die Wahrscheinlichkeit, da8 eine von einer Lawine 
ausgehende Kette die s-te Generation nicht erreicht. Damit m Lawinen sich genau 
bis zur s-ten Generation fortpflanzen, darf keine der m Einzelketten die (s+ 1)-te 
Generation erreichen ((1—K,,,)”). Davon abzuziehen sind die Falle, wo alle 
m Einzelketten nicht bis in die s-te Generation gelangen ((1 — K,)”). 


Mit der in (29) angegebenen Darstellung fiir K, folgt 
rw ps —1 a pet = 1 
w(m,s) = (Fs) ( Pos ) (4 = 1) 
s2m — (s? — 1)™ 


eas leat 


(35) 


w(m,s) = 


Nach den obigen Ausfiithrungen folgt fiir die gesuchte Wahrscheinlichkeit 


Wr, k, mi) =m) elmer pine (36) 


CO 
x Wk, m) + w(m,r +1—f) 


m=1 


1 Siehe Anhang 2. 
Ge 
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Um eine Aussage tiber den mittleren Verlauf einer r-gliedrigen Kette 
zu erhalten ware dann gemaB 

in(k, r) =>. m- Wr, k, m) (37) 

m=1 

der Mittelwert von m zu bilden. Setzt man die allgemein gultigen Dar- 
stellungen von W(k, m) und w(m, r+1—k) [s. Gl. (34) und (35)] in (36) 
bzw. (37) ein, so erhalt man zundchst recht uniibersichtliche Formeln. 
Da die Bildung endlicher Lawinenketten praktisch nur in der Umgebung 
von 4=1 gut zu beobachten ist, bedeutet es keine wesentliche Ein- 
schrinkung, wenn man sich auf diese Umgebung beschrankt. Setzt man 


| 
08 ALITA 
| 4% Tp > 
M 94 ef 
ah. m2 
rr a ae er toi. 0 dae 
i—- j ae 


Fig. 9 au. b. Darstellung des Ablaufs einer Lawinenkette von r = 10 Generationen a) durch die Wahrscheinlich- 
keiten W(r, k, m) fiir das Auftreten von m Lawinen in der &-ten Generation; b) durch die mittlere 
Lawinenzahl m(&,r) in der &-ten Generation. 


@=1-+ e und entwickelt W(r, &, m) nach Potenzen von e, so verschwin- 
det das in ¢ lineare Glied. Dies bedeutet insbesondere, daB m(h, 7) fiir 
j#@=1 ein Extremum annimmt (es handelt sich um ein Maximum) und 
folglich in der Umgebung von «=1 nur geringe Anderung zeigt. Fiir 
=1 findet man 

: r+ti)-r k—1\"(jr+i—k\m / y—k \m 

W(r, k, m) = aan (4 (4) =(=4 7} 68) 

und 

m(k, 7) = (2k —1-+ k/r) (1 — Rj +4). (39) 
Bei festem 7 stellt m(k,r) eine Parabel dar. Bei nicht zu kleinem 2 
liegt das Maximum etwa bei k=7/2+-1 und nimmt rund den Wert 7/7 
an. In Fig. 9 sind als Beispiel fiir eine zehngliedrige Lawinenkette die 
Wahrscheinlichkeiten W(10, &, m) tiber k und m aufgetragen. AuBerdem 
ist die Kurve m(k, 10) eingezeichnet. Dieses Schaubild zeigt deutlich 
die starke Streuung von m um den Mittelwert m. In der letzten Gene- 
ration (= 10) findet man eine starke Haufung bei m—=1 und 2, denn 
die gréBeren Lawinenzahlen haben nur eine sehr geringe Wahrschein- 
lichkeit, keine neuen Elektronen auszulésen. 


"| 
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Beobachtungen des Verfassers! an einer unselbstindigen Entladung 
in unmittelbarer Nahe der Ziindspannung kénnen als Mikroketten der 
hier beschriebenen Art gedeutet werden. Dabei handelt es sich um Ketten 

von etwa 100 Generationen Dauer. 


V. Die ziindende Townsenp-Entladung. 


Es wurde gezeigt, daB im Falle “> 1 ein einzelnes Elektron mit der 
Wahrscheinlichkeit 1—1/u die Ziindung der Entladungsstrecke herbei- 
fihrt. Betrachten wir nur solche Ketten, die zur Ziindung fiihren, so 
k6nnen wir wieder bestimmen, mit welcher Wahrscheinlichkeit Z (k, m) 
in der &-ten Generation m Lawinen auftreten. Damit die Entladung 
zundet, mu8 mindestens eine dieser Lawinen Ausgangspunkt einer un- 
endlich langen Kette sein. Dies ist mit der Wahrscheinlichkeit z(m) 
der Fall. 

z(m) =1— (1/u)”. (40) 


Beweis. Eine Lawine fiihrt mit der Wahrscheinlichkeit 1— Koo = 1/ nicht zur 
Zindung [s. Gl. (30)]. Die Wahrscheinlichkeit, daB von m Lawinen keine ziindet 
ist (1/u)”. Mit der Wahrscheinlichkeit 1— (1/u)" wird also mindestens eine der 
m Lawinen die Ziindung herbeifiihren. 


Analog wie bei den endlichen Ketten wird nun nach der Formel 


Z(k, m) ea WUE oa) = z(m) _ (41) 


berechnet. Setzt man hierin (31) und (40) ein, so erhalt man als End- 
ergebnis 


Zhan) = EB) (yp 4) (H Lanett L) . (42) 


Der mittlere Verlauf der ziindenden Entladung wird durch den Mittel- 
wert von m als Funktion von k beschrieben: 


eT At IP (43) 


m=1 


Die Auswertung mit Hilfe von (42) ergibt 


£0—— (44) 


Fiir groBe & (genauer fiir y?*~*>>1) kann diese Formel noch etwas ver- 
einfacht werden: 
m(k) ~ wu" (yw — 1). (45) 


1 Wird in Kiirze veréffentlicht. 


+ 


236 WERNER LEGLER: 


Ohne Beriicksichtigung der hier dargelegten Verteilungen rechnet man 
in der k-ten Generation einer von einem Elektron ausgehenden Entladung 
mit im Mittel *~! Lawinen. Der tatsidchliche Mittelwert ist nach (45) 
um etwa den Faktor yu/(—1) groBer. Dieser Faktor ist aber nur dann 
wesentlich von Eins verschieden, wenn yz nur wenig groBer als Eins ist, 
wie es bei sehr kleinen Uberspannungen der Fall ist. Unter diesen Um- 
stinden ist auch ein Einflu8 auf die Berechnung der Aufbauzeit der 
TowNnsEND-Entladung méglich. In Fig.10 ist der mittlere Verlauf der 
Lawinenzahl pro Generation einer durch ein Elektron gestarteten 
TownseEnD-Entladung fiir ~=1,1 nach 
Gl. (44) dargestellt. AuBerdem ist noch 
der Verlauf fiir ~—-1 eingezeichnet. Der 
Vergleich mit der Kurve p*~? fiiru=1,1 
zeigt den groBen EinfluB der statisti- 
schen Verteilung auf das Anwachsen der 
Entladung. Eine experimentelle Priifung 
dieser Zusammenhange ist relativ einfach 
durchfiihrbar, indem man den Strom 
der ziindenden Entladung beim stati- 
schen Durchbruch oszillographiert. Da- 
mit nur ein einzelnes Elektron die 


Fig. 10. Mittlerer Verlauf m(&) einer durch a ne 
ein einzelnes Elektron geziindeten Town- Ziindung herbeifihrt, ist nur extrem 


seND-Entladung. Zum Vergleich die Kurve — a 
sr deeb Gln toamaeitaig es me schwache Bestrahlung der Entladungs 


Statistik ergeben wiirde. strecke anzuwenden. 


Die Verteilung der Aufbauzeit. 

Als Aufbauzeit der Entladung bezeichnen wir die Zeit, die vom Start 
der Primarlawine vergeht, bis in einer Generation die Lawinenzahl m 
einen kritischen Wert N tiberschreitet. Bei diesem kritischen Wert NV 
wird die Entladung auf Grund der Raumladungswirkung entscheidend 
umgebildet. N liegt je nach den Entladungsbedingungen in der GréBen- 
ordnung 10° bis 10” [10}. Die Aufbauzeit wird hier in Vielfachen des 
zeitlichen Abstandes zweier Generationen ausgedriickt. Die Wahr- 


scheinlichkeit, daB die Ziindung bis zur k-ten Generation einschlieBlich 


erfolgt, wird durch ps 


W, (8) = Z (bm) (46) 

m= N 
gegeben. Diese Formel erfaBt von allen Ketten, die einmal mehr als 
N Lawinen in einer Generation fiihren nur solche, die auch ohne die 
ziindungsfordernde Wirkung der Raumladung unendlich lang wiirden. 
Die hierbei vernachlissigten Falle sind aber extrem unwahrscheinlich!. 


1 Die Wabrscheinlichkeit ist kleiner als u~%, bei N = 108, 7 = 1,001 also kleiner 
als 107454, 
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Setzt man nun Z(, m) nach Gl. (42) ein1, so erhalt man 


W.( )= ers | 1 — t/a" ; . ai 


_ Da W,(k) erst dort wesentlich von Null verschieden wird, wo yu* etwa 


die GréBenordnung von N erreicht, also groB gegen Eins wird, kénnen 
wir die Naherung (1—14/u*) = e~¥“ benutzen und erhalten 


W,(#) = exp (— z-*-N). (48) 


Fir eine iibersichtlichere Darstellung 
cer Verteilungen ist es zweckmabBig, 
in der Skala der Generationenzahl 
einen Bezugspunkt ky von der un- 
gefahren GréBe der Aufbauzeit zu 
definieren : 


In N- (4 —1 y 
ky = — 42 (49) 70 8 90 00 110 120 130 
J 4199 — 
Setzt man dies ky in (48) ein, so Fig. 11. 


Die Haufigkeitsverteilung der Aufbauzeit einer 
TownseEnp-Entladung. Die Skala ,,k/k- 100“ 
gilt fiir den Fall, daB In(N-(u—1)) =18 ist. 


erhalt man 
W,(k) = exp (—u-*-). (50) 


Die Verteilung der Aufbauzeit wird also unabhangig von den speziellen 


_ Werten von WN und w, wenn man (k—,) - Inw als unabhangige Variable 


wahlt. Die zugehérige Haufigkeitsverteilung (differentielle Verteilung) 


_ der Aufbauzeiten (s. Fig.11) lautet: 


é 


_ gung der Statistik, sosetzt man einfach w= N und erhalt ky =— 


: 


d 


w, (#) = W, (2) = In p~ exp (— pF), (54) 


Berechnet man den Mittelwert k der Aufbauzeit gem&B der Definition 


oes {k-dW,(k), so findet man die Beziehung 
0 


Dabei ist C die EuLersche Konstante. Die mittlere Aufbauzeit der 
TowNnsEND-Entladung aus einem einzelnen Elektron ist also durch 


ea In i Ha A (53) 
ny 
gegeben. Berechnet man dagegen die Aufbauzeit k, ohne Beriicksichti- 
InN: 
ip e- 


1 Dabei wird in dem Faktor (w”—1) die 1 gestrichen, da in der Summe (46) 
stets yw" >pN>1 ist. So ist z.B. in dem ungiinstigsten Fall N= 10°, w= 1,001 


me = 1084, 


4 


4 
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Eine wesentliche Abweichung von dem Mittelwert k ist nur in der Nahe 
von 4=1 zu erwarten. Nimmt man N =108 an, so ist bei pe=1,1 bzw. 
w=1,01 der Wert von k, etwa 11 bzw. 30% groBer als k. Da aber die | 
bisher gemessenen und berechneten Aufbauzeiten um etwa einen Fak- 
tor 2 differierten, kommt dieser Abweichung noch keine groBe Bedeu- 
tung zu. 

AbschlieBend soll noch die mittlere quadratische Streuung der Auf- 
bauzeit bestimmt werden. Diese ist definiert durch 


of =f (kk)? dW, =f RdW, -P=R_—P. (54) 
0 


0 
Fiir den Mittelwert von #2 erhalt man 


2 = kh + 2C k,/ln wu + B/(Iny)?. (559 


Dabei ist B eine Konstante, die durch 
p= f(n x)?-e-*dx 1,98 
gegeben ist. Nach (52) ist : 
R= k2 + 2C k,/ln w + C2/(In w)?. 
Somit finden wir fiir die Streuung o die Beziehung 


_ B-@ 1,28 
ee Ing Dy oe (56) 


Die relative Streuung der Aufbauzeit o/k betragt dann etwa 5 bis 8% 
und ist ziemlich unabhangig von N und ys. 


Anhang I. 


Im Abschnitt II wurde gezeigt, daB fiir »)>1 die LawinengréBe » der Vertei- 
lung 
1 ; 
1(n, d) = -e@ 1% — x4 5 
uv (n, d) =" e ‘ Nm =e (57) 
unterworfen ist. Welche Verteilung v(m, N, d) gilt nun fiir die Gesamtzahl N der 
Elektronen von m Lawinen ? 
Fiir die Verteilung (+) der Summe von zwei unabhangigen statistischen Va- 
riablen mit den Verteilungen /,(*) und f,(*) gilt nach CRAMER [12] 


+00 
f(a)= ff fla) tele — 4) day. (58) 
4,=—0 
Mit Hilfe dieser Formel findet man z.B. fiir zwei Lawinen 
N 


v(2, N,d) = [2 Lent, 1. (N= ay 
No No 
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Durch vollstandige Induktion ergibt sich 


iN silt ai 
v(m, N,d) =. E) ge DT (59) 


Ny “(m — 1) ! 7 Noy 
Ersetzt man mdurch»,, N durch x und d durch d—x, so erhalt man die in Abschnitt IT 
angegebene Formel (7) 
—a(d—x) (”° e—%(d—x))m—1 


n,,n,a — = 
vim asia (n, — 1)! 


-exp (— n+e—%(d—2)), 


Anhang 2. 


Lést eine Lawine im Mittel « Nachlieferungselektronen aus, so war die Wahr- 
scheinlichkeit, daB & Elektronen ausgelést werden, gegeben durch 


pe 
=- —_“_. 16 
Op a1 + pet - ( ) 


Laufen nun statt einer Lawine y Lawinen in der Entladungsstrecke, so werde mit 
U(v, m) die Wahrscheinlichkeit bezeichnet, daB insgesamt m Nachlieferungselek- 
tronen erzeugt werden. Da sich die Lawinen gegenseitig nicht beeinflussen, ist 


U(v, m) = ¥ U(1, ky) - U(A, Ry) -...- UA, By) (60) 


v 
Dabei ist iiber alle Produkte zu summieren, fiir die > k;=m gilt. Fihrt man durch 
$=1 


U(v, t) = > U(v, m) - (61) 


die erzeugende Funktion ll(y,#) der Verteilung U(y, m) ein, so folgt aus (60) 


U(v, t) = (W(1, 8)”. (62) 
Nach (16) ist 
> ; pun 
ON. 8) = tn = yee 
und damit 
1 1 
u 4-5 = ~-———_—— « —_-—___— , 63) 
(1, t) i ew ia ( 
1 ft 


Fiir U(y, t) ergibt sich nach (62) und (63) 


1 


Le = (1 eh a ae ) 


3 


Der Koeffizientenvergleich mit (61) ergibt 


U(v, m) = (—1)"(— i f 


\ Mm 
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a ; ; ‘— t, = 


. [ J _ : a 7 2. y i A . 7 q 


oe , i (~ = en 
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2 / a 

i= 

folgt die bei (33) angegebene Darstellung 

y—1+ ™) | 


0m = ("ae ina 
fiir (v +0). Die dort noch angegebenen Werte 
CHG, Op te U(0, m) = 0 (m += 0) ex 


folgen als Trivialitaéten aus der Definition von U(y, m). 


Diese Arbeit wurde als Diplomarbeit am Institut fiir An r 
Physik der Universitat Hamburg angefertigt. Dem Leiter des Institutes 
Herrn Professor Dr. H. RAETHER, danke ich fiir das standige Interesse 
an dieser Arbeit und fiir anregende Diskussionen. % 
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